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Summary 



Globally-propagating bright fronts in the solar corona are a poorly 
understood phenomena despite almost 15 years of research. Here, the 
kinematics and morphology of these disturbances are studied using 
Extreme Ultra Violet (EUV) observations from the Solar TErrestrial 
RElations Observatory (STEREO) and Solar Dynamic Observatory 
(SDO) spacecraft. 

The first observations of a coronal bright front (CBF) made using 
STEREO are presented, with the pulse observed in all four EUV 
passbands, namely 171, 195, 284, and for the first time, 304 A. The 
pulse exhibited similar kinematics in all passbands studied, with a 
velocity peak of 238±20 km s _1 within ~28 minutes of launch and 
acceleration varying from 76 m s~ 2 to -102 m s~ 2 observed in the 
304 A passband. Similar kinematics are found in the 195 A passband, 
while lower values were found in the lower cadence 284 A passband. 
The higher cadence 171 A passband shows much stronger variation, 
with the velocity peaking at 475±47 km s -1 within ~20 minutes of 
launch, and the acceleration varying from 816 m s~ 2 to -413 m s~ 2 . 
The derived pulse velocity and acceleration were found to be strongly 
influenced by the observing cadence, implying that previous kinemat- 
ics may have been underestimated. 

Several different image processing techniques for identifying CBFs 
were tested to determine the optimal method for CBF pulse analysis. 
It was found that the techniques traditionally used to identify CBFs 
are fundamentally flawed and prone to undefined user-dependent er- 
rors. The kinematics of a simulated data-set were derived using 
several different numerical differencing techniques. These techniques 



were found to be inadequate for dealing with the typical data-set sizes 
encountered when studying CBFs. A semi-automated identification 
and tracking algorithm is proposed to provide a statistically rigorous 
analysis of CBFs. 

This algorithm was applied to STEREO observations of four CBF 
events and used to determine their kinematics. Following launch at 
initial velocities of ~240-450 km s -1 each of the four events stud- 
ied showed significant negative acceleration ranging from ~ —290 to 
—60 m s~ 2 . The spatial and temporal widths of the CBF pulses were 
also found to increase from ~50 Mm to ~200 Mm and ~100 s to 
~1500 s respectively, indicating that the pulses are dispersive in na- 
ture. The variation in position-angle averaged pulse-integrated inten- 
sity with propagation was examined to identify evidence of dissipation, 
but exhibited no clear trend across the four events studied. 

A CBF pulse was also observed using contemporaneous data from 
both STEREO and SDO. The CBF exhibited a lower initial veloc- 
ity with weaker deceleration in STEREO observations (~340 km s _1 
and —72 m s~ 2 ) compared to SDO (~410 km s" 1 and —279 m s~ 2 ). 
The CBF kinematics from SDO were found to be highly passband- 
dependent, with an initial velocity ranging from 379 ±12 km s _1 to 
460 ± 28 km s" 1 and acceleration ranging from —128 ± 28 m s~ 2 
to —431 ± 86 m s -2 in the 335 A and 304 A passbands respec- 
tively. These kinematics were used to estimate a quiet coronal mag- 
netic field strength range of ~l-2 G. Significant pulse broadening was 
also observed, with expansion rates of ~130 km s _1 (STEREO) and 
~200 km s _1 (SDO). By treating the CBF as a linear superposition of 
sinusoidal waves within a Gaussian envelope, the resulting dispersion 
rate of the pulse was found to be ~8-13 Mm 2 s -1 . 

All of these results indicate that coronal bright fronts are best inter- 
preted as fast-mode magnetoacoustic waves propagating in an inho- 
mogeneous medium. 
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Chapter 1 
Introduction 



Here comes the Sun 
And I say it 's alright 

George Harrison 



No physical object (aside from the Earth) has had more influence on human 
civilization than the Sun, with the heat and light that it provides allowing life to 
flourish on Earth. Every individual culture has worshipped the Sun in its own 
way, telling epic tales of the Sun Gods (such as Helios, Ra and Surya; the Greek, 
Egyptian and Hindu Sun Gods respectively), building temples that align with 
the rising or setting Sun (e.g., Newgrange) and by treating its heroes, kings and 
saints as human incarnations of the divine Sun (cf. religious halos and "Le Roi 
Soldi"). 

While science and technology have removed much of the mystique that has 
traditionally surrounded the Sun, humanity has not stopped its worship of our 
closest star. Every year millions of people travel to warmer parts of the globe to 
bask in the healing powers of the Sun, while the suntan has become symbolic of 
health and vigour. Although the ability of the Sun to burn is a cause of some 
caution and alarm, humanity still finds it impossible to live without the heat and 
light that only the Sun can provide. 
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The Sun may have been present throughout humanity's stay on the Earth, 
but it is only recently that it has become a source of endless investigation and 
countless unanswered questions. The first inkling that the Sun may have an effect 
on Earth that extends beyond the provision of light and heat came from obser- 
vations of a flash on the Sun made by Carrington (1859). Carrington observed a 
flash of white light in a particularly large sunspot region when making his daily 
observation of the Sun. It was also noted that the magnetic field experiment at 
Kew Observatory in London exhibited a strong reaction to some unknown driver 
at about the same time. Several days later, electrical and magnetic equipment 
all around the world began behaving erratically with no apparent cause, with a 
telegraph operator on the transatlantic telegraph from Valentia to Newfoundland 
thrown across the room by a large electric shock coming from the wires while 
aurorae were observed as far south as Cuba. Although the contemporary nature 
of these events was interesting, it was not suggested that they were linked; "one 
swallow does not make a summer" . 

Since then the Sun-Earth connection has been the source of detailed inves- 
tigation, with the Sun shown to have a major influence on the Earth (beyond 
just heat and light). The advent of larger and more powerful ground- and space- 
based telescopes has shown that the Sun, once thought of as a serene globe around 
which the world revolves (or which revolves around the Earth), is instead a tur- 
bulent roiling ball of plasma and relatively powerful magnetic fields. The outer 
solar atmosphere, once only visible during the pseudo-night of a solar eclipse, is 
now observed regularly using artificial means to block the face of the Sun, re- 
vealing magnetic structures and vast ejections of plasma that are fired into the 
heliosphere. 

Scientifically, the Sun is classified as a G2V yellow dwarf star found on the 
main sequence of the Hertzsprung-Russell diagram (a graphical representation of 
stars obtained by plotting the absolute magnitude against the spectral class; see 
Figure 1.1) with a mass of ~2xl0 30 kg and a radius of ~6.96xl0 8 m. The Sun 
has a chemical composition consisting primarily of hydrogen (~75 %) and helium 
(~25 %) with some trace amounts of heavier elements. The high concentration 
of hydrogen, combined with its large size, provides the conditions necessary for 
the production of energy via nuclear fusion deep within the heart of the Sun. 
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Figure 1.1: The Hertzsprung-Russell diagram. The Sun is a G2 star located 
on the main sequence, with a surface temperature of ~6000 K. Note that stellar 
luminosity values are measured relative to the solar luminosity (L = L s j ar /L ; 
picture courtesy heasarc.gsfc.nasa.gov). 



1.1 The Solar Interior 

The solar interior consists of three main regions; the core, the radiative zone 
and the convective zone (see Figure 1.2). Each region is determined using the 
variation in physical values such as density and temperature, with the result that 
their absolute limits are not well defined. 

In the core, the temperatures and densities are sufficiently high (~10-15 MK 
and ~ 10-150 g cm~ 3 respectively) that proton-proton chain nuclear fusion can 
occur. This typically occurs through the pp I chain 1 according to the steps, 

1 H+ 1 H ->■ 2 H + e + + v (1.1) 
2 H+ 1 H ->■ 3 #e + 7 (1.2) 

Tt is also possible to get fusion through the pp II chain using Li VII as an intermediate. 
However this process typically requires higher temperatures than those available and as a result 
makes a minimal contribution to the energy produced in the core. 
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Figure 1.2: Cartoon showing the different regions of the Sun along with 
some of the features observable in the solar atmosphere (picture courtesy so- 
howww.nascom.nasa.gov). 



3 He+ 3 He -> 4 #e + 2 x # + 7 (1.3) 

where l H is a hydrogen nucleus, 2 H is the deuteron isotope of hydrogen, 3 He and 
He are the 1- and 2-neutron isotopes of helium respectively, e + is a positron, v 
is a neutrino and 7 is a gamma-ray photon. The overall reaction therefore looks 
like, 

4 1 H ->■ 4 i/e + 2e + + 2z/ + 7. (1.4) 

The positrons (e + ) are then annihilated through interaction with electrons after a 
very short period of time, producing 1.02 MeV. One set of these interactions pro- 
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duces a total energy of ~26.7 MeV, emitted primarily as high-energy 7-ray pho- 
tons. Although a temperature of ~10 10 K is required to overcome the Coulomb 
repulsion between the protons through purely thermal interactions, this reaction 
can occur at lower temperatures comparable to those in the core by including 
quantum mechanical effects (cf. the Gamow factor; Gamow, 1928). 

This process occurs within ~0.25 R of Sun centre where the temperature 
and density are conducive to fusion. At ~0.25 R , the temperature and density 
drop sufficiently low that fusion can no longer occur. This decrease in density 
makes the plasma optically thin, allowing the transfer of energy outwards from 
the core through radiative diffusion. However, the temperature is still too high to 
allow neutral atoms to form, with the result that the energetic photons produced 
by fusion in the core are scattered by the free electrons, neutrons and protons in 
the radiative zone. This scattering results in the photons performing a random 
walk through the radiative zone that can take ~10 4 years to complete. 

At ~0.7 R (Christensen-Dalsgaard et al, 1991), the temperature drops suffi- 
ciently to allow the formation of atoms, with the plasma becoming opaque. The 
constant influx of energetic photons from the radiative zone produces a tempera- 
ture gradient, causing the plasma to become convectively unstable. This results 
in the formation of giant convective cells (akin to those found in a boiling kettle) 
that carry energy towards the solar surface. The density continues to drop with 
distance from Sun centre until the plasma is no longer opaque. The point at 
which this occurs for visible light defines the surface of the Sun, where the energy 
carried by convective cells from deep within the Sun can be radiated away as 
visible photons. 

As a massive body containing constantly moving charged ions, the Sun is 
magnetically active with a magnetic field that can be treated as a dipole similar 
to that of the Earth. However, unlike the Earth's magnetic field, it does vary 
from a dipolar to quadrupolar configuration on an 11 year timescale. This is 
a result of the variation of the solar rotation rate from that of a solid body in 
the core and radiative zone to a differential rotation in the convective zone. The 
interface between these two rotation types is called the tachocline, and it is at this 
point that the solar magnetic field is wound up by the differential solar rotation. 
The magnetic field lines then become buoyant, and are carried to the surface 
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by convective motion before emerging into the solar atmosphere (the a-Q effect; 
Babcock, 1961). 



1.2 The Solar Atmosphere 

The solar atmosphere, similar to the Earth's atmosphere, is defined as the region 
starting above the visible solar surface and extending out into the heliosphere 
(defined as the region in space directly influenced by the solar wind), technically 
to the heliopause at the edge of the solar system, where the solar wind pressure 
is equal to the pressure of the ambient interstellar medium. It consists of a 
tenuous, highly structured, ionised plasma that is threaded and modified by the 
solar magnetic field. Although difficult to see due to the dominant photospheric 
emission, the structure of the outer atmosphere can be discerned during a solar 
eclipse. 

The different regions of the solar atmosphere are shown in Figure 1.2, while 
the variation in the density and temperature with height is shown in Figure 1.3. 
Each of the regions of the solar atmosphere are defined according to variations in 
density and temperature and these are noted in Figure 1.3. However, the dynamic 
nature of the solar atmosphere means that these are not firm definitions with 
clearly defined layers but general trends defined by different plasma properties. 

1.2.1 The Photosphere 

The photosphere is the visible surface of the Sun, with a black-body temperature 
of ~6000 K. As it is situated at the top of the convective zone, the photosphere 
is home to small-scale granulation arising from small convective flows resulting 
from temperature gradients below the surface. These features are visible in Fig- 
ure 1.4, which shows a close-up image of a sunspot taken using the Swedish Solar 
Telescope (SST) in the Canary Islands. The dark sunspot is seen in the centre of 
the image, with smaller dark pores around it while the mottled "bubbles" show 
the granulation resulting from the small-scale convective flows. The solar mag- 
netic field also begins to play a significant role here, with the buoyant magnetic 
field lines from deep within the Sun emerging through the photosphere. Upon 
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Figure 1.3: Plane parallel model showing the variation in solar temperature and 
density with height. The different regions of the solar atmosphere are marked, 
with both temperature (T e ; K) and density (N e ; cm~ 3 ) variation shown using 
logarithmic scales. Figure from Gallagher (2000) after Gabriel &; Mason (1982). 



breaching the photosphere, these regions of concentrated magnetic field (which 
are cooler than the surrounding photospheric regions) appear as dark sunspots. 

The photosphere comprises a very thin portion of the solar atmosphere at a 
temperature of ~6000 K. This region of the Sun is best observed using visible 
light such as that observed by the Solar Optical Telescope (SOT) onboard the 
Hinode spacecraft or by ground-based telescopes such as the SST. 



1.2.2 The Chromosphere 

Above the photosphere, the temperature and density both continue to drop with 
height. Although the density continues to drop into the outer atmosphere (see 
Figure 1.3), the temperature reaches a minimum at a height of ~200 km above 
the photosphere before beginning to increase once more. The region of the solar 
atmosphere where this occurs is called the chromosphere, taking its name from the 
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Figure 1.4: A close-up of a sunspot in the photosphere observed by the Swedish 
Solar Telescope in the Canary Islands (figure courtesy The Royal Swedish Academy 
of Sciences/The Institute for Solar Physics). 



coloured appearance it has during a solar eclipse (from the Greek word khroma, 
meaning colour). 

Despite the gradual increase in temperature, the steady drop in density with 
height in the chromosphere results in a drop in the gas pressure of the chromo- 
sphere (as p gas = n e k B T e , where n e is the density, ks is the Boltzmann constant 
and T e is the temperature). This drop in the plasma pressure produces a situation 
where the magnetic pressure (in cgs units, ps = B 2 /8tt, where B is the magnetic 
field strength in gauss, or in SI units, ps = B 2 /2fi, where B is the magnetic field 
strength in tesla and fi is the permeability of free space) approaches and begins to 
dominate the gas pressure, with the result that the magnetic field begins to dom- 
inate the structure of the chromosphere. This is illustrated by the variation of 



8 



the plasma-/? value from the chromosphere into the corona, where the plasma-/? 
is defined as, 

The gradual domination of the solar magnetic field can be seen in observa- 
tions of the solar chromosphere, which exhibits a very dynamic structure in all 
observable passbands including Ha and EUV emission (such as 304 A; panel c 
of Figure 1.5). The relatively cool nature of the chromosphere means that it has 
been possible to observe it from ground-based telescopes using Ha filters for many 
years, with many observatories taking daily synoptic Ha observations of the solar 
chromosphere (a practice that still continues today, cf. Kanzelhohe Observatory 
in Austria). 

It was using Ha observations of a solar flare that the first globally-propagating 
disturbance in the solar atmosphere was observed by Moreton (1960). Called the 
"Moreton wave" , this disturbance was observed to have a velocity of ~1000 km s _1 
and be strongly associated with solar flare eruptions and Type II 1 radio bursts 
in the solar corona, indicating a moving shock. There were some issues with the 
observations, particularly regarding the high velocity of the pulse, which greatly 
exceeds the estimated Alfven speed for the solar chromosphere. To overcome this 
problem, Uchida (1968) proposed that the Moreton wave was the chromospheric 
signature of a fast-mode magnetoacoustic wave propagating in the solar corona; 
this is discussed in more detail in Section 2.1. 

The region between the chromosphere and the corona is known as the tran- 
sition region (see Figure 1.3), with the temperature showing a dramatic increase 
over a height range of ~100 km. The reason for this sudden temperature increase 
(the coronal heating problem) remains uncertain, and is one of the major topics 
of investigation in modern astrophysics. However, there are two main branches of 
theories that have been proposed to resolve this issue. The first deals with the dis- 
sipation of magnetic stresses caused by the gradual displacement of coronal loop 
foot-points by photospheric motions, producing direct heating (the direct current 

1 There are five types of radio bursts typically observed in solar physics. Type I are as- 
sociated with solar flares and active regions, Type II indicate moving shocks in the corona, 
Type III indicate particle emission, Type IV correspond to stationary shocks and Type V 
indicate short-lived continuum emission 
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(a; 21 1A) 



(b; 171A) 




Figure 1.5: Close-up images showing the dynamic nature of the solar atmosphere 
observed in the Extreme UltraViolet (EUV) by SDO / AIA. Panel (a) shows the 
corona (211 A; T ~2 MK), panel (b) shows the upper transition region and low 
corona (171 A; T ~0.7 MK), while panel (c) shows the upper chromosphere and 
lower transition region (304 A; T ~0.06 MK). Each panel has been coloured using 
the relevant colour tables for clarity. 
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Figure 1.6: The structure of the solar atmosphere based on an energy balance 
model proposed by Gabriel (1976). The magnetic field lines emerging through the 
photosphere fan out in the chromosphere and expand radially outwards from the 
Sun through the corona. 

or DC models). The other branch of theories proposes that the displacement of 
coronal loop foot-points occurs on a rapid timescale compared to the end-to-end 
Alfven travel time, generating waves which dissipate to heat the local coronal 
plasma (the alternating current or AC models). For more information, see the 
review of coronal heating processes by Klimchuk (2006). 

The anomalous nature of the upper chromosphere and transition region has 
important implications for a thorough analysis of the low corona, with the diffuse 
structures found in the corona having their foundation in this region. The con- 
centrated magnetic field lines emerging from the photosphere fan out through the 
chromosphere into the corona as a consequence of the decreasing plasma-/? value. 
This is shown in Figure 1.6, which provides a graphical outline of the structure 
of the solar atmosphere based upon an energy balance model (Gabriel, 1976). 
Although this has since been superseded by a more dynamic model, it remains a 
useful guide to the basic structure of the solar atmosphere and the relationships 
between the various atmospheric regions. 
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1.2.3 The Corona 



The solar corona comprises the outer part of the solar atmosphere, extending from 
~2000 km above the solar photosphere out into the heliosphere. The intensity of 
the corona is several orders of magnitude less than that of the Sun itself, with the 
result that it is best observed by blocking out emission from the solar disk. The 
best example of this occurs during a total solar eclipse, when the near-perfect 
alignment of the Moon and the Sun leads to the solar disk being blocked, thus 
allowing the corona to be observed (see Figure 1.7). 

As previously shown in Figure 1.3, the temperature of the corona is ~10 6 K, 
much hotter than the photosphere and chromosphere below it, while the density 
is ~10 6 -10 9 cm -3 . As a result, much of the emission from the low corona comes 
from highly ionised heavy elements such as iron. This emission is primarily in the 
extreme ultraviolet range of the electromagnetic spectrum, and consequently must 
be observed using space-based instruments. This is apparent in Figure 1.8, which 
shows the EUV corona as observed in the 171 A passband which is dominated by 
emission from the Fe IX line. 

The solar corona is a low-/? plasma, and as a result is dominated by the 
solar magnetic field, with the plasma "stuck" to the magnetic field lines. This 
is apparent in Figure 1.7, where the magnetic structure of the inner and outer 
corona is quite visible, particularly above the bright active regions and at the 
north and south poles. The fine structure in Figure 1.8 is also a result of the 
\ow-/3 nature of the corona. The dominance of the magnetic field is such that it 
defines the different regions within the corona; active regions, coronal holes and 
the quiet corona. 

Coronal active regions (seen as the brighter regions in Figures 1.7 and 1.8) are 
regions of intense magnetic field that are associated with photospheric sunspot 
regions. Here, magnetic field strengths can reach up to several hundreds of Gauss, 
while the interaction of different field lines can result in magnetic reconnection, 
producing solar flares and coronal mass ejections (CMEs). Coronal holes lie at 
the opposite end of the spectrum, with open magnetic field lines that radiate out 
into the heliosphere (this is most apparent at the poles of the Sun in Figure 1.7). 
This allows plasma to propagate away from the Sun along the open field lines as 
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Figure 1.7: A combined EUV and white light image of a total solar eclipse from 
2010 July 11 reproduced from Pasachoff et al. (2011). The EUV image of the 
Sun has been processed using a radial gradient filter, while the white-light corona 
has been processed using the technique proposed by Druckmiiller et al. (2006) to 
highlight features. 

the fast solar wind (cf. Parker, 1958), with this continuous loss of plasma a reason 
why coronal holes are darker than the surrounding corona in EUV passbands. 

The quiet corona is the intermediate region between active regions and coronal 
holes. Here, the topology consists of a combination of open and closed magnetic 
field lines and small scale structures that are changing constantly on a variety 
of timescales. This "magnetic carpet" can produce small solar flares and CME 
eruptions, with a density and temperature that is higher than a coronal hole, but 
lower than an active region. 
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Figure 1.8: An EUV image from the 171 A passband showing the solar corona as 
observed by SDO /AIA. The image has been displayed using the appropriate colour 
table for clarity. 



The concentrated magnetic fields and increased densities in active regions com- 
pared to the quiet Sun produce the increased emission apparent in Figures 1.5, 
1.7 and 1.8. Similarly, the low densities and open magnetic field lines result in 
coronal holes being much darker than the surrounding quiet Sun regions when ob- 
served using EUV and X-ray images. The different passbands typically observed 
by space-based observatories allow these different regions of the solar atmosphere 
to be studied by examining different emission lines. This is most apparent in 
Figure 1.5, where the decrease in temperature from panel (a) to (c) highlights 
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different structures and regions of the atmosphere. 

The heights of the different regions of the solar atmosphere observed by each 
passband can be estimated given the temperature of the plasma observed. This 
is done following the example of Chapman & Zirin (1957) and assuming that the 
solar corona is in hydrostatic equilibrium. In this case, 



dpgas 
~dh~ 



-pg (1.6) 



which indicates that the change in gas pressure p gas with height h is dependent 
on the mass density p and solar gravity g. The solar gravity can be calculated 
using Newton's Law as, 

GM Q 

where G is the gravitational constant (G = 6.67 x 10~ n m 3 kg -1 s~ 2 ), M is the 
solar mass (M = 2 x 10 30 kg) and R is the solar radius (R©= 6.96 x 10 8 m). 
This produces an acceleration due to gravity of ~270 m s~ 2 at the solar surface. 
The mass density p can be determined using the ideal gas law p gas = nksT, 
where n is the particle number density, ks is the Boltzmann constant (ks = 
1.38 x 10 -23 m 2 kg s -2 K _1 ) and T is the temperature. The ideal gas law can 
then be rewritten as, 

Pgas = nk B T = > p = -j-^Pgas, (1-8) 

where m is the proton mass, 1.67 x 1CT 27 kg. This can then be substituted into 
Equation 1.6 to give, 

rlr> m n 

Pgas (1-9) 



dp gas _ rng 



dh k B T 

If isothermal equilibrium is assumed, this can be solved by rearranging the equa- 
tion and integrating as, 



P9as dp gas mg r h 



PC 



Pgas ksT Jq 



dh, (1.10) 



where p is the surface gas pressure and Pgas is the gas pressure at the height h 
above the surface. Solving this equation gives the estimate of gas pressure p gas 
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at a height h above the surface as, 



P gas (h) =p e 



h/H 



(1.11) 



where H is the scale height and is defined as, 



H = 



(1.12) 



Using the values defined above for a typical coronal temperature of ~l-2 MK 
produces an estimated scale height of ~30-60 Mm, while typical chromospheric 
temperatures of ~10 4 -10 5 K produce scale heights of 0.3-3 Mm. These values 
correspond quite well with those estimated in Figure 1.3, indicating that the 
passbands shown in Figure 1.5 can be used to study the upper chromosphere and 
low corona. 

Although the temperature and density of the corona can be estimated using 
emission lines, the magnetic field in the corona is more difficult to study. The 
effects of the magnetic field are quite apparent in the solar structures shown in 
Figure 1.8, but the strength of the magnetic field must be inferred using photo- 
spheric magnetograms or by calculating the speed of coronal shocks using Type II 
radio bursts. These techniques produce an estimated magnetic field strength of 
~10 G in the quiet corona, with values of the order ~100 G in active regions. 

While the magnetic field strength, temperature and density variation in the 
corona can be estimated using various techniques, these estimates do have large 
associated errors. A technique that could potentially be used to probe the physical 
structure of the corona directly is coronal seismology, which uses the oscillation of 
coronal structures such as loops to estimate the magnetic field strength and den- 
sity (e.g., Ballai & Erdelyi, 2004). However, this technique requires high temporal 
cadence observations of the solar corona with very high spatial resolution. 

An understanding of the magnetic field strength variation is vital for under- 
standing the structure of the solar corona as a direct consequence of the fact that 
the corona is a low-/3 plasma. Although the structure of the corona is undeniably 
complex (as shown in Figures 1.5, 1.7 and 1.8), it is possible to derive a mathe- 
matical framework that can be used to study it. This is done by combining the 
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electromagnetic equations of Maxwell (1865) with the fluid equations, producing 
the magnetohydrodynamic equations. 

1.3 Magnetohydrodynamics 

The concept of examining the electromagnetic-hydrodynamic waves resulting 
from the motion of a conducting liquid in a constant magnetic field was first pro- 
posed by Alfven (1942). This involves the use of magnetohydrodynamics (MHD), 
which combines the theories of electromagnetism and fluid dynamics, allowing the 
evolution of magnetic plasma to be studied. At its most basic level, MHD theory 
quantifies the interplay of forces involved in the motion of a conducting plasma 
in the presence of a magnetic field. As a result, MHD is of vital importance for a 
thorough analysis of the solar corona, which is a highly ionised plasma threaded 
by the solar magnetic field and in a constant state of flux. In this section, the 
MHD equations will be defined and used to examine the types of waves that can 
result in this environment. Note that a number of different books were used as 
reference sources for this section; these may be found in the Bibliography. 

Maxwell's equations constitute the bulk of the electromagnetic contribution to 
the set of MHD equations and allow the relationship between the electric currents 
and magnetic fields to be quantified. They are typically given as, 

VxE = (1.13) 

V-E = -p (1.14) 
VxB = /i j (1.15) 
V-B = (1.16) 

where E is the electric field vector, B is the magnetic field vector, e is the 
permittivity of free space, p is the charge density, po is the magnetic permeability 
and j is the current density. Faraday's law (Eqn. 1.13) indicates that a spatially 
varying electric field can induce a magnetic field, while Ampere's law (Eqn. 1.15; 
simplified here by assuming that the plasma velocity in an MHD system is much 
smaller than the velocity of light) shows that gradients in a magnetic field induce 
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electric currents. Gauss' law (Eqn. 1.14) covers the conservation of electric charge, 
while the solenoidal constraint (Eqn. 1.16) indicates that magnetic monopoles 
cannot exist. The generalised Ohm's law is also included in addition to the above 
equations, 

j = a(E + vxB) (1.17) 

where a is the electrical conductivity and v is the velocity of the plasma. This 
allows the fluid and electromagnetic equations to be coupled through the shared 
plasma velocity parameter. 

When considering MHD in solar physics, it is often useful to replace the 
electric field E with the magnetic field B. This can be done using Ohm's Law 
(Equation 1.17) and Ampere's law (Equation 1.15) as, 

1 1 

E = -j-vxB = VxB-vxB. (1.18) 

a fx a 

This can then be substituted into Faraday's Law (Equation 1.13) to give, 

r)B 

— = V x (v x B) - V x (r]V x B), (1.19) 

where the magnetic diffusivity i] = Xj^a. This is assumed to be constant, and 
vector identities are then used to reduce this equation to, 

— = V x (v x B) + r]V 2 B. (1.20) 
at 

This is the induction equation, and is an important equation in MHD as it allows 
the magnetic field B to be calculated if the velocity v is known. In addition, it can 
be used to determine the behaviour of the magnetic field through the Magnetic 
Reynolds number R m . This is defined as, 

_ V x (v x B) 
Rm ~ V V*B ' (L21) 

and indicates the ratio of the advection and diffusion terms of the induction 
equation. If R m is larger than unity, then the magnetic field is frozen into the 
plasma and the motion of the magnetic field is defined by plasma motion, whereas 
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if R m is smaller than unity, then the magnetic field tends to diffuse through the 
plasma. R m is analogous to the Reynolds number used in fluid dynamics, and 
can be estimated using, 

Rm ~ ^, (1-22) 

rj 

where 77 = 10 9 T~ 3//2 m 2 s -1 . Using typical coronal values of T~10 6 K, l ~10 5 m 
and vo ~10 3 m s" 1 gives a value of R m ~10 8 . The magnetic field in the solar 
corona is therefore dominated by plasma motion. 

As well as the electromagnetic equations, the fluid equations are also used in 
MHD to study the mass motion of the plasma. The conservation of mass of the 
plasma is defined as, 

f + V-(pv) = (1.23) 

where p is the plasma density. This indicates that the density at a point increases 
if mass flows into and decreases if mass flows out of the surrounding region. 
Equation 1.23 can also be rewritten using vector identities to give, 

|^+v- Vp + pV-v = 0. (1.24) 

However, for an incompressible fluid along a path with velocity v the convective 
time derivative D/Dt = d/dt + v ■ V = 0. As a result, Eqn. 1.24 reduces to, 

V-v = 0, (1.25) 

indicating that there are no sources or sinks in an incompressible fluid. 
Next, the equation of motion of the fluid must be considered, 

P^ = -Vp + F (1.26) 

where p is the pressure and F is the external force acting on the fluid. In this 
case, the external force consists of the Lorentz force j x B and the gravitational 
force pg. Equations 1.26 can therefore be written as, 

<9v 

9 dt + piy ' V)v = ~ Vp + j X B + pg (1 ' 27) 
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where the electromagnetic and fluid equations have been combined through the 
Lorentz force and D/Dt has been rewritten as shown above. 

As well as the equations outlined above, the MHD equations also include the 
ideal gas law, 

p = nk B T (1.28) 

where n is the density, ks is the Boltzmann constant (1.38 x 1CT 23 m 2 kg s -2 K _1 ) 
and T is the temperature. This allows the plasma to be treated as an ideal gas, 
simplifying analysis. The energy equation for the plasma is also included, 

where 7 is the ratio of specific heats (7 = c v jc v x ] usually taken as 5/3) and the 
plasma has been assumed to be adiabatic in order to simplify calculations. This 
can then be rewritten as, 

| + ( v .V)p=- 7 pV-v (1.30) 

using vector identities and Equation 1.23. 

With the MHD equations defined, they can now be used as a basis from which 
to examine the evolution of a magnetised plasma and identify the different waves 
that result from this evolution. 



1.3.1 Acoustic Waves 

To determine the different MHD wave types possible in a magnetised plasma, 
the above equations are linearised by considering small perturbations. This is 
achieved by assuming that the wave amplitude is small and that each parameter 
consists of an equilibrium value (with subscript zero in the following equations; 
Bo) and a perturbed value with a spatial and temporal component (with subscript 
one; Bi). The different parameters can therefore be described as, 

B = Bo + Bxfot) (1.31) 
1 c p is the specific heat at constant pressure while c v is the specific heat at constant volume 
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v = v + vi(r,t) = vi(r,t) (as v = 0) (1.32) 
P = Po + pi(r,t) (1.33) 
p = p + Pl (r,t) (1.34) 

With the aid of these assumptions, the MHD equations presented above can be 
linearised to give, 



9Bi 



= V x (vi x B ) (1.35) 



dt 

V-Bi = (1.36) 



dp 



d I + V • ( Po vi) = (1.37) 

Po^ = -Vpi + -(V x Bi) x B 
at n 

+ -(VxB ) xBx + pxg (1.38) 

^ + vi • Vp = -7PoV • vi (1.39) 



The first step in determining the different wave types present in a magnetised 
plasma is an identification of the properties of acoustic waves. To do this, it 
is assumed that the plasma is not magnetised (B = 0) and that there are no 
gravitational effects (g — 0). The equilibrium pressure gradient is also assumed to 
be zero, giving a constant equilibrium pressure (po = constant and p = constant). 
The linearised MHD equations then become, 

(1.40) 
-Vpi (1.41) 
-7PoV • Vl (1.42) 

Fourier analysis is then used to search for plane wave solutions to these equations 
by assuming a solution of the form Ce 1 ^^^, where k is the wavenumber, u> is 
the angular frequency and C is a constant. This form of analysis is used as the 



°Pi „ 

dpi 
dt 
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mathematical operators may be simplified as, 



d_ 
di 



d 



= ik 



dt 2 

d 2 



= —LU 



= -k 2 



dx dx 2 
ik- Vx = ik x V = ik 



1.43) 



Equations 1.40 to 1.42 can then be solved using this approach, giving the 
equations, 

vi = — k (1.44) 

k- Vl = ^ (1.45) 
Po 

k- Vl = ^. (1.46) 

7Po 

These equations indicate that the result is a longitudinal wave, with the velocity 
vector v parallel to the wave vector k. The resulting dispersion relation for this 
plane wave solution can then be found by rearranging equations 1.44 to 1.46 to 
give, 

oo 2 = k 2 c\ (1.47) 
where c s is the sound speed of the plasma and has a value given by, 

4 = TSa. (i.48) 

Po 

It is possible to estimate the sound speed in the solar corona using typical coro- 
nal values for this equation. From the ideal gas law in Equation 1.28 above, 
Equation 1.48 can be rewritten as, 




c s = J-^-. (1.49) 



For T~10 6 K, 7 ~1, k B = 1.38 x 1(T 23 m 2 kg s~ 2 Kr 1 and m p ~1(T 27 kg, this 
produces an estimated sound speed of c^ ~100 km s -1 for the solar corona. 
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This indicates that there is a wave solution to the MHD equations that de- 
pends entirely on the pressure and density of the plasma, i.e., an acoustic wave. 
Analogously, it is possible to ignore the effects of pressure and density and deter- 
mine if there is a purely magnetic wave solution to the MHD equations. 



1.3.2 Magnetic Waves 

To identify pure magnetic waves, density and pressure variations must first be 
ignored, setting p = 0, po 7^ 0, p\ = const, p\ = const and g = 0. For simplicity 
of calculations, the magnetic field is also assumed to be uniform and directed in 
the z direction. The linearised MHD equations (equations 1.35 to 1.39 above) 
then become, 



V-Vi = (1.50) 
dB 1 



Po 9 ^ = ^(VxBOxBo (1.51) 



V x (vi x B ) (1.52) 

V-Bi = (1.53) 

Once again, plane wave solutions to these equations may be found using 
Fourier analysis, giving, 

ik-vi = (1.54) 

-iup vi = -(ik x Bi) x B (1.55) 

-iwBi = ik x (vi x B ) (1.56) 

ik-B 1 = 0. (1.57) 

Here, the plane wave solution has a velocity vector Vi that is perpendicular to the 

wave vector k, and also to the direction of the magnetic field vector B^ These 
equations may then be manipulated using vector identities to produce, 

«• = (1.58) 
Wo 
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which can then be used to produce the dispersion relation for a purely magnetic 
wave, 

u = ±v A k z (1.59) 

where k z is the wave-vector in the z-direction. The wave velocity in this case (va) 
is called the Alfven speed and is derived as 

va = -^=, (1.60) 

in S.I. units, where fj, is the magnetic permeability. It is typically assumed that 
the plasma is fully ionised, with solar physics also using cgs nomenclature by 
convention, allowing Equation 1.60 to be rewritten as, 

va= y= , (1.61) 

where -Bo is in units of gauss, rii is the ion number density (cm~ 3 ), and m p is 
the proton mass (g). Once again, it is possible to estimate a value for the Alfven 
speed in the solar corona using typical coronal values. For rii ~10 8 cm -3 , m p = 
1.67 x 10~ 24 g and B ~10 G, this produces an Alfven speed of v A ~2000 km s" 1 
for the corona. 

These waves are called Alfven waves and are pure magnetic waves that result 
from the oscillation of magnetic field lines. They are anisotropic, incompressible 
waves that propagate along magnetic field lines. 



1.3.3 Magneto-acoustic Waves 

As the solar corona contains variations in magnetic field, density, and pressure, 
any waves present will contain a contribution from both acoustic and magnetic 
wave modes, giving magneto-acoustic waves. Identifying these waves requires the 
combination of magnetic and gas pressure, although gravitational effects can be 
ignored as they are negligible in comparison. The linearised MHD equations (par- 
ticularly Equations 1.35, 1.38 and 1.39) are then used to determine the dispersion 
relation for a magnetoacoustic wave. 

Both the magnetic and gas pressure are assumed to be uniform and of the form 
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B = B z for the magnetic pressure and po = C (where C is a constant) for the 
gas pressure. By combining Equations 1.35, 1.38 and 1.39 with the mathematical 
operators shown in 1.43, it is possible to derive the simplified equation, 

oo 4 - oo 2 k 2 (c 2 s + v\) + c 2 v 2 A k 4 cos 2 6 B = 0, (1.62) 

where 6b is the angle between the wave- vector k and the direction of the magnetic 
field B. This is a quadratic equation in terms of ou 2 /k 2 and may be solved using 
the quadratic formula to produce the dispersion relation for magnetoacoustic 
waves, 

i 

UJ ' ' 



k 



Q(c 2 + v 2 A ) ± ^ V / cf + <-2cMcos2^ 2 . (1.63) 



There are two distinct solutions to this equation. Solving for the negative sign 
produces the slow-mode MHD wave, while solving for the positive sign produces 
the fast-mode MHD wave. The equation also reduces to the dispersion relations 
for an acoustic and magnetic wave when magnetic and pressure effects respectively 
are ignored. 

To examine the differences between the fast-mode and slow mode MHD waves, 
first consider the case where 6b = (i.e., k || B). This gives the solutions, 

^ = V fast-mode = V A and — = V s l ow - mo de = C s (1.64) 

for the fast-mode (left) and slow-mode (right) MHD waves. This means that 
the fast-mode wave travels at the Alfven speed while the slow-mode wave travels 
at the local sound speed when propagating along the magnetic field. For the case 
where 6 B = tt/2 (i.e.: k _L B) the solutions are, 



^fast-mode = ^s + V A &nd ^ = V slow _ mode = (1.65) 



for the fast-mode (left) and slow-mode (right) respectively. This means that the 
fast-mode wave velocity when propagating across the magnetic field is dependent 
on both the sound speed and the Alfven speed, while the slow-mode wave cannot 
propagate across the magnetic field. 

The solutions to equation 1.62 for all angles 6 are shown in Figure 1.9 for the 
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Figure 1.9: Polar plot showing the phase velocities of magnetoacoustic waves 
for a ratio c s /va = 0.7. The sound and Alfven speeds are shown as solid lines, 
with the slow-mode, Alfven and fast-mode wave speeds shown as dashed lines (cf. 
Aschwanden, 2004). 

ratio c s /va = 0.7 with the Alfven and sound speeds also shown for comparison. It 
is clear that the Alfven and sound velocity (va and c s ) do not show any variation 
as the orientation of the magnetic field changes. In contrast, the fast-mode, 
slow-mode and Alfven wave velocities are highly dependent on the magnetic field 
orientation, a result that has major implications for the study of the solar corona. 

It should be noted that in ideal MHD as discussed here the different wave ve- 
locities are constant, with any variation between events dependent on the plasma 
parameters. These waves are also non-dispersive and non-dissipative, with the 
result that the MHD equations must be extended to include resistive terms when 
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considering phenomena such as flares or magnetic reconnection in the solar atmo- 
sphere. However, recent work by Murawski et al. (2001) and Nakariakov (2001) 
suggests that the randomly structured nature of the solar corona may result in 
the deceleration, dispersion and dissipation of a propagating MHD wave. This 
behaviour was observed using simulations of MHD waves and may be applicable 
to observations of the solar corona. The variation of the morphology of larger- 
amplitude waves in a plasma has been studied by Vrsnak & Lulic (2000a,b). 
In these examples, the pulse exhibits non-linear behaviour, with a propagation 
velocity that is dependent on the amplitude, while deceleration, dispersion and 
dissipation of the pulse follow as a natural consequence. 

1.4 Thesis Outline 

The work outlined in this thesis deals with the physical properties of large-scale 
propagating disturbances observed in the low solar corona, with particular at- 
tention paid to their kinematics and temporal variation in morphology. The 
kinematics of these disturbances remain a source of much debate as a conse- 
quence of conflicting observations and the relatively low observing cadence of the 
instruments used to study them. It is shown here that this low observing cadence 
may have produced an underestimation of the kinematics previously estimated 
for these pulses. The different techniques traditionally used to derive the pulse 
kinematics are also examined and compared to a bootstrapping technique, with 
the bootstrapping technique shown to produce a better estimate of the true pulse 
kinematics. A semi-automated technique for identifying, tracking and analysing 
this phenomenon is presented and used to show that the disturbances exhibit 
clear deceleration and pulse dispersion with propagation. Finally, the dispersion 
relation of a pulse is calculated using high cadence observations of a CBF pulse. 

An overview of previous observations of globally-propagating disturbances 
in the corona is presented in Chapter 2, with the different theories proposed to 
explain the phenomenon also discussed. The different instruments used in this 
work are outlined in Chapter 3, along with a discussion of the software used to 
process the data. The first observations of a globally-propagating disturbance 
made by the STEREO spacecraft are presented in Chapter 4, with some of the 
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issues resulting from the analysis of this event discussed in Chapter 5. A semi- 
automated technique for studying these disturbances outlined in Chapter 5 is used 
in Chapter 6 to provide an unbiased analysis of four separate events observed by 
the STEREO spacecraft. This technique is also applied to an event observed by 
the SDO spacecraft in Chapter 7 and used to examine the physical properties of 
the pulse and the corona itself. Finally, the major conclusions of this thesis and 
some possible directions for future work are provided in Chapter 8. 
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Chapter 2 

Observations and Theory of 
Coronal Bright Fronts 

/ am not accustomed to saying anything with certainty after only one or two 
observations. 

Andreas Vesalius 

In this chapter, a brief history of the observations of globally-propagating 
disturbances in the solar atmosphere is presented along with an overview of the 
myriad theories proposed to explain them. The first observation of a globally 
propagating disturbance in the EUV solar corona was made by Moses et al. 
(1997) using the Extreme ultraviolet Imaging Telescope (EIT; Delaboudiniere 
et al, 1995) onboard the SOlar and Heliospheric Observatory (SOHO; Domingo 
et al, 1995) spacecraft. They noted the appearance of "a localized increase in 
emission propagating across the solar disk" associated with a Coronal Mass Ejec- 
tion (CME) eruption from 1996 December 23. The observed pulse was labelled 
an "EIT wave" , and became a source of much interest as a result of its association 
with other solar phenomena such as flares and CMEs. 

"EIT waves" have been studied in detail for almost 15 years and while ad- 
vances have been made with regard to their physical interpretation, many unan- 



29 



Figure 2.1: Running difference images showing the 1997 May 12 "EIT wave" 
event studied by Thompson et al. (1998). The time of the leading image is given 
in the bottom left of each panel. 

swered questions still remain. The various observations of these pulses are pre- 
sented in Section 2.1, with the two main branches of theories proposed to explain 
them outlined in Section 2.2. This section is complemented by recent reviews 
of "EIT waves" by Warmuth (2007, 2010), Wills-Davey & Attrill (2010), Vrsnak 
& Cliver (2008) and Zhukov (2011), with much of the material used here taken 
from the review paper published in Space Science Reviews by Gallagher & Long 
(2011). 

It should be noted at this point that "EIT waves" are so-called for their wave- 
like appearance and their first observation using S OHO /EIT. Since their initial 
discovery, they have been observed using myriad instruments across different 
passbands, while their true physical nature remains ambiguous. Although this 
name is still widely used for historical reasons, in this work we shall refer to "EIT 
waves" as "coronal bright fronts", or CBFs, following the lead of Gallagher & 
Long (2011). 

2.1 Observational Properties 

After the initial identification by Moses et al. (1997), an event was observed on 
1997 May 12 (see Figure 2.1), and studied in much greater detail by Thompson 
et al. (1998). This is traditionally treated as the original CBF identification 
paper, as it classifies the different coronal signatures associated with the event, 
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including the formation of dimming regions close to the event origin, a post- 
eruption arcade and the bright wavefront propagating quasi-radially away from 
the source region. 

The work of Thompson et al. (1998) was subsequently followed by a detailed 
analysis of a CBF observed using the Transition Region And Coronal Explorer 
(TRACE; Handy et al, 1999) spacecraft. Here, Wills-Davey & Thompson (1999) 
studied the higher cadence observations from both the 171 and 195 A passbands 
available using TRACE and found that the pulse propagated non-uniformly, with 
variations in both trajectory and velocity along the front of the pulse. The pulse 
was also noted to have a Gaussian form, indicating a single compressional front. 

With the assistance of the continuous data stream available from SOHO/ETT, 
the number of observed CBF events steadily increased, with numerous events 
observed and catalogued (see Thompson & Myers, 2009, for a catalogue of events 
from 1997-1998, but not published until 2009). These observations were used for 
statistical studies to try and determine the relationship of CBFs to other solar 
phenomena with either CMEs or flares thought to be the initial driver. 

This uncertainty in origin was overcome with the work of Biesecker et al. 
(2002) and Oliver et al. (2005) who used statistical studies of CBFs to study the 
relationship between CMEs, flares and CBFs. After correcting for statistical bias 
and observational effects, Biesecker et al. (2002) found that CBFs were always 
associated with a CME, whereas the converse was not necessarily true. This was 
followed by Cliver et al. (2005), who found that approximately half of the CBF 
sample studied were associated with small-scale flares of GOES intensities below 
C-class. This suggested that flares did not have the requisite energy to drive a 
CBF pulse, indicating that CMEs were the most likely driver. 

Although most apparent in EUV observations as a result of the data available 
from SOHO/EYT, globally-propagating disturbances have been observed in other 
passbands including soft X-ray (SXR; Hudson et al, 2003; Warmuth et al, 2005; 
Khan & Aurass, 2002), radio (White & Thompson, 2005), He I 10830 A (Gilbert 
& Holzer, 2004), and Ha (Moreton, 1960, as mentioned in Section 1.2.2). 

The relationship between CBFs and Moreton waves has been a particular focus 
of investigation, with many authors trying to identify the physical connection (if 
any) between the two phenomena. The first work to investigate Ha and EUV 
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Figure 2.2: Top: Distance-time plot from Warmuth et al. (2001) showing strong 
correspondence between Moreton wave and CBF kinematics for an event from 
1997 November 3. Bottom: Distance-time plot from Eto et al. (2002) showing the 
different kinematics of a Moreton wave and a CBF for an event from 1997 Novem- 
ber 4. 



observations of the same eruptive events (Warmuth et al., 2001) found compelling 
evidence that the Moreton waves and CBFs travelled co-spatially, suggesting 
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that they could be different forms of the same disturbance (see the top panel of 
Figure 2.2). However, this work was quickly repudiated by Eto et al. (2002), who 
found no overlap between the visibility times of a Moreton wave and a CBF for an 
event from 1997 November 4 (bottom panel of Figure 2.2), while a morphological 
and kinematical investigation of the same event showed no similarities between 
the different disturbances. Using a statistical approach, Okamoto et al. (2004) 
found that the majority of CBFs were not associated with a Moreton wave, while 
Warmuth et al. (2004a,b), found a 100 % correlation with CBF events using a list 
of 12 observed Moreton waves as a starting point. This suggests that while CBF 
events may be relatively common, an additional constraint must be met before a 
Moreton wave can occur. 

A possible link between the coronal EUV observations of a CBF and the 
chromospheric Moreton wave was observed by Vrsnak et al. (2002) using the 
chromospheric He I 10830 A emission line (see the top panel of Figure 2.3). 
This emission line provides a way to probe the upper chromosphere, potentially 
allowing it to be used to combine Ha and EUV observations, as initially proposed 
by Vrsnak et al. (2002). The importance of this emission line was emphasised 
by Gilbert et al. (2004), who observed 2 CBF pulses traveling co-spatially with 
He I pulses. These observations implied that the He I pulse was an imprint of the 
CBF pulse in the chromosphere as it propagates in the corona, rather than being 
a distinct wave formed by the same process as the CBF. Although it is tempting 
to view He I observations as the "missing link" between the Moreton wave and the 
CBF, the complex formation mechanism of the He I line, combined with the low 
observation rate in this passband makes this difficult (Gilbert & Holzer, 2004). 
However, these observations must be explained as part of any theory attempting 
to explain the CBF phenomenon. 

While globally propagating disturbances have been observed directly using 
radio data (cf. White & Thompson, 2005), the most apparent radio signature 
of a coronal disturbance is the Type II burst. A strong correlation was noted 
between Type II bursts and chromospheric Moreton waves (Kai, 1970), while a 
statistical survey of observed Type II bursts and Moreton waves found that ~70% 
of Moreton waves were associated with a type II burst (Pinter & Dryer, 1977). 
Klassen et al. (2000) compared a list of recorded Type II radio bursts with data 
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Figure 2.3: Top: Difference images from Vrsnak et al. (2002) showing a pulse 
observed using the He I emission line from 1998 August 24. Bottom: Plain images 
from Khan & Aurass (2002) showing a pulse observed on 1997 November 3 using 
soft X-ray emission from Yohkoh. 



from SOHO /EIT, finding that ~90 % of the observed radio events had an asso- 
ciated CBF, although the converse is not necessarily true (Biesecker et al, 2002). 
Microwave emission at 17 and 34 GHz from the Nobeyama Radioheliograph was 
used by White & Thompson (2005) to derive a velocity of ~830 km s _1 for a 
CBF event, leading the authors to suggest that the low cadence of SOHO /EIT 
may be producing an underestimation of the pulse kinematics. However, the 
observations in this case were poorly constrained with a low signal-to-noise ra- 
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tio. This poor signal-to-noise ratio is a major contributing factor to the sparse 
observations of radio pulses, although it may be rectified using modern radio ob- 
servatories such as the LOw Frequency ARray (LOFAR) or the Frequency Agile 
Solar Radiotelescope (FASR; Bastian, 2004) systems. 

Globally propagating disturbances have also been observed in soft X-ray 
(SXR) emission although they have been difficult to observe due to the rela- 
tively low signal-to-noise of the pulse. The first observation was made using data 
from the Yohkoh spacecraft (Khan & Aurass, 2002, bottom panel of Figure 2.3), 
with the studied event also showing a Moreton wave and a CBF. The X-ray wave 
was found to have a derived velocity of ~550 km s _1 , while the observation times 
of all three pulses were found to be comparable. The extrapolated leading edge 
of the X-ray wave was found to match the Moreton wave and CBF very closely, 
suggesting that all three disturbances were morphologically related. X-ray waves 
have since been observed by both Hudson et al. (2003) and Warmuth et al. (2005), 
with both authors suggesting that the emission is consistent with a propagating 
coronal shock wave. 

This interpretation has been questioned by Attrill et al. (2009), who com- 
bined X-ray observations from the X-Ray Telescope (XRT; Golub et al, 2007) 
onboard the Hinode spacecraft (Kosugi et al, 2007) with EUV observations from 
the Extreme UltraViolet Imager (EUVI; Wuelser et al, 2004) onboard the So- 
lar TErrestrial RElations Observatory (STEREO; Kaiser et al, 2008) spacecraft. 
The X-ray emission corresponding to an X-ray pulse was shown to match emis- 
sion from an EUV pulse, with both matching the edge of the erupting CME. 
These observations lead Attrill et al. (2009) to suggest that the X-ray and EUV 
emission was a result of the CME edge shocking the coronal plasma as it erupts 
rather than being due to a freely propagating wave. The rare nature of X-ray 
waves means that this remains an open question. 

2.2 Theories of Coronal Bright Fronts 

The multitude of observations across many diverse passbands outlined in the pre- 
vious section underlines the difficulties associated with elucidating "EIT waves" , 
and has resulted in the evolution of two sets of theories that attempt to explain 
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them. On one hand, CBFs are interpreted as waves using either magnetohydro- 
dynamic wave theory or shock wave theory. Alternatively, they are visualised 
as brightenings resulting from the reorganisation of the coronal magnetic field 
during the eruption of the associated CME. In this case, the "EIT wave" is not 
a true wave but instead a brightening due to small-scale magnetic reconnection 
or Joule heating. Both of these sets of theories are discussed below. 

2.2.1 Wave Interpretations 

A globally-propagating wavefront in the corona was initially proposed by Uchida 
(1968) to explain the Moreton wave phenomenon first observed by Moreton & 
Ramsey (1960). This was a propagating front moving in the chromosphere at 
velocities much greater than the predicted local Alfven velocity. Uchida (1968) 
proposed that the observed signatures could be explained by the propagation 
of a fast-mode MHD wave in the solar corona, expounding on this theory in 
later work (Uchida, 1970; Uchida et al., 1973; Uchida, 1974). As a fast-mode 
MHD wave propagates in the corona, the variation of the Alfven speed at the 
corona/chromosphere intersection should result in a circular intersection line that 
is observed as a Moreton wave. This matches the "up-down" swing in the Ha 
line observed by Dodson & Hedeman (1968) corresponding to the passage of a 
Moreton wave; indicative of a pulse compressing the chromosphere as it passes. 

With the discovery of the "EIT wave" in S OHO /EIT data (Moses et al., 1997; 
Dere et al, 1997) it was thought that this could be the coronal counterpart to 
the chromospheric Moreton wave originally predicted by Uchida (1968). The 
first thorough analyses of the coronal bright front were made by Thompson et al. 
(1998) and Thompson et al. (1999), who studied events from 1997 May 12 and 
1997 April 7 respectively. These papers were the first to document in detail 
the different observational features associated with the CBF event, including the 
formation of dimming regions close to the event origin, a post-eruption arcade 
and the bright wavefront which propagated quasi-radially away from the source 
region. The pulses are noted to be quite diffuse, anisotropic, and inconsistent 
with the propagation of a shock, while the estimated velocities are quite low 
(-250 km s" 1 ). 
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These observations of the so-called "EIT wave" did not match the expected 
kinematics of a fast-mode wave that could explain the Moreton wave. In particu- 
lar, the low velocities (typical CBF velocities are ~2-3 times lower than that of a 
Moreton wave) and morphology (CBFs are observed to be broad and diffuse with 
Moreton waves thin and pronounced) are quite different. Despite these disagree- 
ments, Moreton waves and CBFs do show some similarities, with Veronig et al 
(2006) outlining evidence of refraction at a coronal hole boundary for a Moreton 
wave and an associated CBF. 

A detailed examination of a CBF pulse was made by Wills-Davey & Thomp- 
son (1999) using observations of an event made by TRACE, which allowed a 
much higher cadence, higher resolution observation to be made using a smaller 
field-of-view in both the 171 and 195 A passbands. The pulse was observed to 
propagate non-uniformly, with variations in both trajectory and velocity along 
the front of the pulse, indicating Alfvenic behaviour with a strong compressional 
component. An intensity slice across the pulse revealed a Gaussian form, indi- 
cating a single compressional front, while the intensity ratio between the 171 and 
195 A passbands was used to show the presence of plasma heating and constrain 
the temperature through which the pulse propagated to ~1— 1.4 MK. 

The wave interpretation of CBFs was supported by an analysis of multiple 
events performed by Warmuth et al. (2004a,b), who found that the kinematics of 
CBFs lay on a similar kinematical curve to those of a Moreton wave. This result 
strongly implied that the observed CBF could be directly related to the Moreton 
wave and is supported by observations of CBFs and Moreton waves traveling 
co-spatially (Thompson et al, 2000; Okamoto et al, 2004). 

The wave interpretation of CBFs has also been supported by observations 
of reflection and refraction (Long et al, 2008; Veronig et al, 2008; Gopalswamy 
et al, 2009, see Figure 2.4) at active region and coronal hole boundaries. The 
kinematics of CBFs have also been determined to a high degree of accuracy using 
the higher cadence available from the STEREO spacecraft (cf., Long et al, 2008; 
Veronig et al, 2008; Kienreich et al, 2011; Long et al, 2011a), with the derived 
kinematics showing a decelerating pulse that also appears to exhibit dispersion 
with propagation. These observations show good agreement with simulations of 
fast-mode wave propagation (Wang, 2000; Ofman & Thompson, 2002) and are 
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Figure 2.4: Distance-time stack plot showing the variation in intensity along a 
defined line from the erupting active region (N01W08; bottom) towards a nearby 
coronal hole (S10W45; top). The propagating CBF pulse is indicated, along with 
the reflected pulse at ~13:00 UT (Gopalswamy et ai, 2009). 



also consistent with the simulated propagation of a fast-mode wave through a 
randomly structured magnetic field (Murawski et ai, 2001). 

A slight modification to the fast-mode interpretation has been proposed by 
Grechnev et al. (2008), who modeled a CBF event using self-similar dimensional 
analysis (cf. Taylor, 1950a,b; Sedov, 1959). They considered a strong point-like 
explosion in a region of constant density and a region with a radial density drop- 
off from the source (i.e., p oc r~ a ). For a region of constant density, the energy 
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of a strong shock can be determined using dimensional analysis to be 



E oc Po rV, (2.1) 

where r is the radius of the shock, v is the shock velocity and po is the density. 
This produces a shock velocity, 

v oc (—A oc r- 3 / 2 , (2.2) 

with r oc t 2 / 5 . For a region with a radial density fall-off p = br~ a , the velocity v 
is given by, 

uocr [-(3-«)/2], (2.3) 

where r oc £[ 2 /(5-<*)]_ As a result, a strong spherical shock should show deceleration 
for a < 3 and acceleration for a > 3. 

Grechnev et al. (2008) derived kinematics for an event from 2004 July 13 
(see Figure 2.5), arguing that the CBF could be best explained as a coronal blast 
wave (i.e., a shock wave with a low Alfvenic Mach number resulting from a sudden 
release of energy in the corona). This event was quite distinct, with the erupting 
CME breaking apart and some material falling back onto the Sun. A Moreton 
wave and a CBF were both observed for this event, with the authors suggesting 
that their derived kinematics are consistent with an initially strong or moderate 
intensity shock that later damps to a moderate wave pulse. However, Grechnev 
et al. (2008) do note that this mechanism may not be universally consistent across 
all observations of CBFs. 

The fast-mode wave solution has been decried as insufficient by multiple au- 
thors. Wills-Davey et al. (2007) proposed that the physical properties of CBFs 
are more consistent with an MHD soliton than a fast-mode MHD wave. In par- 
ticular, they claim that the single pulse nature, stable morphology and varying 
kinematics can be best explained by a soliton approximation (cf., Russell, 1844). 
As the velocity of a soliton is dependent on the amplitude of the pulse, this would 
account for the large range of derived CBF velocities (Thompson & Myers, 2009), 
while the inherent physical cohesiveness of a soliton would allow it to propagate 
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Figure 2.5: Distance-time plot with derived kinematics for a Moreton wave and 
CBF event from Grechnev et al. (2008). Panel a shows positions of Moreton 
(triangles) and "EIT wave" (crosses) measured in the same direction from the 
eruption source. Panel b shows the derived kinematics for constant p (dotted line) 
and p oc r~ 2 (solid line). 



across large regions of the solar atmosphere relatively coherently. 

This interpretation has failed to gain much traction in the literature, as more 
recent results suggest that CBF pulses may not be as stable as previously thought, 
with clear evidence of dispersion noted by multiple authors (Veronig et al, 2010; 
Kienreich et al., 2011; Long et al, 2011a; Muhr et al, 2011). Widespread pulse 
deceleration has also been noted (e.g., Warmuth et al, 2004a,b; Long et al, 
2011a) that may help to account for the wide range of derived pulse velocities. 

A number of simulations have also been performed indicating that CBFs can 
be best described as fast-mode MHD waves. In particular, the simulations of 
Wang (2000) and Wu et al. (2001) produce pulses that propagate across the quiet 
Sun, avoiding coronal holes and active regions. In both photospheric 
magnetogram was used to extrapolate a realistic coronal magnetic field through 
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which the pulse was allowed to propagate. However, both sets of simulations 
have been criticised for using plasma-/? values that were unrealistic (Chen et al, 
2002). A more recent simulation by Schmidt & Ofman (2010) used the BATS-R- 
US simulation software combined with photospheric magnetograms to simulate an 
event from 2007 May 19 which was observed by the STEREO spacecraft. In this 
case, the simulation is in good agreement with the observed eruption, producing 
a CME and a CBF pulse, which exhibited reflection and refraction from coronal 
hole and active region boundaries similar to the actual event. The advancement 
of simulation software combined with the higher spatial and temporal resolution 
data available from the SDO spacecraft will allow greatly improved simulations 
of CBF events. 

2.2.2 Pseudo-wave Theories 

Several authors have raised issues with the interpretation of CBFs as waves, 
instead choosing to explain them in terms of the magnetic field restructuring 
during the eruption of a coronal mass ejection. This is supported by the results 
of Biesecker et al. (2002), who found that every CBF is associated with a CME, 
although the converse is not necessarily true. A number of observations incon- 
sistent with the behaviour of a magnetoacoustic wave have also been reported, 
including the presence of stationary bright fronts (Delannee & Aulanier, 1999) 
and a dependence of the rotation of the CBF on the helicity of the erupting 
active region (Attrill et al, 2007a). Numerical simulations of CME eruptions 
have also produced evidence of co-existing CBFs and Moreton waves, where the 
propagation of increased emission due to successive opening magnetic field lines 
is observed as the CBF (Chen et al, 2002). 

The first non-wave theory for the interpretation of CBFs was promoted by De- 
lannee & Aulanier (1999). This resulted from the analysis of an event observed 
by SOHO /EFT to contain both moving and stationary bright fronts. Here, the 
CBF was explained as plasma compression near the footpoints of opening field 
lines adjacent to a separatrix during the eruption of a CME. As the CME erupts 
and field lines open farther from the source, the CBF would be expected to prop- 
agate, but crucially, would stop and remain stationary when faced with regions 



41 




Figure 2.6: Simulation illustrating the current shell CBF model proposed by 
Delannee et al. (2008). Left-hand image is a top-down view, right hand image is a 
side-on view. The current shell is visible as the yellow transparent surface. 

containing vertical field lines (such as loop footpoints and coronal holes). 

The theory was readdressed and extended in further analyses (Delannee et al, 
2007, 2008), with the bright CBF ascribed to Joule heating at the edges of a 
current shell containing the erupting CME flux rope (see Figure 2.6). A 3- 
dimensional MHD code was used to simulate the CME eruption corresponding 
to two well-observed CBF events (from 1997 April 7 and 1997 May 12), with 
a twisted flux rope in a bipolar magnetic field configuration used as the initial 
condition. This flux rope was then allowed to erupt, compressing the surrounding 
potential field and producing return electric currents that evolve into a large-scale 
current shell. An integration along this current shell produces a structure similar 
in appearance to the observed CBF pulse. 

There are some issues with this interpretation for CBFs (cf. Patsourakos et al., 
2009). The current shell model implies that the bright front observed as the 
CBF should form at the edges of the erupting CME. This would produce succes- 
sively higher emission as the CME propagates out into the heliosphere. However, 
quadrature observations using the STEREO spacecraft have indicated that the 
CBF emission is confined to ~90 Mm above the photosphere (Kienreich et al., 
2009; Patsourakos & Vourlidas, 2009), a region consistent with the scale height 
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of observed EUV emission. In spite of this criticism, the current shell model 
has been invoked to explain multiple CBF events. Most recently, Schrijver et al. 
(2011) explained a CBF/CME/flare event from 2011 February 15 using the cur- 
rent shell model. Here, the CBF is formed by plasma compression and Joule 
heating at the interface between the erupting CME and an overarching helmet 
streamer structure. The observed emission signatures and derived temperatures 
are also noted to show good correspondence with the predicted behaviour of the 
current shell model. 

An alternative to the current shell model was proposed by Chen et al. (2002) 
based on the results of numerical simulations of a CME eruption. Chen et al. 
(2002) found evidence of two disturbances that propagated laterally away from 
the source of the eruption as the CME flux rope expanded out into the heliosphere. 
As the CME flux rope rises, it induces a piston-driven shock forward of the flux 
rope. The legs of this shock extend to and sweep through the low corona and 
chromosphere, prompting Chen et al. (2002) to interpret this as the manifestation 
of a Moreton wave. Behind these shock legs, a region of plasma compression 
caused by successive opening of field lines covering the flux rope is observed 
propagating laterally away from the source at a velocity three times lower than 
the shock front. It is suggested that this moving compression front is the observed 
"EIT wave" , with Chen et al. (2002) pointing to observed discrepancies in derived 
kinematics for Moreton waves and CBFs for confirmation (Smith & Ramsey, 1964; 
Klassen et al, 2000). 

This model is illustrated in Figure 2.7, with the top part showing several steps 
of the simulation and the bottom part showing the variation in density (p) with 
distance (x) for successive times along a cut of constant height. The piston-driven 
shock is visible in the upper part of Figure 2.7, while the CBF is apparent in the 
lower panel as an increase in density following behind the Moreton wave. The 
theory was extended in work by Chen et al. (2005) and Chen & Fang (2005) 
who used MHD simulations to account for the observational features of CBFs, 
including stationary bright fronts and higher velocities in the quiet Sun compared 
to active regions (Foley et al, 2003). The piston-driven shock model outlined by 
Chen et al. (2005) predicts that two pulses should be identifiable in coronal data; 
the CBF and a coronal pulse corresponding to the Moreton wave. Stationary 
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Figure 2.7: Top: The Chen model for CBFs. As the flux tube rises a shock front 
is induced in front of it (red in the figure). The legs of this shock front sweep 
the chromosphere and are observed as Moreton waves. Simultaneously, the CBF 
is observed as a brightening as the field lines covering the flux rope open. Bottom: 
Stack plot showing the evolution of density (/?) with propagation distance (Chen 
et al, 2002). 



bright fronts are explained as a consequence of the opening magnetic field lines 
associated with the erupting flux rope coming into contact with a neighbouring 
active region. In this case, the magnetic field lines connect to the active region, 
producing a long-lived bright region observed as a stationary bright front. 
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Figure 2.8: Cartoon illustrating the Attrill CME model for the 1997 April 07 
event. 



While the piston driven shock model does account for the low observed CBF 
velocities and observations of stationary bright fronts, it has been difficult to 
definitively prove or disprove as a result of the relatively low cadence of avail- 
able data. The model predicts two pulses associated with each CBF event, with 
a coronal Moreton wave moving at a velocity that is ~3 times faster than the 
CBF pulse. For CBF velocities of ~200-400 km s -1 this would require a coronal 
Moreton wave with a velocity ~600-1200 km s _1 , too fast to be observed by 
SOHO /EIT. However, only single pulses have been observed by STEREO /EUVI 
and the Atmospheric Imaging Assembly (AIA; Lemen et al, 2011) onboard the 
Solar Dynamics Observatory (SDO; Pesnell et a/., 2011) spacecraft. This is de- 
spite the fact that SDO /AIA would be expected to observe very fast pulses as a 
result of its high cadence (Kienreich et al, 2011; Long et a/., 2011a). A double 
pulse event was apparently observed by Chen & Wu (2011), although the evidence 
remains unconvincing. 

An alternative proposal for the reconfiguration of the magnetic field during a 
CME eruption, resulting in a CBF was put forward by Attrill et al. (2006). As 
the CME erupts, the magnetic field lines undergo successive reconnection with 
adjacent quiet Sun loops, producing bright regions through small-scale chromo- 
spheric evaporation that are observed as the CBF. The concept was extended in 
later work (Attrill et al, 2007a, 2009), allowing stationary bright regions to be 
explained as the result of continuous reconnection between the overlying magnetic 
field lines and the the open field lines of a coronal hole. 
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A cartoon illustrating the Attrill concept is shown in Figure 2.8. The erupt- 
ing flux rope (black dotted line) pushes the overlying magnetic field line (black 
solid line), resulting in magnetic reconnection with adjacent quiet Sun loops and 
chromospheric brightening at points A, B and C either side of the erupting flux 
rope. This process is then repeated, with the overlying magnetic field line (black 
solid line) reconnecting with the grey solid line, producing brightening at points 
A', B' and C and so on. 

The Attrill cartoon is not without its detractors, with Delannee (2009) in par- 
ticular questioning the validity of the successive reconnection mechanism. They 
examined the 1997 May 12 event previously studied by Attrill et al. (2007a) using 
an extrapolated magnetic field topology to examine the interconnectivity of the 
coronal magnetic field. Whereas Attrill et al. (2007a) claimed that the erupting 
flux rope was magnetically linked to adjacent quiet Sun loops, Delannee (2009) 
found that this was not the case, with the erupting active region linked to mag- 
netic field lines ~300 Mm away. There have also been some conflicts between the 
predicted and observed pulse morphology for this model, with quadrature obser- 
vations of CBFs indicating that the emission is produced too high in the solar 
atmosphere to be a result of small-scale chromospheric brightening (Patsourakos 
& Vourlidas, 2009; Kienreich et al, 2009). 

Several authors have proposed that CBFs may be a combination of both 
wave and pseudo-wave modes rather than a single interpretation. This was first 
suggested by Zhukov & Auchere (2004), who studied several CBF events using 
the 171, 195, and 284 A passbands available from SOHO/EYT (in the process 
making the first CBF observation in the 284 A passband). The events studied 
were found to display both wave and non-wave characteristics, with the isotropic 
nature of the pulses implying a wave nature, while the extent of the coronal 
dimming region behind the pulse and the low pulse velocities indicate that this 
interpretation is insufficient. 

Later work also concluded that this scenario may be the optimum solution to 
the CBF problem. Cohen et al. (2009) used the BATS-R-US MHD code to sim- 
ulate the eruption of a coronal mass ejection from 2009 February 13, comparing 
the simulation with observations of the eruption made by STEREO /EUVI. It was 
found that the CBF front was best represented by a density enhancement result- 
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ing from a combination of both wave and non-wave processes. The CME drives 
the front through lateral expansion in the low corona, with this lateral expansion 
resulting in widespread magnetic reconnection as the coronal magnetic field reor- 
ganises. A similar result was obtained by Downs et al. (2011), who used a similar 
model to study the evolution of an event from 2008 March 25. In this case, the 
signatures of the different EUV passbands available from STEREO /EUVI were 
synthesised, allowing the simulation to be directly compared to observations. 
Two distinct signatures were observed, one identified as a fast-mode wave pulse, 
the other a separate compression front produced by the expansion of the CME 
into the surrounding magnetic structures. Much work remains to be done in this 
regard, and the high cadence data now available from SDO/AIA will allow this 
issue to be resolved. 
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Chapter 3 



Instrumentation 



These instruments have play'd me so many tricks that I have at last found them 
out in many of their humours. 



Sir William Herschel 



In this chapter, the different telescopes used to study "EIT Waves" are dis- 
cussed, with special emphasis placed on the software required to process the raw 
data and prepare it for scientific analysis. The host spacecraft for each instru- 
ment have unique capabilities that are also discussed as these can influence how 
the data is approached. 

Although there have been some indications of large-scale disturbances in ad- 
ditional portions of the electromagnetic spectra (e.g., radio and X-rays), the dis- 
turbances studied in this project were observed using the Extreme Ultra Violet 
(EUV) part of the solar spectrum. The Earth's protective ozone layer makes 
this difficult as EUV radiation does not reach the ground, making it necessary 
to use space-based observatories to study the low corona. The two observatories 
primarily used in this work were the STEREO and SDO spacecraft. 
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3.1 The STEREO /Extreme UltraViolet Imager 
(EUVI) 

The STEREO mission consists of two identical spacecraft, one following the Earth 
(STERE 'O-Behind) and the other leading the Earth (STEREO -Ahead) on its 
orbit around the Sun. STEREO-k has an orbit slightly inside that of the Earth, 
with STEREO-B following an orbit slightly outside that of the Earth. As a 
result, both spacecraft are separating from the Earth and each other at rates of 
~22 and ~44 degrees per year, respectively. This configuration was chosen to 
allow stereoscopic imaging of the Sun and the Sun-Earth line, thus allowing a 
determination of the feasibility of possible future missions. They were launched 
on the same Delta II rocket in October 2006, with science operations beginning 
in December 2006 after a series of maneuvers allowed both spacecraft to escape 
Earth orbit. Quadrature was reached in February 2009, with both spacecraft at 
90 degrees to each other, while opposition was achieved in February 2011 when 
the entire Sun was simultaneously visible for the first time. 

The EUVI telescope is part of the Sun Earth Connection Coronal and Helio- 
spheric Investigation (SECCHI; Howard et al, 2008) instrument suite onboard 
STEREO. SECCHI consists of an EUV imager (EUVI), two coronagraphs with 
different fields-of-view (1.4-4 R for COR-1 and 2.5-15 R for COR-2) and two 
heliospheric imagers observing the Sun-Earth line (at distances of 15-80 R for 
HI-1 and 80-215 R for HI-2) with the fields-of-view of each instrument shown 
in Figure 3.1 for comparison. This allows solar eruptions to be tracked from ini- 
tiation on the Sun through the corona and heliosphere to the Earth. EUVI has 
a similar design to the EIT instrument onboard the SOHO spacecraft, but with 
a higher spatial and temporal resolution. 

The increasing separation of the STEREO spacecraft from each other and 
Earth has resulted in a gradual decrease in the amount of data being sent from 
the spacecraft to Earth and in the observing cadence of the different passbands. 
The scientific emphasis of STEREO has also changed with the increased distance 
from Earth, with the initial focus being on the initiation of CMEs as observed 
using EUVI. The advent of SDO and the present ability of STEREO to observe 
the Sun-Earth line has resulted in an increased emphasis on CME observations 
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Figure 3.1: The field-of-view of the SECCHI instrument suite as observed by 
STEREO-A. Top panel shows EUVI and the coronagraphs Cor-1 and Cor-2. Bot- 
tom panel shows the field of view of the coronagraphs compared to the Heliospheric 
Imagers HI-1 and HI-2 (cf. Howard et al, 2008). 



using the coronagraphs and heliospheric imagers. 
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Figure 3.2: EUVI images showing the 195 A (top left), 171 A (top right), 284 A 
(bottom left) and 304 A (bottom right) passbands. Each passband has been plotted 
with the appropriate colour table for clarity. 

3.1.1 EUVI Instrument Operation 

EUVI is a normal incidence telescope with a Ritchey-Chretien construction, de- 
signed to allow a large field-of-view within a compact space by utilising hyper- 
bolic primary and secondary mirrors. This layout allows an unvignetted, circular 
full-Sun field-of-view extending out to 1.7 R without any requirement for a 
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Figure 3.3: Cross-section of the EUVI telescope onboard the STEREO spacecraft 
(Wuelser et al, 2004). 



focusing mechanism. The telescope images four separate EUV passbands peak- 
ing at 304 A, 171 A, 195 A and 284 A, with an effective temperature range from 
~60 000 K (in the 304 A passband) to ~2.5 MK (in the 284 A passband). Images 
from the four passbands observed by EUVI are shown in Figure 3.2. The different 
panels have been coloured using the appropriate colour tables for clarity. 

The actual operation of the telescope is as follows; light enters the telescope 
at the front (the right of Figure 3.3), where thin metal film filters (aluminium- 
on-polyimide on a coarse nickel grid for 171 and 195 A and aluminium on a fine 
nickel mesh for 284 and 304 A) remove undesired UV, IR and visible radiation, 
and keep heat out of the telescope. The required passband is selected through 
the use of a rotating quadrant selector allowing the radiation to pass through to 
one of the four quadrants of the optics. These are specially coated using narrow- 
band, multilayer MoSi coating of variable thickness to ensure optimal reflectivity 
for each passband, producing peak reflectivity values of ~ 15-39 %. While this 
may appear low, the very high intensities of the observed lines mean that this 
does not greatly affect performance. The telescope is also fully baffled to negate 
the effects of scattered light and charged particles entering the front aperture. 

The CCD used by the EUVI telescopes is a backside-thinned, back-illuminated 
full-frame CCD. This arrangement allows improved low-light performance and 
short wavelength response. The CCD is also passively cooled to below — 60°C 
(~213 K) through the use of an aluminium "cool finger" to reduce dark current. 
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Figure 3.4: Peak-normalised response curves versus temperature for each EUV 
passband observed by STEREO /EUVI. Peak emission temperatures of the different 
passbands are given in the upper right for comparison (cf. Wuelser et al, 2004) 



EUVI has a spatial resolution of 1.6 arcsec per pixel with an initial passband- 
dependent observing cadence of 75-300 s for 171 A, 300-600 s for 195 A and 
304 A and 1200 s for 284 A. However, as the lifetime of the mission has continued 
and the spatial separation of the spacecraft has increased, the observing cadence 
has been greatly reduced to allow for decreased bandwidth. Currently, EUVI 
operates at a typical cadence of 300 and 600 s in the 195 A and 304 A passbands 
respectively, with both the 171 A and 284 A passbands operating synoptically 
(i.e., one image every ~4 hours). 

3.1.2 EUVI Temperature Response 

The four passbands observed by EUVI cover a broad range of temperatures and 
hence different regions of the solar atmosphere. The peak-normalised tempera- 
ture response curves for each of the four passbands are shown in Figure 3.4, while 
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Table 3.1. Temperature response characteristics of EUVI passbands 



Passband (A) Dominant ion T pea k (MK) 



304 He II, Si XI a 0.08 

171 Fe IX, Fe X 0.9 

195 Fe XII 1.4 

284 Fe XV 2.0 



a Although the hydrogen-like He II is the 
dominant ion here, there is a small contribu- 
tion from the coronal Si XI emission line. 



the peak emission temperatures (T pea k) and dominant ions of the individual pass- 
bands are outlined in Table 3.1. 



The range of values for T pea k shown in Table 3.1 show that EUVI is well- 
suited for studying the relatively cool quiet corona, while the 284 A passband 
allows an analysis of the hotter active regions. This is apparent in the different 
panels in Figure 3.2, where the cool 304 A passband shows the structure of the 
upper chromosphere, with cool structure visible in the 171 A passband and hotter 
material apparent in the hotter 195 A and 284 A passbands. 

Although the individual passbands observed by EUVI have been tuned using 
multilayer coatings and filters, the temperature response curves remain quite 
broad with significant overlap between passbands as shown in Figure 3.4. This 
means that EUVI tends to observe plasma at a typical temperature of 1-3 MK, 
consistent with the low solar corona and particularly the quiet Sun. This makes 
EUVI a valuable tool for studying the initiation of CMEs, in line with the original 
science objectives. 
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3.1.3 EUVI Data Preparation Software 

Before the data from EUVI can be used for scientific analysis, the raw FITS 1 files 
must be calibrated and corrected for instrumental and other effects. The raw 
FITS file (classified as a Level 0.5 FITS file) is treated using the secchLprep.pro 
calibration routine contained within the SolarSoftWare (SSW; Freeland & Handy, 
1998) software package. 

The secchLprep.pro routine is designed to treat data from the SECCHI instru- 
ment suite onboard STEREO, calibrating the data according to the instrument 
used. Once the EUVI image has been positively identified, the calibration involves 
the application of several corrections to the data to account for instrumental and 
downlink effects. Figure 3.5 shows the same image from 2007 May 19 at 12:42 UT 
in the 171 A passband for different stages of the calibration process, with panel 
(a) showing the original raw image as taken by STEREO-k. 

The first step in the calibration process involves identifying and removing cos- 
mic rays hits and missing or corrupted pixels from the image which can be pro- 
duced by hardware failure or transmission losses. Next, the image is normalised 
to account for the different passband-specific filters on the telescope, with the 
resulting image shown in panel (b). The next step involves identifying the image 
bias value defined in the file header and subtracting it from the image (panel c). 
Panel (d) shows the result of normalising the image with respect to the exposure 
time. This is a consequence of the different exposure times available for each 
passband and also for certain observing programs. Next, the effects of vignetting 
are removed using a calibration image and a photometric correction applied to 
the image values using a calibration factor; this is shown in panel (e). The fi- 
nal step involves correcting the image for any image processing that may have 
been carried out onboard the spacecraft and ensuring that solar north is aligned 
towards the top of the image. The image can then be scaled and displayed as 
shown in panel (f). 



1 Flexible Image Transport System; a filctype typically used in astronomy that includes an 
informational header in addition to the observed data. 
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All of these calibration processes are applied to the data to produce a scientific 
quality Level 1.0 FITS file. The data may then be processed further as required, 
or plotted (as shown in Figure 3.2). Each of these processes is automatically 
applied, although it is possible to process the data without applying one or more 
of these steps using keywords. This is often required when analysing temperature 
or emission measures. 

3.2 The SDO /Atmospheric Imaging Assembly 
(AIA) 

The SDO spacecraft (launched in 2010 February, with science operations begin- 
ning in 2010 April) is the first part of NASA's "Living with a Star" program, 
designed to study the effects of the Sun on the near-Earth environment. SDO is 
currently positioned in a geosynchronous orbit located above a dedicated ground 
station in White Sands, New Mexico. This is a necessary consequence of the very 
large data output of SDO; ~ 2 terabytes of data per day, every day. The largest 
contribution to this data flow comes from the Atmospheric Imaging Assembly 
(AIA) instrument, which uses four separate telescopes to take one 4096x4096 
pixel image in ten different passbands at a temporal cadence of ~ 12 s per pass- 
band. Eight of the ten different passbands available from SDO /AIA are shown 
in Figure 3.6. 

3.2.1 AIA Instrument Operation 

The Atmospheric Imaging Assembly instrument (Lemen et al, 2011) consists 
of four Cassegrain telescopes designed to observe the low solar atmosphere in 
both UV and EUV emission. The telescopes have an F-number f/20 1 with a 
20 cm primary mirror and a secondary mirror incorporating three piezoelectric 
transducers that offer tip-tilt stabilisation. Metal filters (aluminium for the 171 A 
and longer wavelength channels and zirconium for the shorter passbands) prevent 
stray EM radiation from entering the telescope, while each telescope is fully 

1 e.g., a focal length that is 20 x the pupil diameter 
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Figure 3.6: Panels showing some of the different passbands available from 
SDO/AIA for ~06:15:00 UT on 2011 Jun 7. Each passband has been plotted 
with the appropriate colour table for clarity, with the passbands identified in the 
top left corner of each panel. 

baffled to protect against scattered light. The metal filters are hosted on a 70 line- 
per-inch nickel mesh, which can result in diffraction patterns from high intensity 
point source emission (i.e., flares). 

A cross-section of one of the AIA telescopes is shown in Figure 3.7. Radiation 
enters the telescope through the entrance aperture on the left of the figure and 
passes via the primary mirror and active secondary mirror to the shutter at the 
right of the figure. The shutter is used to control the exposure time of the image; 
this varies according to passband but is of the order ~2 s. AIA includes an 
Automatic Exposure Control (AEC) system, which automatically modifies the 
exposure time if the intensity of the image increases dramatically (as a result of 
e.g., a flare). In this case, the exposure time can drop to ~0.2 s. Each telescope 
in AIA is designed to study two distinct passbands, with the desired wavelength 
channel chosen using a filter wheel. 

The radiation then passes through to the CCD (which is identical for each 
telescope). Each CCD contains 4096x4096 pixels, and is back-thinned and back- 
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Figure 3.7: Cross-section of one of the four AIA telescopes onboard the SDO 
spacecraft (Lemen et at, 2011). 

illuminated to improve signal. The telescope has a full-Sun spatial resolution out 
to ~0.5 R above the solar limb of 0.6 arcsec per pixel, superior to SOHO/EYT, 
TRACE and STEREO /EUVI. 

AIA is designed to operate at a continuous cadence of 12 s in each passband, 
although it is possible to increase the cadence to 10 s if necessary. This com- 
bination of high temporal cadence and large image size produces a data volume 
that exceeds the telemetry allowance, so a lossless Rice compression algorithm is 
applied to the data onboard the spacecraft to ensure continuous data supply. 

3.2.2 AIA Temperature Response 

The large number of passbands observed by AIA allow it to continuously monitor 
a large part of the solar atmosphere and a wide range of phenomena. Table 3.2 
outlines the different passbands studied along with the dominant ions and peak 
emission temperatures in each case. 

STEREO /EUVI and SDO /AIA both observe the 304 A, 171 A, and 195/193 A 
passbands, making these passbands useful for confirming observations. However, 
the 171 A and 195/193 A passbands in particular have different temperature re- 
sponse curves, and as a result do not necessarily share observations of the same 
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Figure 3.8: Peak-normalised response curves versus temperature for each EUV 
passband observed by SDO / AIA. The peak emission temperatures of the different 
passbands are given in the upper right for comparison. This figure is a peak- 
normalised version of that from Lemen et al. (2011). 



phenomenon. This is most apparent when studying "EIT waves", as the pulse 
appears in emission in EUVI 171 A, while simultaneously appearing in absorp- 
tion in AIA 171 A. This observational discrepancy is discussed in more detail in 
Section 7.3.3, along with an analysis of how this affects the interpretation of the 
pulse. 

Although AIA observes multiple passbands using highly tuned niters, the 
broadband nature of the imager is apparent from the temperature response curves 
shown in Figure 3.8. As a result, determining temperature and differential emis- 
sion measure estimates from AIA images is subject to misinterpretation; these 
issues are discussed in Section 7.3.3. 
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Table 3.2. Temperature response characteristics of AIA passbands. 



Passband (A) 


Dominant ion 


log T peak (K) 


UV/ Visible 


Continuum, C IV a 


3.7, 5.0 


304 


He II 


4.7 


171 


Fe IX 


5.9 


193 


Fe XII, Fe XXIV 


6.2, 7.3 


211 


Fe XIV 


6.3 


335 


Fe XVI 


6.4 


94 


Fe XVIII 


6.8 


131 


Fe VIII, Fe XXI 


5.6, 7.0 



a Visible and UV light is observed in the 1600, 
1700, and 4500 A passbands using one half of one 
telescope. 



3.2.3 AIA Data Preparation Software 

Once downloaded to the base station at White Sands in New Mexico, all SDO 
data is sent to the Joint SDO Operations Centre (JSOC), located at Stanford 
University in California where the data is stored permanently. Here, a processing 
pipeline is used to convert the raw Level data to fully-calibrated Level 1 data. 
This involves several distinct processes including, 

• The removal of over and under scan regions of the image 

• The removal of dark current effects 

• The application of a flat-field correction 

• The automated removal of bad pixels and cosmic ray hits 

• Flipping the image to place solar north at the top 

Once this pipeline processing has been completed, the Level 1 FITS files are 
sent to a number of data hosting stations worldwide (including the Harvard- 
Smithsonian Centre for Astrophysics in the U.S.A., the University of Central 
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Lancashire in the U.K. and the Royal Observatory of Belgium) where it is made 
publicly available. 

Once downloaded by the user, the Level 1 FITS files can be further processed 
using the aia.prep.pro processing routine. This applies additional corrections 
to the data that allow an intercomparison of the different passbands. This in- 
volves (i) co-aligning the images to correct for residual offsets between the four 
telescopes; (ii) adjusting the plate-scale to correct for fractional offsets in focal 
length; and (iii) co-aligning the bore-sight pointing of the different telescopes. 
With these adjustments applied to the data, the FITS header is adjusted to 
Level 1.5, and the images are now ready for scientific analysis. 
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Chapter 4 



The Kinematics of a 
Globally-propagating 
Disturbance in the Low Corona 



In natural science the principles of truth ought to be confirmed by observation. 

Carolus Linnaeus 



In this chapter, the first observations of a CBF made using data from the 
EUVI instrument onboard the STEREO spacecraft are used to determine the 
kinematics of the pulse. It is shown for the first time that a CBF has similar 
kinematics in all four EUV passbands studied (304, 171, 195, and 284 A). In the 
304 A passband the disturbance shows a velocity peak of 238±20 km s _1 within 
~28 minutes of its launch, varying in acceleration from 76 m s~ 2 to -102 m s~ 2 . 
This passband contains a strong contribution from a Si XI line (303.32 A) with 
a peak formation temperature of ~1.6 MK, suggesting that the 304 A emission 
may be coronal rather than chromospheric in origin. Comparable velocities and 
accelerations are found in the coronal 195 A passband, while lower values are 
found in the lower cadence 284 A passband. In the higher cadence 171 A passband 
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the velocity varies significantly, peaking at 475±47 km s _1 within ~20 minutes 
of launch, with a variation in acceleration from 816 m s~ 2 to -413 m s~ 2 . The 
high image cadence of the 171 A passband (2.5 minutes compared to 10 minutes 
for the similar temperature response 195 A passband) is found to have a major 
effect on the measured velocity and acceleration of the pulse, which increase 
by factors of ~2 and ~10, respectively. This implies that previously measured 
values (e.g., using EIT) may have been underestimated. It is also noted that 
the disturbance shows strong reflection from a coronal hole in both the 171 and 
195 A passbands. These observations are consistent with an impulsively generated 
fast-mode magnetoacoustic wave. The work presented in this chapter is based 
upon results published in The Astrophysical Journal Letters by Long, Gallagher, 
McAteer & Bloomfield (2008). 

4.1 Introduction 

The launch of the STEREO spacecraft, with its EUVI instruments signalled the 
beginning of a new age of solar physics. As noted in Section 3.1, STEREO 
consists of two separate spacecraft, designed to study the inner heliosphere from 
ahead and behind the Earth on its passage around the Sun. This represents the 
first time that the Sun can be observed stereoscopically, potentially allowing 3-D 
modeling of the inner Heliosphere. 

As well as offering two distinct locations for joint observations of the Sun and 
the inner Heliosphere, the suite of imaging instruments onboard STEREO repre- 
sent a major improvement compared to SOHO. Whereas SOHO carried an EUV 
imager (EIT) and several coronagraphs (LASCO C-l, -2 and -3), STEREO car- 
ries the EUVI telescope, two coronagraphs (Cor-1 and -2) and two Heliospheric 
Imagers (HI-1 and -2) that allow direct imaging of the inner Heliosphere. Each 
of these instruments also operate at a higher cadence to those onboard SOHO, 
while EUVI observes all four passbands at a high temporal cadence, unlike the 
synoptic operations of SOHO /EIT. These instruments were designed specifically 
to allow a deeper understanding of CMEs and the on-disk phenomena associated 
with them. 

However, the launch of the STEREO spacecraft in December 2006 coincided 
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with the beginning of an unexpectedly long period of low solar activity. As a 
result, the first major eruption on the Sun did not occur until ~6 months after 
launch, in May 2007, by which time the STEREO spacecraft were ~8 degrees 
apart. This small separation meant that events occurring close to disk centre 
as observed from the Earth could be observed by both spacecraft with minimal 
projection effects. 

The design of the STEREO spacecraft and the EUVI telescopes make them 
ideal for studying the on-disk propagation of CBF pulses, particularly early in 
the mission when the imaging cadence is high and both spacecraft can observe the 
same event. Traditional analysis of CBFs has been hampered by the low imaging 
cadence available from SOHO/EYT, which produced two or three observations per 
event. This makes it difficult to derive pulse kinematics, and forces a combination 
of the EUV observations with measurements from other passbands. While these 
combinations are vital for a full understanding of the pulse, the sparse nature of 
the data complicates interpretation. The true kinematics of a CBF pulse must 
therefore be determined using purely EUV observations, a situation heretofore 
impossible. 

These difficulties in deriving the kinematics of CBF pulses have resulted in a 
number of alternative interpretations of this phenomenon. Although magnetohy- 
drodynamic wave theory (see Section 1.3) can explain many of the brightenings 
observed in EUV images as well as explaining some of the characteristics of EIT 
waves, including their slow expansion rates (Wang, 2000), it is inconsistent with 
the low velocities derived from SOHO/EYT observations (see Section 2.2.1 for 
more details). However, there is some indication that CBFs exhibit wave prop- 
erties, with Veronig et al. (2006) noting signs of pulse refraction at the boundary 
of a coronal hole. 

The inconsistency between the observed pulse behaviour and that predicted 
by MHD wave theory has resulted in some authors promoting a non-wave in- 
terpretation for CBFs. Delannee (2000), Chen et al. (2005) and Attrill et al. 
(2007b) have each proposed models where the observed CBF is not a true wave, 
but rather a byproduct of the magnetic field restructuring as the associated CME 
erupts away from the Sun. The specific details of each model differs slightly, with 
Delannee (2000) proposing heating in the separatrix between magnetic field lines 
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as the CME erupts, an idea similar to that of Chen et al. (2005) who proposed 
that the reconnecting magnetic field lines either side of the erupting CME pro- 
duced the brightening observed as the CBF. In contrast, Attrill et al. (2007b) 
proposed that the reconnection of the erupting CME with quiet Sun loops either 
side of the CME resulted in small-scale brightenings seen as the CBF pulse. 

Although the sets of both wave and pseudo-wave models are formed in differ- 
ent ways, they were proposed based on observation, with the result that it is dif- 
ficult to distinguish them by observation alone. However, there are discrepancies 
in the predicted behaviour of the theoretical pulses. In particular, the kinematics 
and morphology of the pulses should provide an indication of the physical prop- 
erties of the pulse. For the first time, the high cadence of STEREO /EUVI allow 
the pulse kinematics to be determined to a high degree of accuracy using purely 
EUV observations, while the dual viewpoints available from STEREO provide an 
independent confirmation of the kinematics. 

In this chapter, an analysis of the CBF event from 2007 May 19 as observed 
by the EUVI instrument onboard both STEREO spacecraft is presented. The 
observations made using the STEREO spacecraft are outlined in Section 4.2, with 
the methods used to process the data discussed in Section 4.3. The data are then 
analysed, with the results of this analysis presented in Section 4.4. These results 
are discussed and some conclusions drawn in Section 4.5. Finally, Section 4.6 
discusses the results of this work in the context of follow-up work performed by 
other authors. 

4.2 Observations 

A GOES class B9 flare erupting from active region NOAA 10956 was observed 
on 2007 May 19, beginning at 12:34 UT, before peaking at 13:02 UT and then 
ending at 13:19 UT. A CME was initially observed in the LASCO C-2 field-of- 
view at 13:24 UT, with the SOHO /LASCO CME catalogue classifying it as being 
associated with this eruption. A filament lift-off was also observed to begin at 
~12:31 UT, slightly before the first observation of the CBF pulse. The eruption 
was observed at high cadence in multiple EUV passbands by both STEREO 
spacecraft, although SOHO /EIT was undergoing CCD bake-out (to remove the 
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Figure 4.1: Evolution of the disturbance with time in all four EUVI wavelengths 
from STEREO-A. Wavelengths are arranged from left to right as 304 A, 171 A, 
195 A and 284 A. 

effects of radiation damage) at the time of the eruption and was not available. 

As the first major solar eruption to be observed by STEREO at high cadence, 
multiple aspects of this event have been studied in detail by other authors. In 
particular, the stereoscopic capabilities of STEREO allowed a detailed examina- 
tion of the filament eruption (Bone et al., 2009; Liewer et al., 2009), while the 
location of the erupting active region near disk centre made it possible to analyse 
the evolution of the magnetic field (Li et al., 2008; Conlon et al, 2010). The high 
cadence observations of the erupting CBF and the nature of its free propagation 
also meant that the CBF was further analysed by additional authors including 
Veronig et al. (2008); Gopalswamy et al. (2009) and Attrill (2010). 

Although the imaging cadence of STEREO /EUVI has shown some variation 
over the timescale of the mission, this event occurred soon after launch. As a 
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result, the cadence was very high, at 150 s for the 171 A passband, 600 s for the 
195 and 304 A passbands and 1200 s for the 284 A passband. The full-disk images 
were also analysed using the full spatial resolution of the instrument (1.6 arcsec 
per pixel) to ensure a thorough analysis. 

Figure 4.1 shows the CBF pulse observed using data from STEREO-A and 
highlighted using the techniques outlined in Section 4.3. The passbands have been 
arranged here from left to right as 304 A, 171 A, 195 A and 284 A, with the times 
of the leading image in each case shown in the top right of each image. These 
images were chosen to show the pulse just after formation, during its propagation 
and before it becomes too diffuse for clear identification. 

4.3 Methods 

Before proceeding with identification and analysis of the pulse, each image used 
was first processed using the SolarSoft calibration routine secchLprep.pro (see 
Section 3.1.3 for details) to remove CCD bias, cosmic ray hits and missing or 
corrupted pixels that could affect the the image quality. The images were then 
flat-fielded to ensure uniform intensity variation and normalised to account for 
the different filters used for each passband. The effects of image processing carried 
out onboard the spacecraft were also accounted for. 

Once these initial processes had been applied to the data, each image was de- 
rotated to the same pre-event time to counter the effects of solar rotation. The 
propagating pulse was then highlighted using running difference images. This is a 
simple image processing technique traditionally used to highlight transient solar 
events (cf. Dere et al, 1997; Thompson et al, 1998; Zhukov & Auchere, 2004), 
that allows the identification of motion by subtracting the previous image in time 
from each image (i.e., Al — I t — I t -i). Each difference image was then smoothed 
using a boxcar filter to highlight the diffuse pulse while retaining detail. The 
width of the boxcar filter applied was defined by the passband, with widths of 
15, 9, 3 and 7 pixels for the 304, 171, 195 and 284 A passbands respectively. 

Once the pulse had been visually identified in each individual image, a he- 
liographic grid was applied to each image to allow the position of the pulse to 
be determined. This is shown in Figure 4.2 for both STEREO-B (left) and 
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Figure 4.2: STEREO-B (left) and STEREO-A (right) images from the 195 A 
passband with heliographic grid overlay. The north pole of the grid in each case is 
centered on the flare kernel. 

STEREO-A (right), with the grid oriented in each case so that the north pole 
was located at the flare kernel. The front of the pulse was then identified in sub- 
sequent images using a point-and-click technique with the heliographic grid lines 
used as a guide to ensure that the temporal variation of the pulse position was 
determined along the same great circle line. The distance travelled by the pulse 
could then be determined using the angular separation between the source and 
the identified front position for a given time. The same grid line towards solar 
west was used for each passband and for each spacecraft to ensure continuity. 
The variation in velocity and acceleration of the pulse were then derived using 
the three-point Lagrangian interpolation technique available from the deriv.pro 
routine in IDL (see Section 5.2.1.4 for more details). 

4.4 Results 

The observed CBF pulse is apparent in Figure 4.1, particularly in the second row 
of difference images which clearly shows a bright front followed by a dimming 
region on the solar disk for all four passbands. Although this has been previ- 
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Figure 4.3: Distance-time (top), velocity-time (middle), and acceleration-time 
(bottom) plots of the wavefront. Left-hand panels show data measured by 
STEREO-A, while right-hand panels show STEREO-B. All distances are measured 
to the leading edge of the running difference brightening along the great circle lon- 
gitude line to solar west (see Fig. 4.2). The first and last data point error bars in 
the 171 A data have been divided by two for display purposes. 

ously observed using the 195 A (Thompson et ai, 1998), 171 A (Wills-Davey & 
Thompson, 1999) and 284 A (Zhukov & Auchere, 2004) passbands, this is the 
first observation of the "EIT wave" phenomenon using the 304 A passband. This 
is also the first time that the pulse has been observed in all four EUV passbands 
simultaneously, and as a result is very important for a complete physical analysis 
of the disturbance. As each passband corresponds roughly to a different temper- 
ature and thus a different region of the solar atmosphere, this observation should 
greatly assist with the 3-D modeling of CBFs. 

The measured distance-time data is shown in the top row of Figure 4.3 for each 
passband at its original cadence. The high cadence 171 A measurements appear 
to show a more rapid increase in propagation distance compared to the lower 
cadence 195 A and 304 A measurements. Similarly, the lowest cadence 284 A 
measurements show a more gradual increase in propagation distance compared 
to the 195 A and 304 A passbands. The similarity in the observing cadence of 
the 195 A and 304 A passbands is also reflected in the similar temporal variation 
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Table 4.1. Derived pulse velocity and acceleration estimates. 





Passband 


Cadence 


Vmax 


O'max 


O'min 


Spacecraft 


(A) 


(minutes) 


(km s^ 1 ) 


(m s~ 2 ) 


(m s~ 2 ) 


STEREO-k 


304 


10 


238 


76 


-102 




171 


2.5 


475 


816 


-413 




195 


10 


262 


93 


-109 




284 


20 


153 


23 


-32 


STEREO-^, 


304 


10 


276 


99 


-105 




171 


2.5 


500 


881 


-339 




195 


10 


284 


105 


-117 




284 


20 


158 


23 


-38 



Note. - Mean velocity errors are <50 km s , mean acceleration 
errors <200 m s~ 2 

in propagation distance. 

The difficulty in directly comparing measurements of unequal cadence is high- 
lighted by the derived velocity-time graphs shown in the middle row of Figure 4.3, 
where the variation in velocity appears to be strongly cadence-dependent. The 
curve corresponding to the 171 A images shows a significant peak in velocity (in- 
creasing from rest to 475±47 km s" 1 within approximately 20 minutes), which 
is distinctly different from the more gradual peaks observed in the 195 A, 304 A 
and 284 A velocities. The 195 A and 304 A passbands show similar velocity 
variations; a peak velocity of 262±4 km s" 1 within ~28 minutes for 195 A and 
a velocity peak of 238±20 km s _1 within ~28 minutes in 304 A. However, the 
lowest cadence 284 A passband exhibits the smallest velocity variation, with a 
peak in velocity of 153±5 km s _1 within ~27 minutes of launch. 

The highly variable nature of the kinematics is also apparent in the derived ac- 
celerations, with variations from 816 m s~ 2 to -413 m s~ 2 , 93 m s -2 to -109 m s~ 2 , 
76 m s~ 2 to -102 m s~ 2 and 23 m s~ 2 to -32 m s" 2 for the 171 A, 195 A, 304 A 
and 284 A passbands respectively. The derived kinematics are summarised in 
Table 4.1 for the observations from both STEREO-k and STEREO-B. The 
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variations in velocity and acceleration can be seen in data from both STEREO 
spacecraft, although the peak values differ slightly between the two spacecraft as 
a consequence of the different spacecraft positions. The errors associated with 
the derived velocity and acceleration values are the standard deviation of the 
derivative at each point. As a result, the error bars for the first and last value 
were larger than for the middle values, although the typical velocity errors were 
<50 km s _1 with the typical errors for the derived acceleration <200 m s~ 2 . 

To overcome the difficulties in directly comparing measurements made using 
different image cadences, the cadence of the 171 A passband was degraded from 
2.5 minutes to 10 minutes to match that of the 195 A and 304 A passbands. This 
was done by taking every fourth image from 171 A, with the effective times of 
the images chosen to be as close as possible to that of the 195 A images. The 
284 A data was not examined here as its very low 20 minute cadence meant 
that it was not comparable to the other three passbands. The same running 
difference analysis previously described was then applied to these images, yielding 
the results shown in Figure 4.4. Note that the reduced cadence had the effect of 
increasing the signal of the pulse for the 171 A passband, allowing the pulse to 
be identified after the time at which it could not longer be identified using the 
normal cadence images. This is a result of the longer time-step between images, 
allowing the pulse to travel beyond its previous extent and making it easier to 
track over a larger distance. 

It is apparent that the derived kinematics for the 171 A passband measure- 
ments made at lower cadence are comparable to the 195 A and 304 A mea- 
surements. The peak in velocity for the 10 minute cadence 171 A measure- 
ments is much more gradual than for the 2.5 minute cadence measurements and 
the acceleration-time plot, while still showing acceleration followed by decelera- 
tion, is not as pronounced as in Figure 4.3. The discrepancy between the higher 
2.5 minute and lower 10 minute cadence 171 A measurements suggests that the 
derived kinematics may be influenced by the observing cadence of the instrument. 
Although 2.5 minutes is the highest available cadence for this event, STEREO 
is capable of observing at up to ~75 s cadence, while the forthcoming Solar Dy- 
namics Observatory will have an observing cadence of ~12 s. This will be vital 
for examining this effect in more detail. 
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Figure 4.4: Same as Fig. 4.3, but with 171 A at 10 minute cadence and excluding 
284 A. All data have the same cadence and running difference images have as close 
as possible effective times in each of the three passbands. The first and last data 
point error bars in the 171 A and 304 A data have been divided by two for display 
purposes. 

4.5 Conclusions 

The first major solar eruption with an associated coronal propagating front to 
be observed by the STEREO spacecraft was studied using EUV images available 
from the EUVI instrument. The disturbance was observed in all four passbands 
covered by EUVI, namely the 171 A, 195 A, 284 A and for the first time, the 304 A 
passband. Despite the varying cadences, it was possible to derive the kinematics 
of the pulse in each passband, allowing the temporal variation in distance, velocity 
and acceleration to be examined. 

The kinematics derived from the higher cadence 171 A measurements indi- 
cate an average front velocity of 196 km s" 1 , peaking at 475±47 km s -1 , with 
an initially positive acceleration that peaked at 816±220 m s~ 2 before decreasing 
to -413±320 m s~ 2 ; i.e., it initially accelerates away from the flare point before 
decelerating further away. This peaked velocity and positive to negative varia- 
tion in acceleration was also noted for the other passbands studied, albeit with 
decreased strength. The observed rise and fall in velocity is in contrast to the 
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model proposed by Delannee et al. (2008) which indicates a constantly increasing 
expansion velocity. The peak in velocity derived for the 171 A passband is also at 
the higher end of the previously measured range of "EIT wave" velocities, while 
the peak velocities determined for the lower cadence 195 A and 304 A passbands 
(262±4 km s _1 and 238±20 km s -1 respectively) are comparable to the typically 
measured values. The results from STEREO-k are qualitatively corroborated 
by those from STEREO-B, with minor differences in the peaks probably due to 
geometric effects. 

The cadence of the 171 A passband was artificially degraded from 2.5 minutes 
to 10 minutes to allow a direct comparison with the 195 A and 304 A pass- 
bands, with the resulting kinematics observed to decrease giving a peak velocity 
~300 km s _1 and a variation in acceleration from ~100 m s~ 2 to ~ —100 m s~ 2 . 
These values are comparable to those estimated for the 195 A and 304 A pass- 
bands, suggesting that the results obtained for these passbands are most likely 
due to temporal undersampling, analogous to the apparent smoothing with de- 
creasing sampling rate exhibited in signal processing. This indicates that the 
effects of temporal sampling must be accounted for when analysing observations 
of CBFs. These effects may explain the lower velocities previously measured in 
SOHO /EIT observations of "EIT waves". The relatively low cadence available 
from EIT (~12 minutes) may have resulted in short time-scale variations of the 
propagation distance being missed. 

Although the pulse is quite faint in the plain images, this event was visible 
in the 304 A passband, marking the first observation of a CBF pulse in this 
passband. This is a result of the consistently high cadence of the 304 A passband 
on STEREO /EUVI and provides a new opportunity to understand the physical 
nature of the disturbance. The 304 A passband is used to observe the upper 
chromosphere using He II emission so a pulse observed using this passband would 
appear to be chromospheric in nature. However, it should be noted that there 
can be a strong contribution from the coronal Si XI 303.32 A emission line in 
this passband. This Si XI line has a peak formation temperature of ~1.6 MK 
(Brosius et al, 1996), which is consistent with the observations made using the 
coronal 171 A, 195 A and 284 A passbands. 

The 304 A data required extensive processing to make the disturbance visible 
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above variable small-scale features present in the running difference images. These 
features are expected due to the dynamic region of the solar atmosphere sampled 
by 304 A, although the wavefront was discernible after filtering. This very small 
signal-to-noise is consistent with the small signal expected from the coronal Si XI 
emission line, and indicates that the observed pulse is coronal in nature. 

The properties of this event should allow clarification as to the physical nature 
of the disturbance, in particular whether it was a fast-mode magnetoacoustic 
wave or the propagation of perturbation sources resulting from CME lift-off. The 
derived velocity of the pulse is higher than previous estimates, and may now be 
consistent with the fast-mode magnetoacoustic wave or propagating MHD shock 
interpretation. As the pulse moves through the quiet corona, it is propagating 
perpendicular to the predominantly radial magnetic field lines of the quiet Sun, 
indicating a fast-mode wave nature. This allows the characteristic wave velocity 
to be calculated using the equation, 

v ms = {cl + v\)^ ) (4.1) 

where the sound speed c s and the Alfven speed va have been defined in Equa- 
tions 1.48 and 1.61 respectively (see Section 1.3 for more details). It has been 
noted by Wills-Davey et al. (2007) that the range of typical Alfven speeds in the 
quiet solar corona is of the order 200-1000 km s" 1 . The pulse velocities found 
in this work are higher than previously estimated values as a result of the higher 
image cadences available, and are now more consistent with the fast magnetoa- 
coustic wave interpretation. 

It should also be noted that this event was associated with a partial halo 
CME and filament lift-off which began at ~12:31 UT. The timescale over which 
the pulse accelerated to its peak velocity is comparable to that for the CME, 
possibly indicating a common launch mechanism. The pulse was also observed 
in the 171 A and 195 A passbands to exhibit some evidence of reflection and 
refraction when interacting with the boundary of an adjacent coronal hole. This 
is consistent with the interpretation of the pulse as a fast-mode magnetoacoustic 
wave pulse. 
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4.6 Further Comments 



The analysis of the 2007 May 19 event outlined above was the first discussion of 
this event to be published (Long et al, 2008) and as a result has been analysed 
in turn by other authors examining this event (such as e.g., Veronig et al, 2008; 
Gopalswamy et al, 2009; Attrill, 2010). Although the major conclusions of Long 
et al. (2008) have been upheld in these re-examinations, some issues have been 
raised with regard to the techniques used. 

The use of a running difference technique in particular, while common and 
historically used for the identification of "EIT waves", has been questioned by 
Attrill (2010). The inherent nature of the running difference technique highlights 
relative motion, and may result in a misidentification of a moving pulse. 

The form of the kinematics derived by Long et al. (2008) are also a source of 
concern. The highly variable nature of the acceleration is inconsistent with all 
previously proposed models, and invites further consideration. The kinematics 
derived by Veronig et al. (2008) for this event, while comparable, do not show the 
same variability. However, both Veronig et al. (2008) and Ma et al. (2009) noted 
the effects of cadence on the temporal undersampling of pulse kinematics (albeit 
for a different event in the case of Ma et al, 2009). The form of the kinematics 
curve obtained for the 171 A passband in particular does not conform with results 
obtained by Veronig et al. (2008), who note that this apparent increasing acceler- 
ation is most likely due to a misinterpretation of expanding coronal loops as the 
initial propagation of the pulse. This highlights the difficulties in discriminating 
between rising loops and the initial propagation of the pulse, particularly when 
the pulse may be initiated by the rising loop. In addition, the initial velocity 
derived for the pulse in the 171 A passband is not consistent with a fast-mode 
magnetoacoustic wave, which should consistently have a velocity that exceeds the 
local Alfven speed. 

Once the results of these analyses had been taken into consideration, it was 
decided to examine in detail the effects of running difference image processing and 
to identify an alternative technique for studying "EIT waves". Simultaneously, 
it was decided to test the capabilities of numerical differencing techniques when 
faced with small data-sets and determine a technique that could be relied upon to 
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produce statistically accurate kinematics estimates regardless of sample size. This 
examination of image processing techniques and numerical differencing schemes 
is discussed in Chapter 5. 
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Chapter 5 



Identification and Analysis 
Techniques 

You see, but you do not observe. The distinction is clear. 

Sir Arthur Conan Doyle (A Scandal in Bohemia) 



Coronal bright fronts are typically observed as broad diffuse features in coronal 
emission lines as previously discussed in Chapter 2. This diffuse nature makes it 
difficult to identify CBF pulses in individual images, complicating analysis of the 
phenomenon. These issues have resulted in the development of several different 
image processing techniques for identifying CBFs, each with advantages and dis- 
advantages. It is shown here that the techniques traditionally used for identifying 
CBFs are fundamentally flawed and prone to undefined user-dependent errors. 
A new semi-automated technique is proposed to allow the rigorous identification 
of CBF pulses. 

Once the pulse has been identified in subsequent images a thorough analysis 
requires a determination of the pulse kinematics. These can be derived using 
several different methods, either by direct numerical differentiation of the data 
or alternatively through the application of a model fit to the data. Again, each 
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technique has drawbacks, particularly when dealing with the typical data-set sizes 
that are encountered when studying CBFs. 

In this chapter, the issues traditionally associated with identifying CBF pulses 
are discussed. The different image processing techniques used are presented in 
Section 5.1.1 with a semi-automated technique for identifying CBFs outlined in 
Section 5.1.2. The different numerical techniques that can be used to derive the 
pulse kinematics directly from distance-time measurements are then presented 
in Section 5.2.1 with an examination of the effects of image cadence and pulse 
uncertainty outlined in Section 5.2.2. It is shown that the errors associated with 
the identification and analysis of CBFs can be minimised through the use of a 
percentage base difference image processing technique allied to a model fit using 
a bootstrapping approach. 

5.1 Pulse Identification 

The diffuse nature of CBFs has made them difficult to observe in single images, 
with identification improved in movies as a consequence of their motion. To 
overcome this issue, CBFs have traditionally been identified using difference im- 
ages, an approach which involves subtracting one image from another in a bid to 
highlight the moving pulse. There are several different types of difference image 
including the running difference (RD), base difference (BD) and percentage base 
difference (PBD) techniques. 

5.1.1 Difference Imaging 

Difference imaging is a simple image processing technique that is designed to 
highlight relative motion in images, making it particularly useful for identifying 
CBF pulses in subsequent images. Some examples of a differencing approach can 
be seen in Figures 2.1 and 5.1, both of which show a CBF pulse highlighted using 
differencing. 

The most common and simplest technique used to accentuate the pulse signa- 
ture is the running difference image (see Section 4.3). This involves subtracting a 
following image from a leading image (i.e., AI — It — h-i), highlighting relative 
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Figure 5.1: Left: Running difference image from 2007 May 19. White (black) 
is a positive (negative) intensity change relative to the previous image. Right: 
Simulated intensity profiles showing a pulse at To (top panel) Ti (middle panel) 
and Ti - To (bottom panel). 

motion with an increase in intensity seen in white and a decrease in intensity 
seen in black (as shown in the left-hand panel of Figure 5.1 and in Figure 5.2). 
While this technique clearly highlights the moving pulse, there are several issues 
that must be kept in mind when examining the resulting images (see also Attrill, 
2010, for some additional issues with the misinterpretation of RD images). 

As shown in the right-hand panel of Figure 5.1, the profile of the pulse seen 
in a RD image does not necessarily conform to that of the original pulse. The 
technique highlights the relative intensity change between images, which allows 
the general rather than specific motion of the pulse to be discerned. The nature 
of the data being analysed means that the resulting RD image will also have to 
be scaled to make the pulse apparent. This scaling is entirely user dependent and 
strongly influences the apparent size of the observed pulse. 

An alternative technique proposed to negate the influence of image-to-image 
variation is the base difference (BD) technique (see Figure 5.2). Here, a pre- 
event image is subtracted from each subsequent image, highlighting the intensity 
changes relative to the pre-event image (i.e., AI — It — Iq). This is designed to 
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Figure 5.2: Plain, running difference (RD), base difference (BD) and percentage 
base difference (PBD) images for the event from 2007 May 19 in the 195 A passband 
from STEREO-A. The times used to produce the images are shown in the top of 
each panel. The plain image is shown using logarithmic scaling, the RD and BD 
images using user-defined intensity scaling and the PBD image using its intrinsic 
intensity range (i.e., -100% to 100%). 

remove only the coronal background, allowing an easier identification of the pulse 
itself. Although this is a better approach for enhancing the true pulse signature 
than running-difference, the small-scale background structure of the corona is also 
highlighted using this technique. Consequently, this technique becomes ineffective 
over long timescales with the resulting image difficult to interpret. 

The BD technique is also subject to user-defined scaling which can over- 
emphasise the pulse signature, while the diffuse nature of the pulse compared 
to the small-scale motion also highlighted by the technique means that image 
smoothing is often required. While this makes the pulse easier to identify, it does 
introduce an additional error that must be accounted for when identifying the 
pulse position. 

A variation on this technique was originally proposed by Wills-Davey & Thomp- 
son (1999) and involves dividing the BD image by the pre-event image (i.e., 
AI = (I t — Iq)/Iq x 100 (see Figure 5.2). This produces images with a relative 
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intensity nominally ranging from —100% to 100%, eliminating the possibility of 
arbitrary image scaling (although it should be noted that maximum intensities 
greater than 100 % can be observed as a result of e.g., flares). The peak pulse 
intensity is also then given as a percentage of the pre-event intensity, allowing 
a better comparison from frame to frame. While this variation eliminates the 
issues with regard to the arbitrary image scaling, the assumption about a quiet 
pre-event coronal background remains. 

One additional approach to minimising the errors associated with the individ- 
ual differencing techniques is to account for the effects of the differential rotation 
of the Sun. This can be done by de-rotating all images to the same pre-event time 
before applying the differencing techniques. While the timescale of a CBF pulse 
is quite short, the de- rotation allows the errors associated with the small-scale 
coronal features to be minimised. This has the effect of reducing the additional 
error at the pulse edges due to small-scale coronal movements and solar rotation, 
allowing a better identification of the pulse. 

5.1.2 Automated Pulse Identification Techniques 

Once identification of the pulse has been confirmed in a difference image, the 
next step in the pulse analysis is to determine the distance of the pulse from a 
defined source point at a given time. The traditional approach for this has been 
to use a point-and-click identification of the edge of the pulse using a RD or BD 
image with a user-defined scaling. However, as shown in Figure 5.3, the scaling 
used for the image defines the position and width of the pulse. 

It is apparent from Figure 5.3 that the pulse is most apparent when a harsh 
(i.e., ±1) scaling is applied to the image. However, while this highlights the full 
extent of the pulse, it also results in an increase in the small-scale variation of the 
background corona. An additional problem with point-and-click pulse identifi- 
cation is defining the location on the pulse to be identified. This is generally the 
front edge of the pulse. However, it is clear from Figure 5.3 that this varies with 
the scaling chosen. Consequently, point-and-click pulse identification is highly 
user-dependent, and not an efficient technique for a detailed analysis of CBFs. 

A more effective approach is to take an intensity profile across the pulse and 
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Figure 5.3: Running difference images with variable photon flux intensity scaling 
(with units of detected photon difference) of (a) ±1, (b) ±10 and (c) ±50. 



identify the location of the peak intensity of the pulse, effectively allowing the 
centre of the pulse to be tracked rather than the edge. Whereas the edge of the 
pulse is loosely defined and highly variable, CBFs have been noted as having 
a roughly Gaussian shape (Wills-Davey & Thompson, 1999), meaning that the 
pulse centre is a more statistically significant indicator of the pulse position. 

A technique was devised to allow semi-autonomous identification of the CBF 
pulse using derotated percentage base difference images. This approach was cho- 
sen as the nominal inherent scaling of the technique used to produce them mit- 
igated against any user-defined assumptions, while the derotation allowed the 
small-scale coronal variation to be minimised. To ensure that the signal-to- 
noise of the pulse was maximised, a great-circle sector (i.e., an area on the sphere 
bounded by two great circles) projected onto the Sun was defined with the source 
point of the pulse defined as the crossing point of the great circles. The intensity 
within this sector could then be averaged across the arc (in annuli of increasing 
radii with 1 degree width on the surface of the sphere) to produce an intensity 
profile as a function of great circle distance along the solar surface away from 
the source location (cf. similar techniques proposed by Warmuth et al, 2004b; 
Podladchikova & Berghmans, 2005; Wills-Davey, 2006; Veronig et al, 2010). 

To minimise the errors associated with identification of the source point, the 
first two observations of the CBF pulse in the 171 and 195 A passbands were 
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Figure 5.4: Left panel: PBD image with arc sector over-plotted. Right panel: 
Intensity profile across arc sector. The relative intensity change at a given dis- 
tance corresponds to the mean intensity enhancement across the arc sector at that 
distance. This is then fitted using a Gaussian model and repeated for the next 
image. The vertical lines in the right-hand plot and arc lines in the left-hand plot 
correspond to the 2a limits of the Gaussian fit. 



visually located and their manually identified outline fitted using an ellipse, with 
the mean centres of the four ellipses taken as the source point of the CBF pulse. 
The orientation of the arc sector was then chosen to maximise the number of 
pulse identifications for a given event, with an arc of 30 degrees chosen as this 
provided the greatest signal while minimising errors resulting from any change in 
the propagation direction of the pulse. 

Figure 5.4 shows how the algorithm works for the event from 2007 May 19 
as observed by STEREO-k. The left-hand panel shows the PBD image with 
the arc sector overplotted, while the right-hand panel shows the 1-dimensional 
intensity profile along the arc sector. The pulse is clearly visible in the arc sector, 
and has been fitted using a Gaussian model allowing the 2a limits of the fit to be 
overplotted as arc-lines in the left-hand panel and vertical lines in the right-hand 
panel. 

This process is repeated for each observation of the pulse, allowing the varia- 
tion with time of the distance of the pulse from the source to be examined. This 
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can then be used to determine the kinematics of the pulse, while the temporal 
variation in the pulse width can be studied for evidence of pulse broadening. The 
next step is to identify the optimal technique for deriving the pulse kinematics for 
a small data-set such as that typically produced by STEREO /EUVI observations 



As well as minimising the errors associated with the actual identification of the 
CBF pulse, it was important to identify the optimum method for deriving the 
kinematics of the pulse. To do this, a number of different numerical differencing 
techniques and a residual resampling bootstrapping approach were tested and 
compared. This was required as a result of the sparse nature of the data ob- 
tained from STEREO /EUVI in addition to concerns with regard to the approach 
previously used to derive the pulse kinematics. 

5.2.1 Numerical Differencing 

The simplest way to determine the pulse kinematics given distance-time measure- 
ments is to use a numerical differencing technique to directly derive the velocity 
and acceleration of the pulse. Here, four different numerical differencing tech- 
niques were identified; a forward, reverse and centre difference method, all of 
which can be derived using a Taylor expansion, and a three-point Lagrangian 
interpolation technique. To allow a direct comparison between the different tech- 
niques, a simulated data-set was created (shown in panel a of Figure 5.5), with 
the numerical differencing techniques used to estimate the known kinematics. 
This simulated data-set was constructed using an equation of the form, 



with r = 50 Mm, v = 400 kms -1 and a = —150 ms~ 2 . The 5r term here 
corresponds to ~5% noise which was added to the simulated data to reproduce 
observational variations. This results in a relative error that increases with prop- 



of a CBF. 



5.2 Statistical Analysis 



r(t) = r + v t + -at 2 + Sr 
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Figure 5.5: Panel (a); Simulated data set (asterisks) with the model kinematics 
shown by the red line and given in the bottom right. Panels (b)-(e); Velocity 
derived using forward-difference (b), reverse-difference (c), centre-difference (d) 
and Lagrangian interpolation (e) techniques. The kinematics determined by the 
fitted dashed line are shown in the bottom left in each case. 
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agation distance and was chosen as it is representative of the increased difficulty 
in identifying an "EIT Wave" pulse at increasing distance from the source. For 
typical observations of an "EIT Wave" , the generally high pulse intensities close 
to the source produce small positional uncertainty that increases with distance 
as the pulse intensity and spatial extent of the pulse decrease and increase re- 
spectively. An observing cadence of 600 s was also used as this is the typical 
observing cadence of STEREO. The resulting simulated data-set is shown as the 
crosses in panel (a) of Figure 5.5, with the known kinematics shown as the red 
line. This simulated data set was then used to compare the different numerical 
differencing techniques. 



5.2.1.1 Forward Differencing 

The first numerical differencing technique to be tested was the forward-difference 
technique. This involves the computation of the derivative at the point r[t + At] 
by extrapolating forward from the point r[t]. The derivative can be determined 
using the Taylor series, 

rltV rltV" 
r[t + At] = r[t] + r[t]'At + -^-(At) 2 + -^j-(At) 3 + ... (5.2) 

By rearranging this equation, it is possible to find the derivative at time t, 

Although this series continues to infinity, the correction values are dominated by 
the term of order At. As a result, this is usually written as, 

r[M] = i±M^iM + 0(Ai) (54) 

where O(At) is the truncation error term. This can also be used to derive the 
acceleration in a similar fashion by replacing r with v to give, 

, W ' = !*±Mz!!W +0(At) . (5.5) 
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However, v[t] is already known from Equation 5.4 above, so by substituting this 
into Equation 5.5, it is possible to derive the acceleration a[t] as, 

alt] = ±tmz^±M±m + (At). (5.6) 

Equations 5.4 and 5.6 can then be used to determine the instantaneous velocity 
and acceleration at the point r[t] for each value of t. 

This approach does have some issues which make it unsuitable for analysing 
CBF observations. The technique calculates the derivative assuming a straight 
line gradient between points, while the final data point is removed with each 
iteration as an inherent consequence of the method. This approach is also par- 
ticularly susceptible to point-to-point variations, with the result that noisy data 
can produce strongly varying velocity and acceleration estimates. 

The results of the forward-difference technique are shown in panel (b) of 
Figure 5.5. The effects of the point-to-point variation are very apparent, while 
the inherent nature of the technique has removed the final data point. The 
kinematics obtained by fitting the data points also show a strong underestimation 
of vq, while there are large errors associated with the estimates of both v and a. 



5.2.1.2 Reverse Differencing 

The reverse differencing approach uses a similar methodology, although in this 
case the derivative at the point r[t] is calculated by extrapolating backwards to 
the point r[t — At]. This can be determined using the Taylor series expansion, 

r\t\" rltV" 
r[t - At] = r[t] - r[t]'At + -^(At) 2 - ^(At) 3 + ... (5.7) 

Once again, this equation can be rearranged to find the derivative at time t, 

ri . r[t]-r[t-At] r[t]", A . r[t]"' . 
r[t]' = U A_ + -y-(Af) - ^r(At) 2 + ... (5.8) 
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The correction values here are again dominated by the term of order At, allowing 
this equation to be rewritten as, 



r[t]' = v[t] = r[t] f t At] +0(At) (5.9) 

where O(At) is the truncation error term. As before, the acceleration term can 
also be derived by replacing r with v to give, 

v[tY= m^_At] +om (510) 

Again, v[t] is already known from Equation 5.9 above, so by substituting this into 
Equation 5.10, the acceleration a[t] can be derived in terms of r[t] as, 

r[ t ]-2r[ t -A f ] + r[ f -2A t ] + 
At 2 

Both equation 5.9 and 5.11 can then be used to determine the instantaneous 
velocity and acceleration at the point r[t] for each value of t. 

Panel (b) in Figure 5.5 shows the result of this analysis for the simulated data- 
set. The reverse-difference approach is very similar to the forward-difference 
technique, and this is reflected in the derived data-points which have identical 
velocity values albeit at a different time. Once again, the point-to-point variation 
is significant, while the first data-point has also been removed. The fitted kine- 
matics show a higher estimated initial velocity but identical acceleration, with 
the initial velocity not corresponding to that of the original model. 

The same caveats that apply to the forward difference technique also apply 
to the reverse difference technique due to their inherent similarities in using two 
adjacent points (r[t\ k, r[t + At] and r[t] & r[t — At] for the forward and re- 
verse difference respectively) to determine the derivative at a given point. Both 
techniques are quite simplistic and as a result are only applicable in simplistic 
cases. Both techniques also remove points from either edge with each iteration 
as an inherent consequence of the technique. As a result, neither technique is 
appropriate for determining the kinematics of a CBF pulse. 
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5.2.1.3 Centre Differencing 



A more accurate approach that utilises the Taylor series methodology is the 
centre difference technique. This uses the points either side of the point under 
examination r[t] (i.e., r[t — At] and r[t + At}), smoothing the data and producing 
a more accurate estimate of the numerical derivative at that point. The Taylor 
series at the points r[t — At] and r[t + At] respectively are defined as, 

rffl" r \tV» 
r[t - At] = r[t] - r[t]'At + -^(At) 2 - ^(At) 3 + ... (5.12) 

rltV rltV" 
r[t + At] = r[t] + r[t]'At + ^-(At) 2 + "^-(A*) 3 + - (5.13) 

as shown previously in Equations 5.2 and 5.7 above. The centre difference ap- 
proach subtracts Equation 5.12 from Equation 5.13 to produce, 

rltV" 

r[t + At] - r[t - At] = 2r[t\ At + 2^-(At) 3 + (5.14) 

which can then be rearranged in terms of r[t]' to give, 

r[t] , = r 1±M _^- M + rJ^ m2+ (515) 

Following the notation from before, this can be rewritten as, 

r[tY = v[t] = r 1± Att^t^ + 0{At2) (516) 

where 0(At 2 ) is the truncation error term. The instantaneous acceleration term 
a[t] can then be determined as, 

r[tr = a[t] = r lt + A t] -Mtl 2+ r lt -^ ] +0{At2) (517) 

This allows the instantaneous velocity and acceleration at the point r[t] to be 
determined. 

The results of applying this technique to the simulated data-set are shown 
in Panel (d) of Figure 5.5. The centre-difference technique produces a smoother 
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estimation of the kinematics than the forward- and reverse-difference and the 
fitted kinematics are closer to the model values. However, the centre difference 
technique does remove both of the edge points as an inherent consequence of its 
operation. 

Despite producing a better estimate of the pulse kinematics, the inherent is- 
sues of the centre-difference technique are a cause for concern. The removal of 
edge points makes analysis of the kinematics difficult, particularly when consid- 
ering that a given STEREO /EUVI observation of a CBF event may have only 
3 or 4 distance-time estimates. Consequently, the centre difference technique is 
usually unsuitable for this analysis. 



5.2.1.4 Lagrangian Differencing 

A more advanced technique for determining the derivative at a point numeri- 
cally is the three-point Lagrangian interpolation technique used by the built-in 
deriv.pro routine in IDL (see Bevington & Robinson, 2003, for details). This 
method uses three adjacent points (r[t — At], r[t] and r[t + At]) to fit a La- 
grangian polynomial function of the form, 

(x - t)(x - (t + At)) 

F[X) ~ {{t-At)-t){{t-At)-{t + At)) r[t ~ Atl 

(x-(t-At))(x-(t + At)) 

(t - (t - At))(t - (t + At)) 11 

(x - (t - At))(x - t) , . . 

+ 77 A 7 t — 7777 -r[t + At] 5.18 

((* + At) - (t- At))((t + At) -t)) 11 K J 

where x is the abscissa of the point r[t]. The derivative of this is then given by, 

, _ 2x-t-(t + At) 

{X} ~ ((t-At)-t)((t-At)-(t + At)) r[t ~ 1 
2x-(t-At)-(t + At)) 
+ (t-(t-At))(t-(t + At)) 11 
2x-(t-At)-t) 

+ 77 7 a 7777 -J77 77^ [t + At 5.19 

((t + At - (t-At ))({t + At )-t)) 11 K J 
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By substituting t for x, the polynomial used to determine the derivative P'(t) is 
given as, 



This approach does not assume a straight line gradient between points (unlike 
the forward, reverse and centre difference techniques above), while also compen- 
sating for edge points by increasing the errors and retaining all of the original 
data points. The Lagrangian technique also smoothes the data, removing the 
spiky appearance produced by the Taylor series expansion techniques. This is 
apparent in panel (e) of Figure 5.5, where the Lagrangian technique has retained 
all data-points. 

Despite this, there are a number of issues associated with this approach, in 
particular as a result of the small data sets associated with these disturbances. 
The cadence of STEREO /EUVI means that there are typically only 4 or 5 ob- 
servations of a CBF pulse per event, of which three are required to determine 
the numerical derivative of one data point. The combination of this with the 
interpolation at the edges has the effect of skewing the edge points, giving the 
impression of kinematic variation that may not be real. This can be seen in the 
fit to the derived Lagrangian kinematics in panel (e) of Figure 5.5, where the 
edges points have strongly influenced an underestimation of the kinematics. 

This skewing effect strongly suggests that the Lagrangian interpolation tech- 
nique is inappropriate for the small data-sets typically obtained for CBF analysis. 
The problems associated with each numerical differencing approach studied here 
indicates that an alternative method must be used to determine the kinematics 
of a CBF pulse. A possible alternative is the fitting of a specified model to the 
distance-time measurements, although once again the small data-sets available 
can make this difficult. 

5.2.2 Effects of Image Cadence and Positional Uncertainty 

In the absence of a numerical technique to determine the kinematics of a CBF 
pulse, the cadence of the observing instrument and the degree of uncertainty 
in the data must also be considered. Here, the influence of both cadence and 
uncertainty are examined using the simulated data discussed in Section 5.2.1. 




(5.20) 
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Figure 5.6: Simulated data (crosses; derived using the model shown in the top 
panel of Figure 5.5) for cadences of 720 s (top; comparable to EIT), 300 s (middle; 
comparable to EUVI) and 12 s (bottom; comparable to AIA). In each case, the fit 
to the data is shown in red with derived kinematics in the bottom right of each 
panel. 



The effects of varying image cadence were the first to be examined; this is shown 
in Figures 5.6 and 5.7. 

The same data-set was used for each of the panels shown in Figure 5.6, with 
the cadence varied to best reflect data from S OHO /EIT (top panel at 720 s), 
STEREO /EUVI (middle panel at 300 s) and SDO/AIA (bottom panel at 12 s). 
The derived kinematics are given in the bottom right of each panel with the fit 
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Figure 5.7: Derived kinematics for varying image cadence with ±5% uncertainty. 
Distance, velocity and acceleration are shown in the top, middle and bottom panels 
respectively, with the x-axis shown using logarithmic scaling to highlight the effects 
of varying the image cadence. 



to the data shown by the red line. It is immediately apparent that the higher 
cadence data allows a better fit to the data despite the uncertainty in the data 
(which was kept at ±5% of the model for each data-set). This implies that higher 
cadence data is required to derive the true kinematics of a CBF pulse. 

The variation in derived kinematics with cadence is shown in Figure 5.7 (again 
for ±5% uncertainty). As the cadence decreases, the derived velocity and acceler- 
ation approach the model values, with the scatter showing a dramatic reduction 
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Figure 5.8: Simulated data (crosses) for noise distributions with widths of ±20% 
(top), ±10% (middle) and ±5% (bottom) of the data value. In each case, the fit to 
the data is shown in red with derived kinematics in the bottom right of each panel. 



below ~50 s cadence. These results are consistent with the observations made by 
both Long et al. (2008) and Ma et al. (2009) and show that the effects of image 
cadence must be accounted for when trying to derive the true kinematics of a 
CBF pulse. 

The effects of varying uncertainty in the data were also examined for constant 
image cadence with the results of this analysis shown in Figures 5.8 and 5.9. 
Figure 5.8 shows the derived kinematics for the simulated data-set with the top 
panel showing uncertainties of ±20%, the middle panel showing uncertainties of 
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Figure 5.9: Derived kinematics for varying uncertainty at 150 s image cadence. 
Distance, velocity and acceleration are shown in the top, middle and bottom panels 
respectively. 



±10% and the bottom panel showing uncertainties of ±5%. The variation in 
uncertainty has a noticeable effect on the derived kinematics in each case, with 
the acceleration in particular exhibiting significant variation with uncertainty. 

The variation in the derived kinematics with uncertainty is shown in Fig- 
ure 5.9, with the variation in offset distance (r ), initial velocity (v ) and acceler- 
ation (a ) shown in the top, middle and bottom panels respectively. Some vari- 
ation is apparent in each case, with the acceleration again showing the strongest 
variation. The derived kinematics are seen to approach the model kinematics as 
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the noise reduces to 0% as expected. 

It is clear that the variation in both noise and imaging cadence can strongly 
influence the derived kinematics of a CBF pulse, while the different numerical 
differencing techniques do not return accurate estimates of the pulse kinematics. 
However, it is possible to use a bootstrapping approach to overcome these issues 
as this is a statistically significant technique that has been optimised for small 
data-sets such as those typically obtained when studying CBF pulses. 

5.2.3 Bootstrapping 

Bootstrapping was first proposed by Efron (1979) as a general interpretation of 
the jackknife, allowing statistical quantities to be estimated from a limited sample 
to a high degree of accuracy. In principle this allows a statistical quantity (such 
as the mean, variance etc.) of a distribution to be estimated using a small sample 
of measurements from that distribution. The resulting estimates are statistically 
significant, with statistically defined associated error estimates. 

While there are several different bootstrapping techniques, each of which is 
useful for a different purpose, in this case the residual resampling technique will 
be used. Whereas the majority of bootstrapping techniques involve removing one 
point at random and recalculating the fit to the data, this is inappropriate in this 
case due to the small sample set available. The residual resampling technique 
however, does not require data-point removal, making it more appropriate for 
small data-sets such as those available here. 

The residual resampling technique works as follows. The distance-time mea- 
surements (i/i] i = 1,2, ... ,n) are first fitted using a quadratic equation of the 
form 

r(t) = r + v t + X -at 2 , (5.21) 

where r is the offset distance, v is the initial velocity, a is the constant accel- 
eration of the pulse, and t is the time elapsed since the first observation of the 
disturbance. This yields the fitted values (yi) and residuals (e = yi — jji) for each 
point. The residuals are then randomly ordered and randomly assigned a sign 
(1 or —1), before being added to the original fit values to produce a new data 
set. These new data are then fitted using the same model, with the resulting 
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Figure 5.10: Top: Simulated data corresponding to ~5% noise added to the 
model (asterisks) with bootstrapped fit shown in red and parameters in bottom 
right. Bottom: Histograms showing the distance (a), velocity (b) and acceleration 
(c) derived using the bootstrapping technique. 



re-fitted parameters recorded. The process is then repeated a large number of 
times (~10 000). 

This has the effect of repeating the fit to slightly varying data a large number of 
times, allowing a distribution of fit values to be obtained. The mean and standard 
deviation of this distribution then provide the estimated value and associated 
error of the parameter. This technique is statistically rigorous and produces a 
more accurate result than a simple model fit to the given data. 
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This technique was applied to the simulated data-set used in Section 5.2.1 
(corresponding to the model with ~5% noise added) to allow a comparison of 
the effectiveness of the bootstrapping approach with the numerical differencing 
techniques. The results of this analysis are shown in Figure 5.10. The upper panel 
here shows the simulated data (indicated by the crosses) with the bootstrapped fit 
to the data shown by the red line. The kinematics obtained by the bootstrapping 
approach are given in the bottom right. The bottom panels (a), (b) and (c) show 
the histograms for the distance, velocity and acceleration parameters respectively. 
These histograms returned by the bootstrapping approach for each parameter 
allow a statistical determination of the fit parameters, a fact reflected in the error 
term for the derived kinematics. 

The kinematics returned by the bootstrapping approach match the model 
kinematics within one standard deviation, while also producing an accurate es- 
timate with well-defined errors, unlike those of the numerical differencing tech- 
niques. It should also be noted that none of the data-points have been removed, 
unlike the majority of the numerical differencing techniques tested. 

As a result of the simulations carried out to test the different approaches to 
determining the kinematics of a CBF pulse, the residual resampling bootstrap- 
ping technique was chosen as the best method. This fits a model directly to the 
distance-time measurements, allowing the kinematics of the CBF pulse to be de- 
termined to a high degree of accuracy. This is a statistically significant approach 
that removes much of the error inherent in the Lagrangian interpolation technique 
used in the previous analysis discussed in Chapter 4. The kinematics for each 
future event studied here therefore reflects the mean value of the bootstrapped 
distribution, with the associated error given by the standard deviation. 

Note that the bootstrapping technique outlined here requires the application 
of a pre-defined model to the data, and as such does not solve the inherent prob- 
lems associated with numerical differencing. However, bootstrapping does provide 
a statistically significant approach to determining the errors associated with the 
pulse kinematics without introducing any additional errors into the data. This 
makes it a very useful technique for determining the kinematics of CBF pulses, 
particularly considering the small data-sets available from STEREO /EUVI. 
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5.3 Conclusions 



As outlined in this chapter, there are many different issues beyond the simple 
identification and analysis of the pulse that must be examined and accounted 
for when studying CBFs. The choice of technique used to identify the pulse 
can have a strong effect on the estimated pulse parameters, while the resulting 
kinematics are influenced by the technique used to estimate them. The cadence 
of the observations and the uncertainty in the pulse position (which is influenced 
by the resolution of the instrument) can also affect the derived kinematics. 

Of the different image processing techniques considered here, the percentage 
base difference images offer the optimum approach for identifying the CBF pulse. 
This technique removes the background corona, highlighting the actual pulse 
rather than an intensity change relative to the previous observation. The inherent 
scaling afforded by the percentage intensity approach also removes any ambiguity 
with regard to the identification of the pulse. 

The combination of PBD images with the semi-automated algorithm pre- 
sented here allows the CBF pulse to be identified in an unbiased way, negating 
the user-defined errors associated with traditional methods of pulse identification. 
The algorithm also uses a consistent approach to identify the pulse, with the use 
of a Gaussian model ensuring that the pulse position and width are comparable 
from both image to image and event to event. This will allow a direct comparison 
of multiple events and represents a opportunity to systematically examine CBFs. 

Once the CBF pulse has been identified and tracked, the next step involves 
determining the pulse kinematics. This has previously been done using numerical 
differencing techniques, but it has been shown here that this approach is incom- 
patible with the typical cadence of CBF observations. The different numerical 
differencing techniques based on a Taylor expansion approach remove data points 
as a direct consequence of their operation, while the Lagrangian interpolation ap- 
proach is strongly affected by the edge points. 

The derived kinematics have also been shown to be strongly affected by the 
cadence of the observations and the associated uncertainty of the pulse position. 
By comparing derived kinematics for a simulated data-set using cadences com- 
parable to those available from SOHO/EYT, STEREO /EUVI and SDO/MA it 
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was shown that the derived kinematics approached those of the model as the 
cadence decreased. This is analogous to the aliasing effect observed in signal 
processing and suggests that historical cadences were insufficient to derive the 
true kinematics of CBF pulses. The uncertainty in pulse position also influences 
the derived kinematics, with a decrease in the positional uncertainty producing 
derived kinematics that approach those of the model. This indicates that it will 
be possible to derive the true kinematics of a clearly identified CBF pulse using 
high cadence data from SDO/AIA. 

The issues surrounding the pulse identification and derivation of the resulting 
kinematics can be negated by combining the semi-automated algorithm presented 
here with a residual-resampling bootstrapping technique. This allows any errors 
associated with the identification of the CBF pulse to be minimised while also 
producing statistically significant estimates of the pulse kinematics with defined 
associated errors. This approach is used to examine a series of CBF events 
observed using the STEREO spacecraft in the next chapter. 
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Chapter 6 



Deceleration and Dispersion of 
Large-scale Coronal Bright 
Fronts 

A man can't prove anything without statistics 

Mark Twain 

One of the most dramatic manifestations of solar activity are large-scale coro- 
nal bright fronts (CBFs) observed in extreme ultraviolet (EUV) images of the 
solar atmosphere. To date, the energetics and kinematics of CBFs remain poorly 
understood, due to the low image cadence and sensitivity of previous EUV im- 
agers and the limited methods used to extract the features. In this chapter, 
work performed to determine the temporal variation of the properties of a sam- 
ple of CBFs is outlined, with particular emphasis placed on the variation in their 
kinematics and pulse shape, dispersion and dissipation. This analysis followed 
the development of a semi-automatic intensity profiling technique to extract the 
morphology and accurate positions of CBFs in 2.5-10 min cadence images from 
STEREO /EUVI. By applying this technique to sequences of 171 A and 195 A im- 
ages from STEREO /EUVI, it was possible to measure the wave properties of four 
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separate CBF events from 2007 May 19 (see Chapter 4 for the original analysis 
of this event), 2007 December 07, 2009 February 12 and 2009 February 13. 

Following their launch at initial velocities of ~240-450 km s _1 each of the four 
events studied showed significant negative acceleration ranging from ~ —290 to 
—60 m s~ 2 . The spatial and temporal widths of the CBF pulses were found 
to increase from ~50 Mm to ~200 Mm and ~100 s to ~1500 s respectively, 
indicating that they are dispersive in nature. The variation in position-angle 
averaged pulse-integrated intensity with propagation shows no clear trend across 
the four events studied. These results are most consistent with CBFs being 
dispersive magnetoacoustic waves. This chapter contains results published in 
Astronomy & Astrophysics by Long, Gallagher, McAteer & Bloomfield (2011a). 

6.1 Introduction 

As the STEREO mission continued, the number of observed CBF events increased 
with 10 wave events observed during the period 2006 Dec 12 to 2008 Nov 4 
(Aschwanden et ai, 2009). The increasing number of events combined with the 
questionable abilities of the typically used difference techniques to identify the 
true position of the propagating pulse (cf. Attrill, 2010) has necessitated the 
development of a more automated approach. This was particularly important 
given the approaching launch of the SDO spacecraft, with its ~12 s cadence 
in multiple passbands implying ~40-50 observations per event rather than the 
~5-10 observations per event available from STEREO /EUVI. 

While the launch of SDO has focused attempts to automate detection of 
CBFs, this is not a new issue, with Podladchikova & Berghmans (2005) originally 
proposing a technique to identify and track "EIT waves" based on their associ- 
ated dimming regions. This technique forms the basis of the online Novel EIT 
wave Machine Observing (NEMO) catalogue 1 which provided real time charac- 
terisation of coronal dimming regions until 2010 August 1. Here, the "EIT wave" 
is identified using the coronal dimming region formed by evacuation of plasma as 
the associated CME erupts. 

1 http : //side . oma. be/nemo/ 



106 



This work was followed by Wills-Davey (2006), who proposed a technique to 
track a CBF pulse in data from the TRACE spacecraft. While this technique does 
provide a way to track a pulse using Huygens tracking, the technique does not 
appear to have been fully automated at this stage. This approach was also never 
applied to data from SOHO/EYT, with only one event described from TRACE. 

Despite these initial attempts to automatically track CBF pulses, not much 
work has been done with regard to applying these techniques to the higher cadence 
data available from STEREO /EUVI. The previously proposed techniques also 
suffer from some shortcomings, with the algorithm proposed by Podladchikova 
& Berghmans (2005) relying on detecting the dimming region rather than the 
actual pulse. While this is preferable for space weather forecasting and detection, 
it is also prone to false-positive detections and does not allow the properties of 
the actual pulse to be examined. 

To eliminate the issues with a point-and-click approach to pulse identifica- 
tion and to overcome the errors associated with running difference images, a 
semi-automated algorithm was devised to analyse CBF pulses. This algorithm 
(described in more detail in Section 5.1) uses de-rotated percentage base dif- 
ference images for pulse identification, neutralising the errors associated with 
running difference images. It should be noted that this algorithm, while effective 
at identifying the pulse, is not yet fully automated and does require user input 
in identifying optimal events for analysis. Once these events have been identified 
however, the actual detection and tracking of the pulse is fully automated. 

Once the issues regarding the identification of the CBF pulse and the deriva- 
tion of the kinematics from distance-time measurements had been resolved (as 
previously discussed in Chapter 5), a sample of four CBF events was chosen for 
analysis. These events were well observed by the STEREO spacecraft, with two 
of them (from 2007 May 19 and 2007 December 07) observed on-disk by both 
STEREO spacecraft. However, the ~90 degree separation of the STEREO space- 
craft meant that the 2009 February events were observed in profile by STEREO-A 
and on-disk by STEREO-B. As a result, only the data from STEREO-B were 
used here for both of these events. 

The algorithm described here is designed to return the estimated position of 
the centroid of the pulse, the pulse width and the integrated intensity of the pulse 
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for each image studied. This allows the kinematics of the pulse to be determined 
using the bootstrapping approach outlined in Section 5.2.3, while the temporal 
variation in pulse width can be examined for evidence of pulse broadening. The 
variation in the integrated intensity of the pulse can be used to determine if the 
pulse is dissipative or non-dissipative. This is important for identifying the true 
physical nature of the pulse. 

The observations used for this work are outlined in Section 6.2, with the algo- 
rithm used to identify and track the pulse summarised in Section 6.3 (although 
a more detailed description can be found in Section 5.1.2). The results of this 
analysis are then presented in Section 6.4, with the 2007 May 19 event discussed 
in detail, allowing a comparison to be made with the analysis previously pre- 
sented in Chapter 4. The application of the algorithm to the additional events 
is outlined in Section 6.5, before some conclusions are drawn about the resulting 
implications for CBFs in Section 6.6. 

6.2 Observations 

The different data discussed here were obtained using the EUVI telescope on- 
board the STEREO spacecraft. Although EUVI observes the Sun in four EUV 
passbands (304 A, 171 A, 195 A, and 284 A), the nature of the passbands meant 
that only the 171 A and 195 A passbands were used here. This is primarily due 
to the high temporal cadence of both passbands (1.5-2.5 minutes for 171 A and 
5-10 minutes for 195 A) and also as CBFs are more readily observed (i.e., of 
higher contrast) in these two passbands. While CBFs have been observed in the 
284 A passband (Zhukov & Auchere, 2004) and the 304 A passband (Long et al, 
2008), the nature of the data and the lower cadence make it difficult to use these 
passbands for more rigorous analysis. 

Figure 6.1 shows percentage base difference images of the four events used 
here, each taken using the 195 A passband onboard STEREO-B. These events 
were chosen as the erupting active region is close to disk centre for each event, 
reducing the errors associated with measuring distances along the curvature of 
the Sun. The events from 2007 May 19 and 2007 December 07 occurred during 
the early phase of the STEREO mission, with the result that the events could 
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Figure 6.1: Percentage base difference images of the 2007 May 19 (top left), 
2007 Dec 07 (top right), 2009 Feb 12 (bottom left) and 2009 Feb 13 (bottom 
right) CBFs, each in the 195 A passband as seen by STEREO-B. The solid lines 
indicate the region in which the disturbance was identified using the intensity profile 
technique, while the cross marks the estimated origin of the event. White (black) 
is an increase (decrease) in intensity from the base image. 



be observed close to disk centre by both STEREO spacecraft. However as pre- 
viously noted, the events from 2009 February occurred later in the mission when 
the STEREO spacecraft were near quadrature. For these events, the erupting 
active region is at disk centre for the STEREO-B spacecraft, but on the limb as 
observed by STEREO-A, with the result that only data from STEREO-B could 
be considered here for these events. 

Each of the events studied was associated with a CME eruption and a GOES 
X-ray flare, with each event also studied independently by other authors. The 
2007 May 19 event (as noted in Chapter 4) has been studied in detail by authors 
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including Veronig et al. (2008), Gopalswamy et al. (2009) and Attrill (2010), with 
the 2007 December 07 event analysed by Ma et al. (2009). The unique nature of 
the February 2009 quadrature events meant that they were also popular choices 
for further analysis, with Kienreich et al. (2009) and Patsourakos & Vourlidas 
(2009) in particular using the events to study the three-dimensional nature of 
the pulses. 

The clarity of the chosen events made them ideal candidates for examina- 
tion using the detection algorithm, while the similar imaging cadence for each 
event minimised the effects of undersampling noted in Chapter 4. The observ- 
ing cadence of the 195 A passband was 600 s for each event, while the 171 A 
passband operated at a cadence of 150 s for each event with the exception of the 
2009 February 13 event where it had a cadence of 300 s. 

In each case, the CBFs were studied using de-rotated base-difference (BD) 
and percentage base-difference (PBD) images. The PBD images were used to 
identify the pulse and determine the variation in pulse position and width with 
time, while the BD images were used to determine the variation in integrated 
intensity with time. The peak intensity of the BD rather than the PBD intensity 
profile was used as this shows the actual intensity increase of the pulse above 
a constant background, while the PBD intensity shows the ratio of the pulse 
intensity increase with respect to the background, and consequently does not 
necessarily produce consistent results for different background values. Both of 
these techniques also indicate the entirety of the CBF, rather than the portion of 
the pulse which has moved beyond its previous extent as is the case with running 
difference images. 

6.3 CBF Pulse Detection 

The algorithm previously described in Section 5.1.2 was used to identify, track 
and analyse the observed CBF pulse for each event studied. In each case, the first 
two observations of the pulse in both the 171 A and 195 A passbands were fitted 
using ellipses, with the means of the four ellipse centres chosen as the projected 
source of the CBF disturbance. This point served as the crossing point of two 
great circles projected onto the Sun, with the sector between the great circles 
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defined as the arc into which the CBF pulse propagated (see Figure 5.4). The 
intensity of the PBD image was then averaged across the position angle (PA) 
within this arc sector for annuli of increasing radii with 1 degree width on the 
surface of the sphere. This was used to produce a 1 -dimensional intensity profile 
as a function of distance away from the source location. This process was repeated 
using BD images to obtain the BD intensity profile. 

The intensity profile produced by this technique was then fitted using a Gaus- 
sian profile (as CBFs were noted by Wills-Davey, 2006, to have a Gaussian form). 
This allows the properties of the pulse to be identified, with the Gaussian fit pro- 
viding the position of the pulse centroid, the width of the pulse and also the peak 
pulse intensity at a given time. The fit also identifies an associated error for each 
of these parameters, allowing the quality of each parameter to be examined. 

The resulting intensity profiles from both passbands studied for the event of 
2007 May 19 are shown in Figures 6.2 (for STEREO-k) and 6.3 (for STEREO- 
'S). In each case, the crosses show the PBD intensity values, with the solid line 
showing the Gaussian fit to the data. Only the first three observations in the 
195 A passband and the first six observations in the 171 A passband from both 
spacecraft were used for this event. Beyond this, the x 2 values of the Gaussian 
fits were too large for accurate analysis. 

The pulse is clearly identifiable in the intensity profiles shown in both Fig- 
ures 6.2 and 6.3, with the model providing a good fit to the data. As previously 
noted in Section 5.1.2, the intensity profile technique is more accurate than a 
normal point-and-click technique as it is semi-automated and reproducible, with 
any associated errors quantifiable as they result from the fitting of a Gaussian to 
an intensity profile across the pulse. In contrast, the point-and-click methodology 
commonly used to identify CBFs is highly user-dependent, with large unquantifi- 
able errors in the identification of the pulse position. Although a pulse may be 
tracked over larger distances using a visual method, increased image processing 
is often required to improve the visibility of the pulse, with the variable scaling 
making the comparison of measurements between images difficult. However, the 
intensity profile technique used here applies the same processing for all images, 
with the result that measurements of the pulse are directly comparable. The in- 
tensity profile technique can also be used to process large amounts of data (such 
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Figure 6.2: PBD image intensity profiles (crosses) for the 2007 May 19 event as 
obtained from STEREO-A in the 195 A (left) and 171 A (right) passbands. The 
Gaussian fit (solid curve) to the positive section of each profile has been overplotted 
on each panel. The time of the leading image (It) in each case is on the upper right 
of the panel. 
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Figure 6.3: Same as Figure 6.2 but for the 2007 May 19 event observed by 
STEREO-B. 
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as those from SDO) rapidly as it does not require the same degree of user input 
as the visual method. 

6.4 Data Analysis 

The CBF identification algorithm discussed in Section 6.3 was used to study data 
from the 2007 May 19, 2007 December 07, 2009 February 12 and 2009 February 13 
events, allowing the temporal variation of the pulse properties to be examined. 
The residual resampling bootstrapping technique outlined in Section 5.2.3 was 
then used to provide statistically accurate estimates of the pulse kinematics for 
each event, while the dispersion and dissipation of the pulse were studied using 
the temporal variation in pulse width and integrated intensity respectively. 

The results of this analysis for the event from 2007 May 19 are discussed here 
in detail, with particular emphasis placed on the methods used to identify each 
parameter. The kinematics, dispersion and dissipation of the pulse are outlined 
in Sections 6.4.1, 6.4.2 and 6.4.3 respectively. This event is described in detail 
as it was the event originally discussed by Long et al. (2008) and outlined in 
Chapter 4. The additional events studied using this technique are presented in 
Section 6.5, with the event from 2007 May 19 used to provide additional context. 

6.4.1 Pulse Kinematics 

The Gaussian model applied to the positive section of the intensity profiles out- 
lined in Section 6.3 was used to plot the temporal variation of the pulse centroid 
distance from the source point. The resulting plot is shown in Figure 6.4 for 
STEREO- A (top) and STEREO-B (bottom), with the errors associated with 
each data-point given by the error associated with the mean of the Gaussian fit 
to the intensity profile in each case. Data from the 171 A and 195 A passbands 
were combined as they have been observed to follow similar kinematical curves 
(Warmuth et al, 2004b; Veronig et al, 2010). The relatively low observing ca- 
dence of the 195 A passband in particular also meant that the combination of 
data from both passbands was necessary to derive accurate pulse kinematics. 
The distance-time data shown in Figure 6.4 were then fitted using a quadratic 
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Figure 6.4: Distance-time plots for STEREO-A (top) and STEREO-B (bottom). 
The 171 A (circles) and 195 A (t riangles) data have been combined as they follow 
similar kinematical curves. The mean offset distance, initial velocity, and accelera- 
tion terms resulting from the bootstrapping analysis (fit indicated by dashed line) 
are also stated in the bottom right of each panel, with errors represented by the 
standard deviation. The errors on each point are given by the error on the mean 
of the Gaussian fit applied to the intensity profile. 



model of the form, 

r(t) = r + v t + l -at 2 , (6.1) 

where r is the offset distance, vo is the initial velocity, a is the constant accel- 
eration of the pulse, and t is the time elapsed since the first observation of the 
disturbance. The data were analysed using the residual resampling bootstrapping 
technique shown in Section 5.2.3 to provide a statistically accurate estimate of 
the fit parameters given the small data-set. 

The values given in the inset of Figure 6.4 show the mean and standard 
deviation of the bootstrapped parameter distributions (cf. the bottom row of 
Figure 5.10). The offset distance (r ), initial velocity (v ) and acceleration (a) 
are given using native SI units of Mm, km s _1 and m s~ 2 respectively for clarity. 
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The estimated kinematics clearly show that the pulse has an initial velocity of 
~444 km s _1 and a deceleration of ~ —255 m s~ 2 as observed by both spacecraft. 

Although the kinematics estimated using both spacecraft are similar, the as- 
sociated errors are slightly different with velocity errors of ±76 & ±87 km s _1 
and acceleration errors of ±136 & ±150 m s" 2 for STEREO- A and STEREO-B 
respectively. While the errors associated with the deceleration are quite large, 
in each case the acceleration is negative within the 1-a error range. The kine- 
matics are also consistent between both spacecraft, indicating that these are the 
true kinematics of the pulse, with the slightly different errors most likely due to 
geometrical effects. 

The initial velocity values given here are similar to previous estimates obtained 
by Long et al. (2008) (see Chapter 4) and Veronig et al. (2008) who both studied 
this event, although in the case of Long et al. (2008) a three-point Lagrangian 
interpolation technique was used to examine the CBF kinematics. The limitations 
of this technique have been discussed in detail in Section 5.2.1.4, with the artificial 
trends introduced to the data sometimes leading to a misinterpretation of the 
derived velocity and acceleration plots. In contrast, the kinematics derived here 
are statistically significant and supported by quantifiable associated errors. 

The reexamination of the kinematics of the 2007 May 19 event presented 
here show a sharp contrast with those presented by Long et al. (2008) and in 
Chapter 4. A comparison with the top row of Figure 4.3 shows strong similarities 
with the positions of the pulse in both the 171 and 195 A passbands, while 
features which were originally identified as the pulse but have since been classified 
as oscillating coronal loops have been removed through the use of the semi- 
automated identification technique. The pulse kinematics derived here using the 
bootstrapping technique exhibit none of the variation observed in Figure 4.3 as 
a result of the model applied to the distance-time measurements. While there 
may be some point-to-point variation in pulse kinematics, it is not possible to 
test for this using the small data-sets currently available from STEREO /EUVI. 
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Figure 6.5: Variation in pulse spatial width (Ar) with distance for the 
2007 May 19 event observed by STEREO- A (top) and STEREO-B (bottom). 
Pulse spatial width here refers to the FWHM of the fitted Gaussian pulse (i.e., 
Ar = 2<7\/21n2). The dashed line indicates the best linear fit to the data. 



6.4.2 Pulse Width 

The parameters returned by the Gaussian model fitted to the intensity profile 
also allowed the variation in pulse width to be examined in addition to the pulse 
kinematics. To do this, the full width at half maximum (FWHM; Ar = 2av / 2m2) 
of the Gaussian fit was determined for each intensity profile studied. The variation 
in pulse width with distance is shown in Figure 6.5 for STEREO-k (top) and 
STEREO-B (bottom), with measurements from the differing passbands indicated 
by the different symbols. It is clear from these plots that an increase in pulse width 
with distance is present for both passbands as observed by both spacecraft. This 
is indicative of pulse broadening and shows that the CBF spreads out spatially 
as it propagates. 
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Figure 6.6: Variation in pulse temporal width (At) with time as observed by 
STEREO- A (top) and STEREO-B (bottom). The dashed line in all panels indi- 
cates the best linear fit to the data. 



The variation in the temporal width of the pulse with time was also considered 
in addition to the variation in the spatial width of the pulse with distance. The 
temporal width of the pulse, At, is denned as, 

Ar 

At = — , (6.2) 

"^pulse 

where Ar is the spatial width of the pulse and v pu \ se is the pulse velocity. While an 
increase in the spatial width of the pulse may indicate pulse dispersion, this could 
be negated by increasing pulse velocity, producing a pulse with constant temporal 
width. By examining the variation of the temporal width of the pulse with time 
it is possible to determine if the pulse is indeed dispersive. The resulting plots 
showing the variation in pulse temporal width with time are given in Figure 6.6 
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for STEREO- A (top) and STEREO-B (bottom). 

The temporal width of the pulse is shown in Figure 6.6 to increase with time 
for observations from both STEREO spacecraft, indicating that the pulse broad- 
ens in both space and time with propagation. The clear increase in spatial and 
temporal pulse width in data from both STEREO spacecraft strongly suggests 
that it is a true feature of the disturbance and not an observational artifact. The 
rate of change of the spatial and temporal width of the pulse with distance and 
time respectively (d(Ar)/dr and d(Ar)/dt) are given in Table 6.1. 

The observed broadening of the CBF pulse in both space and time confirms 
the previous observations of both Warmuth et al. (2001) and Veronig et al. (2010). 
These observations, when taken together with the derived kinematics, are consis- 
tent with the interpretation of a CBF as a dispersive, decelerating pulse. 

6.4.3 Integrated Intensity 

The temporal variation in the PA-averaged integrated intensity of the CBF pulse 
was also examined in addition to the pulse kinematics and the temporal and 
spatial dispersion of the pulse. This provides a verification of the presence of 
pulse dissipation, allowing the physical nature of the disturbance to be identified. 
The PA-averaged integrated intensity of the pulse can be determined using the 
variation in peak intensity with distance {(Ibd)^^ versus r) and the variation in 
the FWHM of the pulse with distance (Ar versus r). This allows the variation 
in PA-averaged integrated intensity to be estimated using, 

hot = Ar(/ B D)max, (6.3) 

and plotted as a function of distance. Figures 6.7 and 6.8 show the variation 
with distance in peak BD intensity (top), FWHM of the fitted pulse (middle) 
and resulting integrated intensity (bottom) for STEREO-k and STEREO-B re- 
spectively for the 2007 May 19 event. 

The variation in the integrated intensity rather than the peak intensity of the 
pulse was examined as the pulse has exhibited broadening (i.e., dispersion). This 
dispersion may influence the peak intensity of the pulse by causing it to decrease 
as the pulse broadens. In this case, the pulse would appear to show a drop in 
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Figure 6.7: Top: Variation in peak base difference intensity ({I bd) max) with dis- 
tance. Middle: Variation in FWHM of the Gaussian pulse with distance. Bottom: 
Variation in integrated intensity {I tot) with distance. All plots correspond to mea- 
surements from STEREO-A, with the dashed (dotted) line in all panels indicating 
the best linear fit to the 195 A (171 A) data. 
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Figure 6.8: Same as Figure 6.7 but showing the STEREO-A data for the event 
from 2007 May 19. 
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intensity but would contain the same amount of energy, giving the impression of 
dissipation without actually dissipating any energy. Similarly, the peak intensity 
of the BD intensity profile was used as this shows the actual emission of the pulse, 
while the PBD intensity shows the ratio of the pulse emission with respect to the 
background. As the coronal background is not uniform across the extent of the 
pulse, similar PBD intensity variations may not correspond to similar increases 
in emission. 

The top panel in both Figures 6.7 and 6.8 shows a large discrepancy between 
the peak BD intensity of the 171 A and 195 A observations. In both cases, the 
195 A data shows a large decrease with time while the 171 A data shows no 
strong variation. As previously shown in Section 6.4.2, the FWHM of the pulse 
follows a linearly increasing trend with distance for both passbands indicating 
pulse dispersion. However, the resulting variation in PA-averaged BD integrated 
intensity is inconclusive in each case shown in Figures 6.7 and 6.8. There does 
not appear to be any significant variation in 195 A measurements with distance, 
although some strong point-to-point variation is noticeable. The 171 A measure- 
ments however show a generally increasing trend in each case. The rate of change 
of I to t with distance (i.e., d(/ tot )/dr) for each event studied are given in Table 6.1. 

The measurements shown here indicate that more analysis is required to un- 
derstand the morphology of the pulse. The higher cadence observations available 
from SDO should allow the variation in peak intensity and PA-averaged inte- 
grated intensity with distance to be determined with a much higher degree of 
accuracy than is possible using STEREO. The variation between the different 
passbands observed by SDO should also allow a better understanding of the tem- 
perature structure of the pulse. 

6.5 Results 

The analysis described in the previous Section 6.4 indicates that the pulse iden- 
tification algorithm described in Sections 5.1.2 and 6.3 is extremely effective at 
identifying the CBF pulse in observations. This is particularly apparent in the 
intensity profiles shown in Figures 6.2 and 6.3, where the model provides a good 
fit to the data in each case. The errors associated with the identification of the 
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pulse are also quantified using this algorithm, unlike the unknown user-defined 
errors associated with the point-and-click approach. 

Once the source point and arc sector into which the pulse propagates have 
been identified, the algorithm provides accurate and reproducible estimates of 
the pulse characteristics, making these results more robust. The re-analysis of 
the CBF event of 2007 May 19 outlined in the previous section indicate that the 
pulse exhibited significant deceleration and dispersion. With the capabilities of 
the algorithm confirmed, it was then applied to three more events, from 2007 De- 
cember 07, 2009 February 12 and 2009 February 13. The results of these analyses 
are discussed in this section, with all of the results summarised in Table 6.1. 

6.5.1 Pulse Kinematics 

The intensity profiles of the events from 2007 December 07, 2009 February 12 and 
2009 February 13 were determined using the arc sectors shown in Figure 6.1, with 
the same orientation for the arc chosen for both STEREO spacecraft to ensure 
that the different results were comparable. The resulting intensity profiles for the 
different events studied are shown in Figures 6.9, 6.10, 6.11 and 6.12. A difference 
in the peak intensity is apparent between the 195 and 171 A passbands, particu- 
larly for the 2007 December 07 event. This is most likely due to the different peak 
emission temperatures observed by the two passbands, and could be explained by 
plasma heating or cooling out of the 171 A passband (cf. Wills-Davey & Thomp- 
son, 1999). It should also be noted that for the event of 2009 February 13, it 
was only possible to accurately detect the pulse in three consecutive images in 
the 171 A passband. After ~05:46:19 UT, the errors associated with the pulse 
identification became too large, with the result that only the first three detections 
were used in this analysis. 

As with the 2007 May 19 event discussed in Section 6.4, the kinematics were 
determined for each pulse studied using the temporal variation of the centroid of 
the Gaussian fit to the intensity profile. The pulses exhibited comparable initial 
propagation velocities ranging from ~249 km s" 1 for the 2007 December 07 event 
as observed by STEREO-B to ~406 km s" 1 for the 2009 February 12 event as 
observed by STEREO-B. The distance-time plots measured for each event are 
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Figure 6.9: Same as Figure 6.2 but for the event from 2007 Dec 07 observed by 
STEREO-k. 
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10: Same as Figure 6.2 but for the event from 2007 Dec 07 observed by 
B. 
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Figure 6.11: Same as Figure 6.2 but for the event from 2009 Feb 12 observed by 
STEREO-B. 
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Figure 6.12: Same as Figure 6.2 but for the event from 2009 Feb 13 observed by 
STEREO-B. 
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Figure 6.13: Same as Figure 6.4 but for the event on 2007 Dec 07. 



shown in Figures 6.13 and 6.14, with the kinematics values for each event shown 
in Table 6.1. 

While the kinematics of the 2007 May 19 event as observed by both STEREO 
spacecraft are very similar, the increased separation between the spacecraft is 
apparent in the kinematics estimated for the 2007 December 07 event, with the 
differences in initial velocity most likely due to geometrical effects. The quadra- 
ture events of 2009 February 12 and 13 were different from each other despite 
originating from the same active region. The 2009 February 12 event exhibited 
a higher initial velocity of ~406 km s" 1 , while the 2009 February 13 event had 
a lower initial velocity of ~273 km s -1 . The 2007 May 19 event appears to have 
been a relatively fast event, while the range of velocities (~240-450 km s" 1 ) is 
consistent with previously estimated values. 

All of the events studied exhibited statistically significant non-zero decelera- 
tion, although the strength of the deceleration varied between events. The de- 
celeration estimated for the 2007 December 07 event was also affected by the 
position of the STEREO spacecraft, varying from -153 ± 73 m s" 2 in STEREO- 
A data to -119 ± 109 m s" 2 in STEREO-B data. The 2009 February 12 was 
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Figure 6.14: Same as Figure 6.4 but for the events on 2009 Feb 12 (top) and 
2009 Feb 13 (bottom), both showing STEREO-B data. 



also observed to exhibit stronger deceleration than the 2009 February 13 event 
(—293 ± 167 m s~ 2 as opposed to —49 ±34 m s~ 2 ). The large errors associated 
with the acceleration terms given here indicate the difficulties associated with 
accurately determining the kinematics of CBFs from low cadence observations, 
despite the minimization of errors through the use of both the intensity profile 
technique and bootstrapping analysis. 



6.5.2 Pulse Broadening 

Next, the variation in both spatial and temporal pulse width with distance and 
time respectively was determined for each event studied. The resulting variation 
in spatial pulse width with distance is shown in Figures 6.15 and 6.16, while 
Figures 6.17 and 6.18 show the variation in temporal pulse width with time. 

It is clear from Figures 6.15 to 6.18 that each of the pulses studied exhibited 
evidence of pulse broadening in both spatial and temporal regimes with distance 
and time respectively. This suggests that pulse broadening is a general character- 
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Figure 6.15: Same as Figure 6.5 but showing the event on 2007 Dec 07. 
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Figure 6.16: Similar to Figure 6.5, now only showing STEREO-B data for the 
events on 2009 Feb 12 (top) and 2009 Feb 13 (bottom). 
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Figure 6.17: Same as Figure 6.6 but showing the event on 2007 Dec 07. 
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Figure 6.18: Similar to Figure 6.6, now only showing STEREO-B data for the 
events on 2009 Feb 12 (top) and 2009 Feb 13 (bottom). 
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istic of CBFs and must be accounted for by all theories that seek to explain this 
phenomenon. Table 6.1 shows the rate of change of both spatial and temporal 
pulse width Ar and At with distance and time respectively, with the d(Ar)/dr 
and d(Ar)/dt parameters positive for both the 171 A and 195 A passbands in all 
events studied. 

These results strongly imply that CBF pulses are dispersive pulses, something 
that may have been previously underestimated as a consequence of the techniques 
used to identify them. The presence of both pulse deceleration and dispersion is 
also suggestive of a wave interpretation for this phenomenon and indicates that 
the pulses are freely propagating. 

6.5.3 Pulse Dissipation 

The derived kinematics and dispersion of the pulse could then be used to calculate 
and examine the position-averaged integrated intensity of the pulse as discussed 
in Section 6.4.3. The resulting plots showing the variation in peak intensity 
(top), FWHM (middle) and PA-averaged integrated pulse intensity (bottom) 
with distance for each event studied are shown in Figures 6.19 (2007 Decem- 
ber 07 STEREO-k), 6.20 (2007 December 07 STEREO-B), 6.21 (2009 Febru- 
ary 12 STEREO-B) and 6.22 (2009 February 13 STEREO-B). 

The variation in PA-averaged integrated pulse intensity with distance is dif- 
ficult to interpret, as apparent from Figures 6.19 to 6.22. The PA-averaged inte- 
grated intensity is flat or decreasing in the 195 A passband for each event, but 
the 171 A measurements are inconclusive. 

For the 2007 December 07 event, the 195 A measurements as observed by 
STEREO-k appear flat, with a slight increase in the corresponding measure- 
ments for STEREO-B. On the other hand in the 171 A passband, an apparent 
increase can be observed in measurements from both STEREO spacecraft. A 
large separation is also evident between the 171 A and 195 A measurements for 
this event. 

The events from 2009 February 12 and 2009 February 13 are also inconclusive. 
While both the 171 A and 195 A passbands exhibit a decrease for the 2009 Febru- 
ary 12 event, the 2009 February 13 event shows an increase and decrease with 
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Figure 6.19: Same as Figure 6.7 but showing the STEREO-A data for the event 
of 2007 December 07. 
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Figure 6.20: Same as Figure 6.7 but showing the STEREO-B data for the event 
of 2007 December 07. 
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Figure 6.21: Similar to Figure 6.7 but showing the STEREO-B data for the event 
on 2009 February 12. 
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Figure 6.22: Similar to Figure 6.7 but showing the STEREO-B data for the event 
on 2009 February 13. 
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distance for the 171 A and 195 A data respectively. The rate of change of the PA- 
averaged integrated intensity with distance for each passband is given in Table 6.1 
for each event studied along with the associated error in each case. 

The inconsistency of these measurements between different events and indeed 
between the different passbands from the same spacecraft for the same event 
indicate that this property of CBFs is not currently well-constrained. This may 
be a consequence of the low observing cadence of the STEREO /EUVI instrument, 
but could also be affected by the inhomogeneous propagation of the pulse and 
possible temperature variations in the front during its evolution. The multi- 
passband, high-cadence observations afforded by the Solar Dynamic Observatory 
(SDO) should allow the true variation (if any) to be determined to a high degree 
of accuracy. 

6.6 Discussions and Conclusions 

A semi-automated technique designed to identify, track and analyse "EIT waves" 
in the solar corona was applied to observations of four CBF events using data 
from the EUVI telescopes onboard both STEREO spacecraft. The technique 
applies Gaussian fits to intensity profiles across the propagating CBFs, return- 
ing the position, width, and PA-averaged integrated intensity of the pulse for 
each observation. This approach is designed to minimise the errors typically as- 
sociated with the identification of CBF pulses, allowing a better understanding 
of the physical properties of the pulse such as the kinematics and rate of pulse 
dispersion. While this method is similar to CBF identification techniques previ- 
ously proposed by Warmuth et al. (2004b), Podladchikova & Berghmans (2005), 
Wills-Davey (2006) and Veronig et al. (2010), the combination of this approach 
with a statistically rigorous bootstrapping method and the high cadence obser- 
vations afforded by STEREO /EUVI, allow the errors typically encountered with 
the analysis of CBFs to be minimised. 

The four events used for this analysis (from 2007 May 19, 2007 December 07, 
2009 February 12 and 2009 February 13) were chosen as the erupting active region 
in each case was close to disk centre, minimising the projection effects of tracking 
a pulse across a spherical surface. The events were also well observed by the 
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STEREO spacecraft, with the close positioning of the STEREO spacecraft for 
the 2007 May 19 and 2007 December 07 events in particular allowing the pulse 
kinematics to be derived using both spacecraft. 

All of the pulses studied were found to exhibit statistically significant decel- 
eration ranging from ~ —300 to —50 m s -2 , with initial velocities ranging from 
~250-450 km s" 1 . This is comparable to previous estimates of CBF velocities, 
and is interesting given that most previous work has used point-and-click tech- 
niques applied to running-difference images. This approach identifies the forward 
edge of the CBF at a given time, using this to determine the kinematics of the dis- 
turbance as a whole. The analyses performed using such techniques have mainly 
returned kinematics that suggest a zero acceleration (i.e., constant velocity) in- 
terpretation of the CBF phenomenon. In contrast, our semi- automated technique 
uses the pulse centroid derived from percentage base difference intensity profiles 
to return a constant, non-zero acceleration (i.e., variable velocity) interpretation. 

This discrepancy can be understood by considering the physical processes at 
work in the propagating pulse. If the pulse is broadening with propagation, as 
found here, the forward edge of a decelerating pulse will appear to move faster 
than the pulse centroid. If the rate of pulse broadening is comparable to the 
deceleration of the pulse, the deceleration may be masked by this broadening. 
This suggests that the true kinematics of CBFs may have previously been dis- 
guised through use of the wrong position within the pulse profile to characterise 
its location. It should be noted that a variable acceleration may be present, but 
this can only be studied with an increased number of data points. 

A positive increase in the variation of the width of the Gaussian fit to the 
CBF pulse in both the spatial and temporal domains was observed for all events 
studied. Although the increase in spatial pulse width is immediately suggestive 
of a dispersive pulse, this property is confirmed by the observed increase in the 
temporal width of the pulse. As shown in Equation 6.2, an increase in the spatial 
pulse width may be negated by an increase in pulse velocity, producing a pulse 
with a constant temporal width. It is apparent from Figures 6.5, 6.6, and 6.15 
to 6.18 that this is not the case and the observed CBF pulses are dispersive in 
nature. Although some previous studies have suggested a dispersive nature for 
"EIT waves" (e.g., Warmuth et al, 2004a), the true extent of this dispersion may 
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have been disguised by the use of point-and-click techniques to analyse running- 
difference images. 

The confirmation of pulse broadening resulted in the variation in PA-averaged 
integrated pulse intensity with distance being studied rather than the variation 
in peak pulse intensity. Although the variation in the peak intensity of the pulse 
with distance is generally considered an indicator of the pulse energy, the dis- 
persion can cause the peak of the pulse to drop without reducing the energy 
contained within it. The PA-averaged integrated pulse intensity was found by 
multiplying the FWHM by the peak pulse intensity at a given time. This was 
found to produce inconclusive results between different events and also between 
different passbands for the same event. While the peak intensity of the pulse 
was generally observed to decrease with distance and the FWHM was observed 
to increase with distance in both passbands, the PA-averaged integrated inten- 
sity typically showed no systematically strong variation with distance in either 
passband, although a strong offset between the 171 A and 195 A data was typ- 
ically observed. Contrasting results were found for the 2009 February 12 and 
2009 February 13 events respectively, suggesting that further investigation us- 
ing simultaneous analysis of multiple passbands at very high cadence is required, 
something that will be routinely available from the SDO spacecraft. 

The analysis outlined in this chapter indicates that the studied CBFs may 
be viewed as dispersive pulses exhibiting negative acceleration. This interpreta- 
tion is consistent with the fast-mode magnetoacoustic wave interpretation for a 
freely-propagating "EIT wave". However, the inconclusive variation in the inte- 
grated intensity of the pulse indicates that more analysis is required to definitively 
determine the physical nature of the disturbance. 

Although only four CBF events were considered here, the consistency of the 
results between the events studied suggests that the conclusions drawn here may 
be applicable to a larger sample of CBFs. The initial velocities of the CBFs 
are comparable to the lower range of estimated Alfven speeds proposed by Wills- 
Davey et al. (2007), suggesting that the randomly structured nature of the coronal 
magnetic field may have an important effect on the propagation of CBFs; an 
effect shown in simulations performed by Murawski et al. (2001). Although the 
kinematics and dispersive nature of the pulses suggest that a magnetoacoustic 
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wave interpretation may be correct, further study is required. The ~10 s temporal 
cadence across multiple EUV passbands available with the launch of the SDO 
spacecraft will allow these results to be studied in much greater detail. 
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Chapter 7 



The Wave Properties of Coronal 
Bright Fronts Observed Using 



In a dispersive medium, the components of a wave group move with different 
speeds, and the phase relations between the components are altered 



Coronal bright fronts are large scale wavefronts that propagate through the 
solar corona at hundreds of kilometers per second. While their kinematics have 
been studied in detail, many questions remain regarding the temporal evolution of 
their amplitude and pulse width. In this chapter, contemporaneous high cadence, 
multi-thermal observations of the solar corona from the SDO and STEREO space- 
craft are used to determine the kinematics and expansion rate of a CBF pulse 
observed on 2010 August 14. 

The CBF was found to have a lower initial velocity with weaker deceleration 
in STEREO observations compared to SDO (~340 km s -1 and -72 m s -2 as 
opposed to ~410 km s" 1 and —279 m s -2 ; cf. Chapter 4). The CBF kinematics 
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from SDO were found to be highly passband-dependent, with an initial velocity 
ranging from 379 ± 12 km s _1 to 460 ± 28 km s _1 and acceleration ranging from 
-128±28 m s" 2 to -431±86 m s~ 2 in the 335 A and 304 A passbands respectively 
These kinematics were used to estimate a quiet coronal magnetic field strength 
range of ~l-2 G. Significant pulse broadening was also observed, with expansion 
rates of -130 km s" 1 (STEREO) and -140 km s" 1 (SDO). By treating the CBF 
as a linear superposition of sinusoidal waves within a Gaussian envelope, the 
resulting dispersion rate of the pulse was found to be ~10-15 Mm 2 s -1 . These 
results are indicative of a fast-mode magnetoacoustic wave pulse propagating 
through an inhomogeneous medium. This chapter contains work by Long, DeLuca 
& Gallagher (2011b) published in The Astrophysical Journal Letters. 

7.1 Introduction 

Although CBFs were first noted by Moses et al. (1997) and characterised by 
Thompson et al. (1998), they remain a source of much debate with many different 
theories proposed to explain the pheomenon. This is partly a consequence of 
their rare nature but is also due to the instruments used to identify and study 
them, with SOHO (Thompson et al, 1999), TRACE (Wills-Davey k Thompson, 
1999) and STEREO (Long et al, 2008) offering a choice between low cadence 
observations with a large field of view (SOHO and STEREO) or high cadence 
observations with a small field of view (TRACE). 

Despite these restrictions, CBFs have been relatively well characterised. They 
are usually observed as diffuse bright fronts that propagate isotropically when 
unimpeded at typical velocities of 200-400 km s _1 across the solar disk (Thomp- 
son & Myers, 2009). Optimum observations occur using the 195 A (Fe XII) 
passband, corresponding to a temperature of ~l-2 MK and a height of ~70- 
90 Mm above the photosphere (i.e., the low corona; Patsourakos & Vourlidas, 
2009; Kienreich et al, 2009). However, they have been observed using other 
coronal EUV passbands including the 171 A (Fe IX; Wills-Davey & Thompson, 
1999), 284 A (Fe XV; Zhukov & Auchere, 2004) and 304 A (He II; Long et al, 
2008) passbands. 

Although the observations of CBFs have been mostly consistent, they have 
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given rise to a multitude of theoretical interpretations. They have been alterna- 
tively interpreted as magnetohydrodynamic waves (Wang, 2000; Warmuth et al, 
2004b; Wang et al, 2009; Schmidt & Ofman, 2010), solitons (Wills-Davey et al, 
2007) and in terms of magnetic field restructuring during the eruption of an as- 
sociated CME (Chen et al, 2002; Attrill et al, 2007a; Delannee et al, 2008). 
This proliferation of theories for interpreting CBFs is a direct consequence of the 
traditional techniques used to study the phenomenon, which produce strongly 
user-defined conclusions. This has resulted in different authors deriving differ- 
ent results for the same event (e.g., Thompson et al, 1998; Attrill et al, 2007a; 
Delannee et al, 2008). 

This combination of limited identification techniques with the low cadence of 
previous observations has produced results that imply a wave interpretation (such 
as pulse reflection and refraction; Veronig et al, 2006; Gopalswamy et al, 2009), 
but with velocities that are lower than those predicted by this interpretation. 
However, observations of decelerating CBFs combined with the effects of low ob- 
serving cadence (Long et al, 2008; Ma et al, 2009; Long et al, 2011a) suggest 
that the low observing cadence and techniques used to identify the pulses may 
influence the derived kinematics. There have also been indications of CBF dis- 
persion with propagation (Warmuth et al, 2004b; Long et al, 2011a), although 
this has been difficult to quantify. While these properties are inconsistent with 
ideal MHD wave theory, they have been shown in simulations by Murawski et al. 
(2001) and Nakariakov et al. (2005) to be a natural result of propagation through 
an inhomogeneous medium. These observations of deceleration, dispersion and 
dissipation are also consistent with the interpretation of CBFs as large-amplitude 
waves as a direct consequence of their inherently non-linear nature (Vrsnak & 
Cliver, 2008; Vrsnak & Lulic, 2000a; Warmuth, 2010). In this case, the rapid 
expansion of an erupting CME or rising flare loops induce a large-amplitude 
disturbance that propagates as a simple wave before steepening and forming a 
discontinuity. As the shock propagates, its shocked edge moves faster than the 
trailing edge, while the amplitude drops as a consequence of geometrical expan- 
sion before the shock ultimately decays to an ordinary small-amplitude pulse. A 
more complete description of the different observations of CBFs and their result- 
ing theoretical interpretations can be found in Chapter 2, as well as in the many 
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recent reviews of the field by Gallagher & Long (2011), Zhukov (2011), Warmuth 
(2010) and Wills-Davey & Attrill (2010). 

The launch of the Solar Dynamics Observatory (SDO) spacecraft in 2011 
marked the beginning of a new age of solar physics, with the Sun now visible 
in more detail than ever before. SDO/AIA observes the Sun continuously at a 
cadence of ~12 s in seven EUV passbands, an improvement on both SOHO /EIT 
(~900 s in one of four passbands) and STEREO /EU VI (~75-600 s in four pass- 
bands). As well as the 171 A, 193 A 1 and 304 A passbands observed by SOHO 
and STEREO, SDO also utilises the 94 A (Fe XVIII), 131 A (Fe VIII), 211 A 
(Fe XIV) and 335 A (Fe XVI) passbands, allowing the Sun to be studied across 
a wide range of temperatures and heights. 

The first observations of a CBF using SDO were made by Liu et al. (2010), 
who noted the presence of the pulse in all EUV passbands monitored by the 
spacecraft, although the pulse was only analysed in detail using the 171 A and 
193 A passbands. This was partly due to the sheer volume of data available 
from SDO (~1.5 TB per day) which is not conducive to manual data analysis. 
This data volume has necessitated the development of both automated and semi- 
automated CBF detection and tracking algorithms (Podladchikova & Berghmans, 
2005; Wills-Davey, 2006; Long et al, 2011a) to allow clarification of the true 
nature of CBFs. 

In this chapter, the semi-automated CBF algorithm discussed in Sections 5.1.2 
and 6.3, and published by Long et al. (2011a) is applied to SDO and STEREO 
observations of the 2010 August 14 CBF event. The observations and data used 
for this analysis are outlined in Section 7.2. This allowed the physical nature of 
the pulse to be studied in detail, with the kinematics of the pulse and variation 
in pulse width discussed in Sections 7.3.1 and 7.3.2. The multi-thermal nature 
of SDO allowed the temperature dependence of the pulse to be examined; this 
discussion is presented in Section 7.3.3. The variation in pulse kinematics, pulse 
width and the multi-thermal nature of the pulse were then combined to study the 
local plasma through which the pulse propagated, allowing the coronal magnetic 

lr The Fe XII passband peaks at 195 A in STEREO /EUVI, but at 193 A in SDO/AIA. As 
a result, all references to the 193 A passband will refer to SDO data, with the 195 A passband 
referring to STEREO data. 
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field strength to be estimated (Section 7.3.4). Finally, some conclusions regarding 
the implications of this work are discussed in Section 7.4. 

7.2 Observations and Data Analysis 

The CBF event observed on 2010 Aug 14 1 erupted from NOAA active region (AR) 
11093, with an associated coronal mass ejection (CME) and GOES C4.4 flare 
which started at 09:38 UT. The position of the active region (N11W65) meant 
that the on-disk CBF evolution was visible from STEREO-k and SDO but not 
STEREO-B due to the relative positions of the STEREO spacecraft. At the time 
of the eruption, STEREO /EUVI-A was operating with an observing cadence of 
300 s and 600 s in the 195 A and 304 A passbands respectively, although the 
171 A and 284 A passbands were both taking synoptic data (i.e., one image every 
two hours). In contrast, SDO/AIA was taking observations with 12 s cadence in 
all seven EUV passbands (94, 131, 171, 193, 211, 304, and 335 A) over the same 
time period. 

As previously discussed in Sections 5.1.2 and 6.3, the semi-automated de- 
tection algorithm used to identify and track the CBF has several steps. The 
algorithm was applied to data from both the STEREO and SDO spacecraft us- 
ing an arc-sector of 30° width in both cases, with the routines used to read in 
and prepare the data the only modifications made to the algorithm. This was 
necessary as the different data require different processing to achieve scientific 
quality (see Sections 3.1.3 and 3.2.3 for more details). 

The source point for the CBF was defined using the mean centre of ellipses 
fitted to the first three observations of the CBF in both 193 A (AIA) and 195 A 
(EUVI) data, giving a source unique to each spacecraft (although both sources are 
comparable when transformed between spacecraft). Percentage base difference 
images were used for this analysis, with the images de-rotated in each case to the 
same pre-event time (09:20:30 UT) to compensate for solar rotation. The base 
image for each passband was defined using a pre-event time of ~09:25:00 UT, with 
a sample of the resulting images shown in Figure 7.1. The arc sector employed by 
the algorithm to identify the pulse was then positioned to allow direct comparison 

1 Solar Object Locator: SOL2010-08-14T09:38:00L353C79 
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of both AIA and EUVI observations and to ensure the maximum number of 
observations using both instruments. 

The operation of the algorithm is similar to that previously described, with 
the PBD intensity of a given image averaged across the position angle of the 
arc sector in annuli of increasing radii, with the standard deviation giving the 
associated error. The width of the annuli used varied between the spacecraft as 
a consequence of the different spatial resolution available with 1 degree annuli 
used for STEREO and 0.5 degree annuli used for SDO. This had the effect of 
increasing the detail of the resulting intensity profiles from SDO, but did not 
affect the results of the analysis. 

The resulting intensity profile was then fitted using a Gaussian function, with 
the centroid and full width at half maximum (FWHM) giving the pulse position 
and width respectively. Each parameter was returned with an associated error 
allowing the ability of the algorithm to detect the pulse to be quantified. Although 
the source point position and orientation of the arc sector is determined by the 
user, all subsequent detection, tracking and analysis of the pulse is automated, 
allowing unbiased identification of the CBF. Once the intensity profiles for each 
image have been produced and fitted, the CBF is identified as a moving pulse, 
with any stationary bright features ignored. 
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7.3 Results 



The actual identification of the pulse was found to be strongly influenced by 
the passband used rather than instrument, with some evidence of a disturbance 
visible in all available passbands. The pulse was observed in the 195 A and 304 A 
passbands from STEREO /EUVI, although only the 195 A data was used here due 
to the low cadence of the 304 A passband. In SDO /AIA, it was possible to track 
the pulse in four of the seven passbands (193, 211, 304, and 335 A). The nature 
of the 94 A and 131 A passbands complicated identification of the pulse, while a 
slight intensity decrease was visually identified in the 171 A passband. However, 
it was not possible to track this using the detection algorithm. A sample of the 
resulting intensity profiles from the four passbands used in this work is shown in 
Figure 7.2. 

7.3.1 Pulse Kinematics 

The temporal variation in the distance of the fitted pulse centroid from the defined 
source point was used to determine the kinematics of the pulse using measure- 
ments from all passbands studied. This was done by using the bootstrapping 
technique described in Section 5.2.3 to fit a model of the form, 

r(t) = r + v t + l -at 2 (7.1) 

to the distance-time measurements. Here r describes the initial distance of 
the pulse from the source point, v is the initial velocity and a is the constant 
acceleration. 

The temporal variation in the distance measurements for each of the four SDO 
passbands studied is shown in the top panel of Figure 7.3, while the STEREO 
195 A and SDO 193 A measurements are compared in the bottom panel of the 
same figure. The kinematics of the individual SDO passbands are seen to be 
comparable, although there does appear to be some separation with propagation. 
The apparent discontinuities in individual SDO passbands seen in Figure 7.3 
are the result of a stationary bright feature close behind the detected pulse which 
influenced the detection of the pulse by the algorithm, while there was a data-gap 
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Figure 7.2: Variation in SDO percentage intensity profiles with Great Circle 
distance (Mm) in the (from left to right) 193, 211, 335 and 304 A passbands, with 
the Gaussian fit applied by the algorithm to the data in each case highlighted in 
red. Time is increasing from top to bottom, with the time of each intensity profile 
given in the upper right of each panel. 



in the 304 A passband which also complicated detection of the pulse. These effects 
were minimal for the distance-time measurements but exhibited a much stronger 
influence on the variation in pulse width with time as discussed in Section 7.3.2. 
The measurements from the similar 193 A and 195 A passbands shown in the 
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Figure 7.3: Top: Distance-time measurements from AIA 304, 193, 211, and 335 A 
passbands for the 2010 August 14 CBF event. Bottom: AIA 193 A and EUVI 195 A 
distance-time measurements for the same event (line shows the best fit to EUVI 
measurements). 

bottom panel of Figure 7.3 also appear homologous, with a slight positional offset 
due to the different spacecraft positions. 

The estimated kinematics derived using the bootstrapping technique are pre- 
sented in Table 7.1. Each of the passbands studied exhibits statistically signifi- 
cant negative acceleration, with the higher cadence of the SDO data contributing 
to the lower associated errors for each parameter. A comparison of the simi- 
lar STEREO 195 A and SDO 193 A passbands shows that the 195 A passband 
exhibits a lower initial velocity and much weaker acceleration relative to the 
comparable 193 A passband. The kinematic estimates from STEREO 195 A are 
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consistent with previous results derived using SOHO/EYT and STEREO /EUVI, 
while the higher initial velocity and acceleration from SDO 193 A for the same 
event suggests a strong influence from the cadence of the observing instrument, 
confirming the results of Long et al. (2008) and Ma et al. (2009). 

The passbands studied by SDO exhibit consistently higher initial velocities 
and stronger negative accelerations than previously estimated for CBF pulses, 
suggesting that the pulses may be consistent with MHD waves. The larger un- 
certainties associated with the kinematics analysis of the 304 A passband may be 
explained by the nature of the passband and also by a data gap, which compli- 
cated detection of the pulse. The discrepancy between the kinematics derived for 
each SDO passband is interesting and will be discussed further in Section 7.3.3. 

7.3.2 Pulse Broadening 

The temporal variation in the FWHM of the fit to the intensity profile was next 
examined for evidence of pulse broadening. This variation is shown in the top 
panel of Figure 7.4, which indicates that the pulse width changes from ~70 Mm 
to ~270 Mm over a time period of ~900 s. This increase was observed for all 
passbands studied, with the STEREO 195 A measurements also showing a similar 
increase, albeit over a longer timescale. 

The data prior to ~09:52 UT (to the left of the dashed vertical line) has been 
corrected to remove the effects of a stationary bright feature close behind the 
CBF (cf. Muhr et al, 2011). This feature initially exerts a strong influence on 
the Gaussian fit, but was negated by subtracting a constant offset value for each 
passband from the FWHM measurements. From ~09:52 UT onward, the CBF is 
sufficiently far from this feature that the fit to the data is no longer affected. The 
effects of this bright feature can also be seen in the bottom panel in Figure 7.4, 
which shows the peak % intensity variation with distance. While the bright 
feature does initially influence the pulse width and peak intensity variation, a 
general increase and decrease is apparent for the pulse width and peak intensity 
respectively. 

To try to understand the physical nature of the observed pulse dispersion, the 
CBF was treated as a linear superposition of sinusoidal waves within a Gaussian 
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Figure 7.4: Top: Variation in FWHM with time for AIA passbands; EUVI 195 A 
measurements show a similarly increasing trend over a much longer time range. 
The black line shows the fit to the combined SDO data. Bottom: Peak PBD pulse 
intensity variation with distance. 



envelope. This can be represented using the equation, 

V>(r, f) ~ exp (- ^2^ ) cos ( k » r ~ u ot) (7.2) 

where ko is the wavenumber, uo is the angular frequency, ay is the characteristic 
width and v g is the group velocity of the pulse (v g = du/dk). The pulse extends 
in Fourier space from k — Ak/2 to k + Ah/2 (Ak ~ a k where a k = l/cr r ), so 
that the velocity varies from v g (k — Ak/2) to v g (k + Ak/2) across the pulse. As 
the pulse propagates, this velocity gradient produces an increase in the spatial 
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extent of the pulse, where the spatial extent (FWHM) is defined as, 



Ar(t) = Ar + 



, , AAA / Ak 
v g[ k o+ ^ J - v g \k - — 



t (7.3) 



and Ar is the initial pulse width. This can be rewritten in terms of the change 
in group velocity v g as, 

Ar(t) ~ Ar + ^^-Akt. (7.4) 

As the group velocity v g = du/dk, this can be redefined to give the width of a 
dispersive pulse at any time t as, 

* / x a d 2 u(k ) t 

Ar( t) = Ar 0+ -^— (7.5) 

where d 2 u(ko)/dk 2 is the rate of change of the group velocity of the pulse with 
respect to k. Equation 7.5 was fitted to the FWHM measurements given in the 
top panel of Figure 7.4, allowing d 2 u(ko)/dk 2 to be determined for each passband. 

The rate of expansion of the pulse Ar(t) and the resulting value of d 2 u(ko) / dk 2 
determined using the top panel of Figure 7.4 are given in Table 7.1 for each 
passband studied. It is clear that the general expansion rate in each case is 
positive within error, indicating statistically significant pulse broadening. This 
implies that CBFs are dispersive pulses, confirming the results of Warmuth et al. 
(2004b); Veronig et al. (2010) and Long et al. (2011a). 



7.3.3 Temperature Dependence 

As a result of the very high cadence of SDO across multiple passbands at varying 
peak emission temperatures, it was possible to derive the accurate kinematics of a 
CBF for each passband studied. These values are given in Table 7.1, with a spread 
in the kinematics estimated for the SDO passbands, from ~380 to ~460 km s _1 
and ~ —128 to ~ —430 m s~ 2 for the initial velocity and acceleration respec- 
tively. Although this variation could be ascribed to measurement errors, there 
is a statistically significant difference between the derived values, indicating that 
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this discrepancy is real. This variation was examined by making a comparison 
with the peak emission temperatures (T pea k) of the different AIA passbands (as 
given in Table 7.1, Figure 3.8 and discussed by O'Dwyer et al, 2010). 

An inverse relationship between the kinematics of the CBF and the peak 
emission temperature T pea k for each passband was noted. As shown in Figure 1.3, 
the general trend in the quiet corona is for temperature to increase with height 
while density and magnetic field strength decrease. This implies that in cooler, 
more dense plasma the CBF has a higher velocity, although it should be noted 
that this is a general trend as plasma elements of different temperatures may co- 
exist at the same height in the corona. While this result appears unusual at first 
glance, it is characteristic of a compressive pulse and combined with the observed 
pulse dispersion and deceleration indicates that the CBF is best described as a 
magnetohydrodynamic wave pulse. The randomly structured nature of the quiet 
corona suggests that any globally-propagating pulse must traverse magnetic field 
lines, indicating a fast-mode rather than slow-mode interpretation for the CBF. 

However, the morphological signatures of the CBF across different passbands 
show some discrepancies that invite further investigation, particularly the simul- 
taneous intensity decrease at 171 A and increase in the cooler 304 A passband. 
While an intensity decrease in the 171 A passband has previously been noted 
for TRACE (Wills-Davey & Thompson, 1999) and SDO (Liu et al, 2010), this 
was not apparent in STEREO data due to the broad nature of the temperature 
response curve for this passband (see e.g., Raftery et al, 2009). This drop in 
171 A emission (visually identified here but not tracked) has previously been 
characterised as evidence of plasma heating from 171 A into the 193, 211, and 
335 A passbands (Wills-Davey & Thompson, 1999; Liu et al, 2010). This heat- 
ing implies that the CBF pulse is coronal, an observation consistent with the 
height measurements made by Kienreich et al. (2009) and Patsourakos & Vourl- 
idas (2009). 

However, this implication is complicated by the observed increase in 304 A 
emission, as this passband is dominated by two chromospheric He II lines at 
303.781 and 303.786 A. Although there is also a coronal Si XI emission line at 
303.33 A (cf. Section 4.5), O'Dwyer et al. (2010) have noted that this line does 
not make a significant contribution to AIA quiet Sun observations. This suggests 
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that the observed intensity increase is most likely due to He II emission. 

The formation mechanism of He II emission has been the subject of detailed 
investigation over many years (see e.g., MacPherson & Jordan, 1999; Andretta 
et al, 2003; Jordan & Brosius, 2007) as a result of its complex nature. Analysis 
has suggested that it is formed by collisional excitation resulting from thermal 
electrons in the quiet corona. This thermal excitation effect could be enhanced by 
the increased temperature gradient caused by passage of a compressive pulse in 
the solar corona, producing the observed 304 A intensity increase. Consequently, 
the CBF would be a coronal pulse, consistent with EUV observations, height 
measurements and the observed drop in 171 A intensity. 

7.3.4 Coronal Seismology 

The passband-dependent kinematics and the observed temporal variation in the 
width of the CBF pulse indicate that the pulse morphology is significantly influ- 
enced by the plasma through which it propagates. While the variation in pulse 
morphology could be used to examine the physical nature of the pulse itself, an 
examination of how the plasma affects the kinematics for each passband can be 
used to directly quantify the characteristics of the quiet coronal plasma. For 
example, as noted in Section 1.3.3, the fast-mode wave speed is defined as, 

Vfm = \jv 2 A + C 2 s1 (7.6) 

where the Alfven speed and sound speed are Va = B / (Annm) 1 / 2 and c s = 
(7/cT/m) 1 / 2 respectively. Here B is the magnetic field strength, n is the par- 
ticle density, m is the proton mass, 7 is the adiabatic index (typically 5/3), k is 
Boltzmann's constant and T is the peak emission temperature (T peak ; the values 
used are given in Table 7.1). 

The pulse kinematics and dispersion are consistent with the properties of a 
fast-mode wave, allowing the CBF to be treated as such. In this case, the final 
pulse velocity (i.e.: the velocity of the pulse when it can no longer be detected 
by the algorithm) can be taken as the fast-mode velocity of the given passband, 
since the pulse can not propagate below this velocity. These values are given in 
Table 7.1 for each SDO passband studied. 
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By taking the peak emission temperature of each passband as the temperature, 
only the magnetic field strength and density are unknown in the above equations 
for the sound speed and Alfven speed. Coronal magnetic field strength estimates 
typically involve extrapolating photospheric magnetic field measurements into the 
corona and are not very well constrained (particularly in the quiet Sun). However, 
coronal densities can be estimated using density sensitive line ratios (Gallagher 
et al, 1999) and are well-constrained. 

The above equations can be rearranged in terms of the magnetic field strength 
B as, 



allowing the quiet coronal magnetic field strength to be estimated using the de- 
rived CBF kinematics. This was done by combining the final velocity values given 
in Table 7.1 with a range of typical quiet coronal densities (~2-6xl0 8 cm -3 ; see 
Wills-Davey et al, 2007, for more details). This produced an estimated magnetic 
field strength range of ~l-2 G for the quiet solar corona. 

This estimated range of values is comparable to the value previously derived by 
West et al. (2011). In this case, West et al. (2011) used detailed STEREO /EUVI 
estimates of the kinematics of a CBF pulse and combined them with accurate 
measurements of the density of the corona through which the pulse propagated 
made using the Extreme ultraviolet Imaging Spectrometer (EIS) instrument on- 
board the Hinode spacecraft. The value of 1-2 G is also comparable to the value 
of ~1.5-3.5 G observed by Kienreich et al. (2011), who used four consecutive 
CBF eruptions originating from the same active region over a period of eight 
hours and observed by STEREO-B. By determining the velocity of each pulse, 
Kienreich et al. (2011) were able to derive the magnetosonic Mach number for 
each pulse and estimate the coronal magnetic field strength using a Saito coronal 
density model. The good agreement of the range estimated here with the work 
of West et al. (2011) and Kienreich et al. (2011) indicates that the assumptions 
made here are correct and CBFs can be used to probe the physical characteristics 
of the plasma through which they propagate. 




(7.7) 
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7.4 Discussion and Conclusions 



The CBF detection algorithm outlined by Long et al. (2011a) was applied to 
contemporaneous EUVI and high cadence AIA observations of the 2010 August 14 
CBF event. This allowed the accuracy of previous CBF kinematics estimates to 
be examined and was used to probe the physical properties of the pulse and the 
plasma through which it propagated. A comparison of pulse kinematics estimated 
using the same passband but at different cadences indicated that the estimates 
of the pulse kinematics are influenced by the observing cadence, supporting the 
conclusion of Long et al. (2008). Significant pulse deceleration was measured 
for this event in both EUVI and AIA data (despite the different cadence and 
spacecraft positions), suggesting that it is characteristic of the phenomenon. 

Analysis of the high cadence observations available for this event also sug- 
gested that the techniques previously used to estimate pulse kinematics may have 
been flawed. These techniques typically involved combining distance-time mea- 
surements from different passbands to calculate the general pulse kinematics for 
an event due to a paucity of data (e.g., Long et al, 2011a; Patsourakos & Vourl- 
idas, 2009; Kienreich et al, 2009; Veronig et al, 2010). While this was necessary 
to derive kinematics from the small data-sets available, the results presented here 
indicate that this approach underestimated the kinematics of the CBF. It may 
have also masked the amount of acceleration experienced by the pulse, and did 
not detail the effect of the plasma on the pulse. 

In addition to the kinematics derived for the pulse, the algorithm used also 
allowed the temporal variation in pulse width to be examined using both AIA 
and EUVI data. Clear dispersion was observed in both sets of data, confirming 
the observations of Warmuth et al. (2004b); Veronig et al. (2010) and Long et al. 
(2011a). These repeated measurements of significant pulse broadening strongly 
indicate that CBFs have a dispersive nature which, allied to the traditional point- 
and-click techniques for identifying them, may have contributed to the uncertainty 
surrounding their acceleration. 

When both the dispersion and deceleration are considered, CBFs may be 
best described using a wave interpretation. Although ideal linear MHD wave 
theory predicts a constant velocity pulse of constant width, deceleration and 
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dispersion are consistent with the propagation of a pulse through a randomly 
structured medium. This behaviour has been observed in simulations performed 
by Murawski et al. (2001) and Nakariakov et al. (2005), and is congruous with the 
randomly structured nature of the quiet corona. The observed deceleration and 
dispersion are also consistent with the interpretation of the CBF pulse as a large 
amplitude simple wave (cf. Vrsnak & Lulic, 2000a; Vrsnak & Oliver, 2008). In 
this case, the pulse is initially driven by the rapid expansion of the erupting CME 
or by rising flare loops, producing a non-linear pulse which steepens to form a 
shock front. Through geometric expansion, this shock front exhibits deceleration, 
dispersion and dissipation with propagation (see e.g., Muhr et al, 2011). 

The observed dispersion also allowed an insight into the physical nature of 
the pulse. This was possible by treating the pulse as a linear superposition of 
sinusoidal waves within a Gaussian envelope. The velocity gradient across the 
pulse could then be quantified and combined with the measured variation in pulse 
width to determine the dispersion relation of the pulse. This is the first time that 
the dispersion relation of a CBF pulse has been determined and provides a unique 
tool for the analysis of CBFs and the solar corona. 

The high cadence observations available from SDO allowed the kinematics of 
the pulse to be determined in each of the four passbands studied. Statistically 
significant differences were noted for both the initial velocity and acceleration esti- 
mated from each passband; a unique result that supports the wave interpretation 
of CBFs. In particular, the pulse exhibited a compressive nature, appearing to 
propagate at a higher velocity with stronger deceleration in cooler, denser plasma. 
This is the first time that this property of CBFs has been observed, and is a direct 
consequence of the very high observing cadence capabilities of SDO/AIA. 

The kinematic variation indicated that the propagation of the pulse was 
strongly affected by the plasma through which it propagated, and provided a 
simple diagnostic of the emitting plasma in each passband. This variation was 
used to determine the physical properties of the plasma through the application 
of coronal seismology (e.g., Ballai & Douglas, 2008). Although this approach has 
previously been proposed as a way of directly probing the structure of the solar 
corona, it has been complicated by the uncertain physical nature of CBFs. The 
results presented here, in addition to recent work by Patsourakos et al. (2010); 
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Kienreich et al. (2011) and Long et al. (2011a) strongly indicate that CBFs are 
fast-mode MHD waves, allowing them to be used to examine the environment 
through which they propagate. 

This analysis relies on the ability to determine typical coronal densities and the 
temperature of the different passbands using alternative techniques. With these 
values known or well-constrained, the magnetic field strength can be estimated 
using the CBF. This was done here for a range of typical coronal density estimates, 
producing an estimated magnetic field strength range of ~l-2 G for the quiet 
Sun. This range of values is comparable to the value estimated by West et al. 
(2011) (and typically assumed for the quiet corona, e.g., Wills-Davey et al, 2007), 
indicating that CBFs can be used to directly probe the plasma through which 
they propagate. 

The results of this analysis are most compatible with the wave interpretation 
of a CBF pulse. The observed dispersion implies that CBFs are not accurately 
described by the soliton model proposed by Wills-Davey et al. (2007), while the 
CBF height range (on-disk near the limb over an extended time period in both 
SDO and STEREO observations) is inconsistent with the progressively higher 
emission predicted by Delannee et al. (2008). The multi-temperature emission 
does not match the low foot-point signature predicted by Attrill et al. (2007a) 
and there was no indication of the additional coronal Moreton wave predicted by 
Chen et al. (2002). 

Despite this, the initial driver producing the CBF remains uncertain and the 
subject of some debate. However, the CBF could be the result of the rapid over- 
expansion of the erupting CME bubble (cf. Patsourakos et al, 2010) low down 
in the solar corona before decoupling and propagating freely. The high cadence 
observations available from SDO will allow this issue to be resolved. 
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Chapter 8 

Conclusions and Future work 



What's next? 

Aaron Sorkin (The West Wing) 



Although "EIT waves" have been studied for almost 15 years, the relatively 
low observing cadence available from SOHO/EYT has traditionally made analysis 
difficult, and has necessitated combination with observations from other instru- 
ments and passbands. While this was and is necessary for a thorough analysis, 
it has resulted in a widely varying description and interpretation of this phe- 
nomenon. Here, data from the STEREO and more recently SDO spacecraft have 
been used to study coronal bright fronts in a more statistically rigorous manner 
and across a wider range of EUV passbands than previously possible. 

The quality of data such as that available from STEREO and SDO has re- 
sulted in a deeper understanding of the inherent physics of the "EIT wave". 
It has been shown that CBF pulses exhibit statistically significant deceleration 
with propagation, while pulse broadening has also been displayed. This pulse 
broadening combined with the techniques traditionally used to study CBFs may 
have contributed to the uncertainty regarding the true kinematics of these dis- 
turbances. 
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The higher time cadence available from more recent spacecraft has also been 
shown to influence the estimated kinematics of the CBF pulses. Although this 
was suggested by the initial observations of Long et al. (2008), the contempo- 
raneous CBF observations from STEREO and SDO (with its higher cadence) 
of the 2010 Aug 14 event have shown that this effect must be accounted for 
when analysing this phenomenon. The availability of higher cadence observations 
from SDO is also necessitating the development of more rigorous, automated ap- 
proaches to the identification of individual disturbances. The technique outlined 
here has been shown to provide accurate determinations of pulse kinematics, 
allowing the physical nature of individual pulses to be examined. 

8.1 Principal Results 

The original aim of this work was to study the "EIT wave" phenomenon, with 
a view to determining its true nature using purely observable basic physics. Al- 
though the traditional theoretical interpretation of this phenomenon is complex 
and divided between alternative point s-of-view, this work was driven by obser- 
vation, with the launch of the STEREO spacecraft at the start of this project 
allowing a new perspective on these disturbances. It should be possible to distin- 
guish between different interpretations by determining and examining the physical 
properties of the disturbances as different theories propose different observable 
properties. The STEREO spacecraft also offered the chance to identify previously 
unknown characteristics and use these as an additional discrimination between 
different interpretations. 

To this end, the project has been successful. Whereas original observations 
of CBF pulses returned kinematics estimates that were inconsistent with a fast- 
mode wave, it has been shown here that these underestimations may have been a 
consequence of the low observing cadence of SOHO/ETT. It has also been shown 
that CBF pulses display clear deceleration and dispersion with propagation, both 
of which are wave properties consistent with the propagation of a pulse through 
a randomly structured medium. 

The observations of dispersion also allowed the physical nature of the pulse 
to be examined. Dispersion in a single Gaussian pulse can be studied by treating 
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the pulse as a linear superposition of sinusoidal waves, allowing the dispersion 
relation of the pulse to be determined. This has important implications for the 
analysis of both CBFs and also seismology of the solar corona as it allows the 
variation in the morphology and kinematics of individual CBFs to be used to 
directly probe the structure of the solar corona. It also allows the basic physics 
of the pulse to be examined theoretically, potentially providing a way to improve 
theoretical modeling of CBFs. 

CBFs have also been shown to be present in the He II 304 A emission line 
for the first time. While this was originally considered to be due to the coronal 
Si XI emission line at ~1.6 MK, it has been suggested that this emission line is 
too weak in the quiet Sun to be a viable candidate. The pulse signature implies 
that the CBF may not be a purely coronal feature, although more recent results 
indicate that the temperature gradient resulting from the propagation of a coronal 
pulse would be expected to produce an increase in He II emission similar to that 
observed. The CBF signature observed in the 304 A passband can therefore be 
taken as an indication that the CBF is a coronal pulse which produces a heating 
effect. This is consistent with the observations of passband-dependent kinematics 
made using SDO. 

8.1.1 Wave Kinematics and Morphology 

This work has shown that the kinematics of CBFs are more complex than orig- 
inally perceived. While they were originally modeled as constant velocity single 
pulses, it has been shown here that CBFs do exhibit deceleration with prop- 
agation, a result consistent with the previous observations by Warmuth et al. 
(2004a,b); Warmuth & Mann (2005); Vrsnak et al. (2006) and Veronig et al. 
(2008). This deceleration is also passband dependent, indicating that the CBF 
is strongly influenced by the plasma through which it passes. This may help to 
explain the wide variation of CBF kinematics noted by multiple authors including 
Thompson & Myers (2009) and Wills-Davey et al. (2007). 

As well as the clear deceleration shown by the pulses studied here, it was 
also shown that the observing cadence used to study the different CBFs strongly 
influences the derived kinematics. The very low pulse velocities determined using 
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data from SOHO /EIT may have been a direct consequence of its low observing 
cadence. Although these low kinematics were a contributing factor in the number 
of theories proposed to interpret CBFs, the results presented here indicate that 
they are not inconsistent with the fast-mode wave model originally proposed to 
explain the Moreton and "EIT wave" phenomenon. The effects of the low observ- 
ing cadence are particularly evident for the event studied in Chapter 7, where the 
same event exhibits different kinematics when observed using STEREO /EUVI 
and SDO/klk. 

The estimation of the CBF kinematics was also greatly assisted by the use of 
a statistically rigorous bootstrapping technique. This represents a huge improve- 
ment on the numerical differentiation techniques typically used to determine the 
kinematics of CBF pulses. It has been shown that the inherent nature of these 
techniques make them inappropriate for use with the small data sets typically 
derived from CBF observations. 

The physical properties of CBFs have also been determined to a high degree 
of accuracy. It has been shown here for multiple events that CBF pulses show 
broadening with propagation. While this had been previously shown in obser- 
vations made by other authors (e.g., Warmuth et al, 2004b; Warmuth & Mann, 
2005; Muhr et al, 2011; Veronig et al, 2010; Warmuth et al, 2001), the rigor- 
ous technique used here to identify the CBF pulse has shown conclusively that 
this is an inherent property of the disturbances. This has major implications for 
the study of CBFs as it implies that the theoretical soliton model proposed by 
Wills-Davey et al (2007) in particular is not appropriate for describing CBFs. 

The deceleration and dispersion of the CBF pulses also raises issues for the 
pseudo-wave interpretations proposed by Chen et al. (2002), Attrill et al (2007a) 
and Delannee et al (2008). None of these proposed models are consistent with the 
observations of pulse reflection observed in the 2007 May 19 event (Gopalswamy 
et al, 2009), while neither deceleration or dispersion is predicted by the authors 
of these theories. The passband dependence of the pulse kinematics is at odds 
with the low coronal nature of the pulse predicted by Attrill et al (2007a), while 
only single pulses were observed in each case studied, with no evidence of the 
coronal Moreton wave predicted by Chen et al (2002). 

Thus, basic physical principles have been used to show that the best represen- 
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tation of the CBF pulse is as a linear superposition of sinusoidal waves within a 
Gaussian wavepacket propagating in a randomly structured medium. This inter- 
pretation is driven by the various observations presented here, and is consistent 
with the observed deceleration and dispersion. This interpretation has also been 
used to quantify the dispersion relation of the pulse and the effect of the pulse 
on the local medium (and vice versa). 

8.1.2 CBF Identification 

Although the kinematics and morphology of CBF pulses have important impli- 
cations for our understanding of the solar corona, these results are predicated on 
an unbiased approach to the identification and tracking of "EIT waves". This 
phenomenon is inherently difficult to identify and track in individual images as a 
direct consequence of the diffuse nature of the disturbance. The result of this has 
been to complicate analysis of these disturbances and produce the wide variation 
in theoretical interpretations discussed in Section 2. 

Whereas previous work has relied on point-and-click identification of a CBF 
in running difference images, it has been suggested that this produces an inaccu- 
rate estimate of the pulse location (Attrill, 2010). This approach is also entirely 
user-dependent, with the scaling required to highlight the pulse and the actual 
identification of the pulse itself highly variable and unquantifiable. However, it 
has been suggested by Muhr et al. (2011) that the kinematics obtained using 
both an intensity profile approach and a visual tracking technique are compara- 
ble. In this case, the front of the pulse was used in both cases rather than the 
pulse centroid, with the pulse positions obtained from both approaches matching 
within error. However, this approach can mask the presence of deceleration in a 
dispersive pulse and give rise to inaccurate pulse properties (cf. Section 6.4.2). 
The use of running-difference images can also complicate pulse identification if 
the pulse position overlaps between images. This issue was overcome in this work 
through the use of percentage base difference images, which remove the coronal 
background to highlight the pulse using a nominally inherent scaling. 

These images were combined with an unbiased, semi-automated technique 
which was proposed for the identification and tracking of CBFs. This allowed 
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the pulses to be identified in a consistent manner, and enabled their deceleration 
and dispersion to be quantified. The technique was subsequently applied to data 
from both the STEREO and SDO spacecraft, exhibiting an ability to return 
accurate pulse kinematics and morphological variation regardless of spacecraft or 
event. However, it should be noted that there are issues with the intensity profile 
technique as currently implemented. The orientation of the arc used to determine 
the intensity profile along with the source point used to define the crossing point 
of the arc boundary are defined by the user. The detection of the pulse can also 
be influenced by stationary bright features as noted in Section 7.3.2. As discussed 
in Section 8.2.4, it is intended to improve the detection algorithm to negate these 
issues. 



8.2 Outstanding Questions and Future Work 

Although this work has provided many answers to the "EIT wave" question, there 
remain many unresolved issues with these disturbances. CBFs remain rare phe- 
nomena, even by comparison with solar flares and coronal mass ejections. Their 
exact relationship with these other solar phenomena also remains uncertain, al- 
though observations from SDO may help to resolve this issue (e.g., Patsourakos 
et al, 2010). From observation, CBFs are strongly associated with specific ac- 
tive regions, with AR 10956 producing ~6 CBFs over the course of ~5 days in 
May 2007 while similar active regions produce no eruptions. Similarly, Kienreich 
et al. (2011) reported observations of four CBF pulses erupting from AR 11067 
over the course of an 8 hour period from 28-29 April 2010. The reasons for this 
selectivity remain unknown, but could potentially have a major role to play in 
predicting space weather. 

8.2.1 What Causes "EIT Waves"? 

Although the evidence continues to mount that CBFs are waves as originally 
suggested, the exact nature of their initiation remains enigmatic. This is partly 
a consequence of the previously low temporal cadence of the various observing 
instruments but also arises from the rare nature of CBFs. The strong relationship 
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of CBFs to CMEs suggests that it is most likely the initial driver, although the 
low rate of occurrence suggests that this is not a simple relationship. As with the 
association of CBFs with active regions, there may be a number of criteria that 
must be met before a CME eruption from an active region produces a CBF. 

These criteria can only be identified through the use of a large sample of CBF 
eruptions. While the analysis of a single eruption may suggest some criteria, 
these may be negated by the analysis of an alternative eruption event. This 
effect can be seen in the previous work done on this topic, where theories have 
been proposed based on single events. A thorough analysis of CBFs must be 
performed on a large sample of events using an automated approach to negate 
individual bias. 

8.2.2 How do "EIT Waves" Affect the Corona? 

The effect of the passage of a CBF on the background corona is most apparent 
in the work outlined in Chapter 7, where the plasma variation produces differ- 
ent derived kinematics for each passband studied. This effect has never been 
previously noted and may provide a link between the coronal "EIT wave" and 
the chromospheric Moreton wave. This has important implications for the study 
of CBFs and also for the study of the solar corona through the use of coronal 
seismology. By quantifying the effect of the plasma on the pulse parameters it 
should be possible to provide a way to study the corona directly using CBFs. 

The different passbands observed by the AIA instrument onboard SDO also 
allow the temperature and density variation of the solar corona and a CBF to be 
studied. By examining these parameters before and after the passage of a CBF, it 
should be possible to identify the processes involved in the propagation of the CBF 
pulse, as well as providing a way to determine the validity of different theoretical 
interpretations. The launch of SDO makes this an area of rich potential over the 
next few years. 

Although SDO provides seven EUV passbands with very high cadence and 
spatial resolution, the true variation in temperature and plasma properties with 
the passage of a CBF pulse can only be studied to a high degree of accuracy 
using spectroscopy due to the broadband nature of EUV imagers. As a result, 
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instruments such as the Extreme ultraviolet Imaging Spectrometer (EIS) onboard 
the Hinode spacecraft will play a vital role in the determination of the basic 
physics at play in the corona and in a CBF pulse. 

8.2.3 How are "EIT Waves" Related to Other Propagat- 
ing Disturbances? 

The historical questions regarding the relationship between the coronal "EIT 
wave" and the chromospheric Moreton wave still remain unanswered. Although 
Moreton wave kinematics are typically estimated to be much larger than those 
of an "EIT wave" , the results in Chapter 7 suggest that this could be influenced 
by plasma variations. However, observations of Moreton waves remain rare and 
as a result the relationship between the two phenomena remains elusive. 

The SDO spacecraft provides an opportunity to observe CBFs with a cadence 
comparable to that used for studying Moreton waves. As a result, a well-observed 
CBF/Moreton wave event should provide a unique insight into the processes 
involved in this phenomenon. 

8.2.4 Automated CBF Detection 

The study of CBFs and the different open questions raised in the previous sections 
is entirely dependent on a rigorous technique for identifying and tracking CBFs. 
Quantifying the basic physical principles at play in a CBF pulse requires removing 
all associated error from the different measurements (where possible). Although 
there are inherent errors associated with the identification of features in images, 
the larger, unquantifiable human errors can be removed through the use of an 
automated approach similar to that proposed here. 

Although rigorous, the algorithm proposed here is only semi-automated and 
does have room for improvement. In particular, the identification of the source 
point and the orientation of the arc sector are currently user-defined. This needs 
to be fully automated, with the source taken as the location of the associated flare 
or some appropriate alternative, and the direction of the arc modified to include 
analysis of the entire Sun from that point. The sheer size of the SDO data also 
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needs to be accounted for to streamline the algorithm and speed up processing. 

In addition to the scientific level analysis, this algorithm could be modified 
to identify and track CBFs as part of a pipeline. This would ensure a consistent, 
bias-free CBF identification process and would allow the SDO data archive (along 
with those from SOHO and STEREO) to be rigorously examined and processed. 
Answering the outstanding issues in this field requires a consistent technique that 
can be provided using this approach. 
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